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COUNTING SLOW ELECTRONS WITH AN 
ELECTRON MULTIPLIER' 


By J. A. Cowan? 


ABSTRACT 


The efficiency of counting slow velocity electrons by a 12-stage activated 
beryllium-copper electron multiplier was investigated. The total electron current 
in the incident beam was measured using an electrometer tube circuit. The 
number of pulses produced in the multiplier by the whole beam was counted 
using high speed scaling circuits. The efficiency of counting was thus determined 
without depending on calibrated slits to divide the beam. 

Efficiency curves for primary energies between 200 and 900 ev. are shown for 
freshly activated surfaces and also for surfaces which had been exposed to still 
air for five months. The curves were similar in shape, the former being a few per 
cent above the latter. A maximum of 63% was obtained for 300 ev. input energies. 


INTRODUCTION 
The electron multiplier is used frequently as a particle detector. One great 
advantage is that it operates in a high vacuum and therefore requires no 


window to separate it from processes occurring in a vacuum. Also, its recovery 
time is so fast that the over-all counting speed is limited only by the associated 
circuits. 

Quantitative work requires that the counting efficiency be known. Very 
little work has been published on this aspect of the electron multiplier. For 
protons with energies of 10 kev., Robson (4) found 100% efficiency with an 
accuracy of 20%. In this case the current of incident protons was measured 
with an electrometer circuit and a calibrated slit permitted the protons to 
enter the multiplier. Also using calibration slits, Allen (1) gave the counting 
efficiency for electrons as 100% for 500 ev. and 40% for 6000 ev. input energies. 
Allen's results differ appreciably from those in this paper. The most probable 
explanation for this discrepancy is that the calibration of his slits was done at 
much higher currents than those actually used in the efficiency measurement. 
The curve shown goes over 100% at 300 ev. so the calibration may have 
changed considerably with the current density. 

Since the usefulness of the electron multiplier as a counter depends on its 
efficiency for detecting the incoming particles, this work was undertaken to 

1 Manuscript received in original form January 28, 1952, and, as revised, October 13, 1953. 
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investigate the possibility of using an electron multiplier for quantitative 
measurements of low velocity electrons which could if necessary be accelerated 
up to a few hundred volts within the detector. 
DESCRIPTION OF APPARATUS 

Electron Multiplier 

The electron multiplier was of similar design to that used by Allen (1), 
that is, a 12-stage multiplier with electrostatic focusing. The multiplier plates 
were made of 0.010 in. sheet beryllium copper. Nickel wires were spot-welded 
to the edges of the plates which fitted into holes in the supporting mica to give 
accurate positioning. Details of the spacing and curvature of these plates are 
shown in Fig. 1. The mica was held by screws to brass plates which in turn 
were bolted to the removable cover of the vacuum container. After being 
cleaned, the plates were activated by heating to a dull red heat in air at a 
pressure of 10~* mm. Hg for 20 min. 


MICA PLATE SUPPORT 






PLATE ELECTRODES 





ELECTRON COLLECTOR MULTIPLIER PLATE— AXIS OF ELECTRON GUN 


ELECTRODES Be Cu SHEET.OIO" 


Fic. 1. Cross section of the electron multiplier. 


The voltage divider which supplied potential to the plates consisted of 
12 2.2-megohm (1 w., 1%) resistors. Since Welwyn resistors were found to 
have excellent degassing properties, this voltage divider was placed inside the 
vacuum container. Thus only one insulated seal was required and the resistors 
were shielded. 

The output pulses from the electron multiplier are due to a number of 
electrons arriving at its collector. To make these pulses a maximum the 
effective output capacity was kept to a minimum in the following ways: 
(1) the collector was made as small as possible; (2) leads were kept short; 
(3) a cathode follower was used as the initial tube in the following circuits. 
The capacity of the collector, including wiring and the input capacity of the 
first tube, was measured to be 23 wyuf. The time constant of the multiplier load 
and all stray capacities in the apparatus used was 0.25 wsec. Making allowance 
for a random distribution of pulses with time, this allows counting at the rate 
of 10° counts per second with a statistical error due to unresolved counts of 
24% (8). 

To obtain the optimum value of the over-all voltage applied to the multiplier, 
a curve was taken of the number of counts, for a constant input, as the H.T. 
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voltage was varied. A value of 5000 v., 417 v. per plate, was chosen as it was 
in the middle of a broad maximum on this curve. The discrepancy between 
this value and the maximum of the efficiency curve (see results) is probably 
due to the fact that the best focusing of the multiplier plates occurs at 417 v. 


Electron Gun 

An electron gun was used as the source of primary electrons for testing 
the multiplier. Currents of the order of 107 amp. were obtained by using 
small limiting apertures in the grid and also by using a temperature limited 
filament. The geometry of the system was such that any electron going 
through both apertures would strike the first plate of the multiplier. The ratio 
of the grid accelerating voltage and the anode voltage was always set at such 
a value that the electrons converged as they approached the multiplier. The 
structure of this gun is shown in Fig. 2. 










GRID ELECTRON GUN 





cre" | WAX SEAL 






WAX SEALS 






OLLECTOR 
ASSEMBLY 


Fic. 2. Diagram of the electron gun, multiplier, and collector. 


The filament was a U-shaped 0.005 in. diameter tungsten wire, requiring 
about four watts input power. Since the gun was temperature limited it was 
necessary to stabilize the heating voltage supplying it. A 200 amp.-hr. storage 
battery was used in parallel with a selenium rectifier charging circuit. The 
charging current was kept within 0.1 amp. of the filament current to assure 
constancy of voltage over a period of an hour. This voltage changed by less 
than 0.0005 v. as measured on a potentiometer. 


Collector 
The collector was designed to pick up the total beam from the electron gun. 
It was a brass cup threaded onto the end of the moving rod of a sliding vacuum 


‘seal which allowed the rod to be pushed in or out. The side arm (see Fig. 2) 


in the glass section was made large enough to house the collector so it could be 
drawn completely out of the electron beam. The glass surfaces were coated 
with dri-film to ensure against electrical leakage, since the resistance to ground 
had to be kept high compared with the electrometer grid resistance (10! ohms). 
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Electronic Circuits 


(a) Amplifier 

The pulses at the output of the multiplier were of the order of 107% v. so it 
was necessary to amplify them in order to trip the scaling circuits. A wide 
band amplifier, 10 kc. to 8 Mc./s., was used to preserve the pulse shape. For 
this frequency response the rise time of 10% to 90% of a step function is 
1/20 usec. as stated by Valley and Wallman (5). Since this was much less than 
the resolving time of the input scaler, the band width was sufficient for this 
work. 

A standard resistance condenser coupled amplifier with peaking coils to 
extend the high frequency response was used. Its total gain was 1.5 X 10°. 
A multiposition switch on a low impedance attenuator was used as a gain 
control to assure that the same setting could be obtained repeatedly. A curve 
of the relative number of counts plotted against the amplifier gain is shown 
in Fig. 3. A gain of 5000 was used for the efficiency measure as the curve 
approaches the horizontal at that point. Above this value of gain, noise caused 
the number of counts to increase. A cascode cathode follower was used as the 
output tube to feed a negative pulse to the 200 ohm input impedance of the 


first scaler. 
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Fic. 3. Graph showing the variation of counting rate with H.T. voltage across the electron 
multiplier. 


(b) Scaling Circuits 

The first scale-of-four, which required a resolving time of } usec. to match 
the output of the multiplier, was designed following the general scheme shown 
in Vol. 19 of the Radiation Laboratory Series (2). 6AIK5 tubes were used for 
the Eccles-Jordon relays. To speed up the circuit, the anode voltages were 
limited at 6743 and 90 v. The circuit used counted regular pulses up to a fre- 
quency of 6 Mc./s. Up to this speed no difference was found between 6AL5 
tubes and IN35 crystals when used as holding diodes. 

The scale-of-four was followed by a scale-of-eight with a resolving time of 
2 usec. This was followed by a commercial scale-of-1000 and mechanical 
counter with a resolving time of 10 usec. This combination of scaling circuits 
was capable of counting random pulses at the rate of 10° per second with a 
statistical error of 23% due to unresolved pulses. 
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(c) Electrometer Circuit 

The current measurement was made by a balanced electrometer circuit 
which is shown in Fig. 4. The variable resistors were used to balance the 
circuit, so that small changes in battery voltages would cause negligible 
current flow through the galvanometer in the plate circuit. The control in the 
filament circuit had the greatest effect on the drift in the value of current 
through the galvanometer. Other variations due to lack of temperature 
equilibrium were eliminated by leaving the circuits connected continuously. 





INPUT 
FROM 
COLLECTOR 

100K 3V 


Fic. 4. Schematic diagram of the electrometer circuit. 


The value of the effective grid resistance, R, of the electrometer tube enters 
into the calculation of current. This resistance consists not only of the actual 
resistor but also of any leakage paths and the effect of grid current within the 
tube. 

The resistance measurement was made by determining the time constant 
of discharge of a condenser through R. To obtain a suitable time constant so 
that the time could be determined with sufficient accuracy with a stop watch, 
i.e. order of four minutes, the grid capacity had to be increased by adding an 
external condenser. Since R is of the order of 10'* ohms the leakage resistance 
of the condenser must be considered. Thus two variable standard air con- 
densers were used first separately and then in parallel. The discharge curves 
obtained showed no effect indicating a change in R throughout the voltage 
range used in this experiment. Several trials were taken for each case and a 
least squares solution applied to the results. R was estimated as 1.07 X 10” 
+ 1% ohms. 

The electrometer circuit was calibrated to find its deflection sensitivity. 
This was done by noting the deflection of the galvanometer in the plate 
circuit for a small known voltage added to the control grid of V1 in Fig. 4. 
With a grid resistance of 1.07 X 10” ohms and a galvanometer with a sensi- 
tivity of 107'° amp. per mm. the sensitivity of the circuit was found to be 
1.9 X 1077 amp. per mm. 

EFFICIENT MEASUREMENT 

The efficiency measurement was done by first measuring the electron beam 
as a current with an electrometer circuit and then counting the number of 
pulses produced by the electron multiplier. The whole of the beam was meas- 
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ured in both cases to eliminate the possibility of the transmission through 
calibrated slits varying with current density. The size of the current which 
could be used was governed by the possible accuracies of the two measure- 
ments. The larger the input currents the easier it would be to obtain an 
accurate current reading. On the other hand there is a limit to the speed of 
counting. A compromise of 107“ amp. or 60,000 electrons per second was 
chosen since currents are readily measurable down to 10~'* amp. and counting 
is almost impossible with present circuits much above 10° counts per second. 
By comparing the number of output pulses per second with the number of 
electrons per second equivalent to the input current, the efficiency could be 
found. Efficiency readings were taken for input electron energies ranging from 
200 to 900 ev. Below 200 ev. the earth’s field and other stray fields affect the 
electron beam so readings were not taken in that region. Since the maximum 
efficiency occurs near 300 ev. and the efficiency decreases gradually for higher 
input energies, readings were not taken above 900 ev. 


(a) Current Measurement 

With all voltages on the electron gun, the collector was placed so the electron 
beam would enter the middle of the cup. This was checked by moving in the 
collector while the electron multiplier was operating. The counting rate of 
the electron multiplier behind the cup was reduced to noise level when the 
collector was in the correct position, which was not critical. To ensure that no 
secondary electrons would be lost the collector was made 45 v. positive to 
the anode. 

The gun filament was turned off and the circuit was left for several hours so 
the electric charges on the glass envelope would come to equilibrium. Then 
the zero of the galvanometer was found by reading its deflection every 15 sec. 
for 20 min. The large number of readings allowed the evaluation of drift and 
also the averaging out of noise fluctuations. 

Next the filament of the electron gun was switched on to give a current of 
about 10~“ amp. At the same time the galvanometer deflection was kept near 
to the zero reading by adjusting the bias voltage in the grid circuit. When 
the deflection became steady the galvanometer deflections were read every 
15 sec. for 20 min. as before. 

The additional bias voltage almost balances the voltage caused by the 
current of the collected electrons flowing through the grid resistance (1.07 
x 10" ohms). The two sets of deflection readings were plotted against time 
and the mean curves drawn. By extrapolating the curve from the first set of 
readings to allow for the steady drift, the change in deflection of the galvano- 
meter was found at the end of the second set of readings. This deflection was 
multiplied by the sensitivity of the circuit and added to the current measure- 
ment (from the bias setting) as a small correction. 


(b) Counting Measurement 
Immediately after the electrometer measurement the collector was with- 
drawn from the electron beam and the 45 v. bias removed from it. The multi- 
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plier voltages were switched on and the reading of the mechanical counter 
taken every minute for 20 min. 

The number of counts per minute was plotted against time and the curve 
extrapolated back to the time at which the difference in electrometer de- 
flection was taken in part (a) above. The statistical correction (3) for the 
number of counts lost as a result of the } usec. resolving time was added to 
the counting reading. This correction was usually less than 23%. 


From these two measurements, (a) and (0), the efficiency can be found by 
dividing the value of the counting measurement by the value of the current 


measurement in the same units. 
An example of an efficiency measurement is shown in Fig. 5. A and B show 
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Fic. 5. Sample of the current and counting readings. 
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the readings of the galvanometer deflections. The large variation is due mainly 
to the thermal noise in the effective grid resistance of the electrometer tube. 
Averaging many readings tends to minimize the effect of this noise. Fig. 5C 
shows the counting readings. The slope of the current and counting curves 
is due to the gradual decrease of the electron emission of the electron gun. 
During a measurement, the current and counting readings decreased similarly. 


RESULTS 
The electron multiplier efficiency curve for different primary energies is 
shown in Fig. 6. From 200 to 275 ev. the curve rises steeply ; from 275 to 350 ev. 
there is a broad maximum at approximately 63% efficiency; above 350 ev. 
the efficiency decreases gradually to 56% at 900 ev. Activated multiplier 
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Fic. 6. Graph showing the efficiency of counting low velocity electrons with a Be-Cu 
electron multiplier. 


plates which had been exposed to still air for five months gave a very similar 
curve. This would indicate that plates activated in the manner described were 
very stable towards exposure to still air. 

The maximum deviation from the efficiency curve is about 5% giving a 
probable error of about 2%. The measured efficiency for primary energies 
below 225 ev. may be too high. This is probably due to electrons not being 
picked up by the collector and still striking the first plate of the multiplier. 
This effect was definitely detected below 200 ev. 

As shown by the efficiency curves, the electron multiplier can be used 
directly for counting electrons with energies down to 250 ev. with efficiencies 
over 50%. The electron multiplier permits very fast counting but its main 
disadvantage is that the output pulses have a large range of amplitude. Since 
the electron multiplier operates in a vacuum and has its maximum efficiency 
at about 300 ev., it is evident that very slow electrons of the order of a few 
volts in energy may easily be accelerated, in the detector itself, to a speed 
permitting counting. 


ACKNOWLEDGMENTS 
The author wishes to express his thanks to Prof. J. M. Anderson and 
Prof. R. W. McKay who gave valuable advice and assistance at all times. 


nr eee En eta 








COWAN: COUNTING SLOW ELECTRONS 109 


During this investigation financial assistance was received from the National 
Research Council by studentships (1947-1950). 


REFERENCES 


1. ALLEN, J. S. Rev. Sci. Instr. 18: 739. 1947. 

2. CHANCE, B., HuGHEs, V., MacNicuot, E. F., Sayre, D., and Wittrams, F. C. Wave- 
forms. Vol. 19, M.I.T. Radiation Laboratory Series. McGraw-Hill Book Company, 
Inc., New York. 1948. p. 607. ; 

3. Lewis, W. B. Electrical counting. MacMillan & Co., Ltd., London. 1942. p. 116. 

4. Rosson, J. M. Rev. Sci. Instr. 19: 865. 1948. 

5. Vattey, G. E. and WaLLtMAN, H. Vacuum tube amplifiers. Vol. 18, M.1I.T. Radiation 


Laboratory Series. McGraw-Hill Book Company, Inc., New York. 1948.  p. 80. 








FINE STRUCTURE OF THE SCHUMANN-RUNGE BANDS NEAR 
THE CONVERGENCE LIMIT AND THE DISSOCIATION ENERGY 
OF THE OXYGEN MOLECULE! 


By P. Brix? AaNp G. HERZBERG 


ABSTRACT 


The Schumann-Runge absorption bands of O2(B*2y — X*Z7Z) have been 
photographed in the fourth order of a 3 m. vacuum spectrograph with a resolution 
of 160,000. Some spectra were taken with the Oz at liquid air temperature. A 
detailed fine structure analysis has been carried out for all bands with v’ > 11. 
In addition to the six main branches (with AJ = AN = +1), for low values of 
the quantum number JN (total angular momentum apart from spin), several 
lines of the six satellite branches (with AN = +1, AJ = 0, #1) as well as of 
the two ‘‘forbidden’’ branches (with AN = +3, AJ = +1) have been identified. 
Values of the rotational constants and the vibrational quanta in the upper state 
have been derived up to v’ = 21. The triplet splitting increases rapidly with 
N and with v’; it cannot be described accurately by the known theoretical 
formulae. 

The origin of the 21-0 band is at 57115 cm™. A very shortextrapolation gives the 
convergence limit at 57128 + 5 cm™. This limit agrees excellently with the one 
derived from the near ultraviolet (A°2% — X%Z7Z) bands if it isassumed that at 
both limits those O atoms that are produced in the *P state are in the lowest 
component level of this state, viz. *P2. A discrepancy pointed out earlier by 
Herzberg is thus removed. The convergence limit just mentioned and certain 
other data derived from the spectrum lead to very precise information about 
the dissociation energy of Os. Without any extrapolation the dissociation energy 
into normal atoms can be given as 41260 + 15 cm™ (or 5.1143 + 0.002 ev. or 
117.96 + 0.04 kcal./mole), which is 0.63% higher than the old value. 


A. INTRODUCTION 


The rotational structure of the Schumann-Runge absorption bands of O» 
has been analyzed by Curry and Herzberg (4) for low values of the vibrational 
quantum number v’ of the upper state (v’ = 1,2,...,7) and by Knauss and 
Ballard (12) for higher v’ values (v’ = 4, 5,..., 15). In emission the Schumann- 
Runge bands have been studied by Lochte-Holtgreven and Dieke (13) for 
values of v’ = 0, 1,2 and v” = 11,12,...,21. While these bands have long 
been known to be due to a *27 — *Z7Z transition (Mulliken (17)), no detailed 
analysis of the triplet structure has been made except that Lochte-Holtgreven 
and Dieke did observe a splitting of the lines of high rotational quantum 
number and Knauss and Ballard found a splitting of the lines of the 14-0 and 
15-0 bands. 

Slightly different values for the convergence limit of the Schumann-Runge 
bands have been obtained by different methods of extrapolation, the most 
recent ones being 56877 cm~! (Gaydon (7)) and 56859 cm~! (Herzberg (10)). 
Since the upper state B*27Z dissociates into one normal (*P:) and one excited 
(Dz) oxygen atom (Herzberg (9)), and since the term difference *P. — 'Dz is 
well known from the OI spectrum, the dissociation energy Do° of the O2 molecule 
in its ground state is immediately obtained from the convergence limit of the 

1 Manuscript received August 24, 1953. 
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Schumann-Runge bands. This is the value usually quoted in the literature. 
However, in a recent paper (Herzberg (11)) it was pointed out that a precise 
evaluation of the convergence limit of the near ultraviolet (A*2{ — X*Z>) 
bands of O, leads to a value for Do®(Oz) that is 220 cm~! higher than that 
obtained from the Schumann-Runge bands. This discrepancy can be explained 
either by the assumption that at the two convergence limits the O atoms are in 
different component states of the *P ground state or by the assumption that the 
extrapolation of the convergence limit of the Schumann-Runge bands is 
inaccurate. 

In order to settle this question and to obtain a definitive value for the dissocia- 
tion energy of oxygen, new measurements of the Schumann-Runge bands near 
the convergence limit were undertaken with the new 3 m. vacuum spectrograph 
of this laboratory. Since near the convergence limit the triplet fine structure of 
the bands becomes of the same order of magnitude as the rotational structure, 
a detailed study of this triplet structure was necessary in order to analyze the 
bands. The problem of the variation of the triplet splitting near the convergence 
limit is, in itself, one that deserves some interest and about which very little is 
known at present. ; 

In the present work, all the Schumann-Runge bands from v’ = 12 right up 
to the convergence limit were analyzed and a new value for Dy°(OQz) was ob- 
tained which we believe is definitive and accurate. This value differs appreciably 
from the previously accepted value but is close to that derived from the near 
ultraviolet bands. 

B. EXPERIMENTAL 

The vacuum spectrograph used for the study of the Schumann-Runge bands 
was of fairly conventional design. A schematic diagram of the apparatus is 
shown in Fig. 1. Slit S, plate holder P, and grating are all contained in a large 
pipe of 26 in. internal diameter. Plate holder and grating can be adjusted from 
outside for various spectral regions without breaking the vacuum. The plate 
holder which carries 2 X 18 in. plates can be moved up and down so that 
several spectra can be taken on the same plate. 





Fic. 1. Diagram of apparatus. 
A: absorption tube, C: cylindrical LiF prism, Fe: iron arc, K: collimator lens, L: Lyman 
discharge tube, P: plate holder, S: slit. The circular arrow indicates how the arm L,K,A is 
rotated for the iron arc exposure. 
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Several gratings were tried as they became available. A grating ruled by D. 
Richardson of the Bausch and Lomb Co. proved to be the best of these. It has 
30,000 lines per inch and a ruled area of 40 X 90 mm. The blaze is at about 
6000 A in first order and correspondingly shorter wave lengths in higher orders. 

In building this vacuum spectrograph one of the main objects was to obtain 
a high resolution by the use of high orders. In order to study spectra in a high 
order, particularly in absorption, it is necessary to eliminate overlapping lower 
orders. This was done by means of a cylindrical LiF lens (C in, Fig. 1) in the 
way described by Curry and Herzberg (4). This lens acts as a fore-prism and 
at the same time forms a line image of the source on the slit. By means of this 
cylindrical lens it has been possible to obtain absorption spectra down to 1200 A 
in the fourth, fifth, and sixth order with dispersions of 0.65, 0.5, 0.4 A/mm., 
respectively, and resolving powers up to 200,000. Because of the much greater 
intensity of longer wave length light, which is far better reflected by the grating, 
it is imperative for a successful operation of the fore-prism to exclude completely 
all light scattered by dust on the prism or small fractures in it. In order to focus 
different wave length regions on the slit S, the whole arm L, K, A, which was 
mounted on a triangular bar, was rotated about C, the position of the cylindrical 
prism. 

The absorption tube A (Fig. 1) of 50 cm. length was separated from the 
spectrograph by the cylindrical prism C and from the light source L by a plano- 
spherical lens K which produced approximately parallel light through the 
absorption tube. For some experiments an absorption tube which could be 
cooled to liquid nitrogen temperature was inserted at A. It consisted of a quartz 
tube of 27 cm. length with fused-on thin quartz windows. This tube was surround- 
ed by a glass jacket filled with liquid nitrogen and this in turn was surrounded 
by an evacuated brass tube. 

The light source was a Lyman discharge tube viewed end-on. A 2 uf. condens- 
er was charged to 16,000 volts and discharged through the tube which was 
lined with a ZrO liner of 3 mm. internal diameter. From 100 to 1000 flashes were 
required for an exposure depending on the slit width used, the type of plate, 
the spectral region, and other factors. After about 8000 to 10000 flashes the liner 
had to be renewed since its diameter had increased to 5 mm. and the spectrum 
became more and more spoiled by many emission lines. 

For most exposures Ilford Q2 plates were used; for a few final exposures 
Ilford Q1 plates were used which have a much finer grain but require much 
longer exposure times. 

In the present work all spectra were taken in the fourth order. Only the 
region from the head of the 11-0 band at 1804 A to the limit of discrete absorp- 
tion at 1750 A was studied in detail. The pressure of oxygen in the absorption 
tube was varied from about 3 to 350 mm. yielding from weak to nearly complete 
absorption. Fig. 2 gives a general view of the whole region studied. With 
increasing pressure many weak lines are brought out in the ‘‘windows”’ left by 
the strong absorption lines (see Fig. 3). Since the strong absorption lines at 
high pressure appear very broad, several spectra at different pressures had to be 
measured in order to include all lines at the best conditions for their measure- 


ment. 
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The spectra taken at liquid nitrogen temperature were mainly used for a 
comparison of intensities. It is seen, for example, from Fig. 3 that most of the 
lines that appear weakly at room temperature are absent at low temperature 
indicating that they originate from the high rotational levels of the ground state 
which are less populated at lower temperature. On the other hand, a few weak 
lines such as those marked 7R in Fig. 3 to the left of the 12-0 head are relatively 
enhanced at low temperature and must therefore originate in low rotational 
levels. 

Some lines were measured on the low temperature plates rather than on the 
room temperature plates because they are less blended. These lines are marked 
by L in the tables. 

All measurements were made against second order iron lines. These were 
obtained by sending the light from an iron arc (Fe in Fig. 1) through the source 
and the absorption tube in a position in which the region 3000-4000 A was 
concentrated by the cylindrical prism on the slit. The wave lengths of the Fe 
lines given in the M.I.T. Tables (8) were converted to vacuum by using Edlén’s 
tables (6). By dividing these vacuum wave lengths by 2, standards were obtained 
in the fourth order from which the wave lengths of the O, lines could be deter- 
mined. Judging from the agreement of measurements made on different plates, 
the relative wave numbers of unblended lines are accurate to about +0.05 cm™!, 
whereas the absolute wave numbers may be in error by +0.2 cm~!. Lines of 
equal and favorable absorption were resolved when they were more than 0.35 
cm~! apart, thus indicating an actual resolving power of approximately 160,000. 


C. ROTATIONAL ANALYSIS 


(1) Assignment of Branches 

Fig. 5 shows an energy level diagram for a *2y — *Z7 transition. The figure 
is drawn to scale for the 15-0 band of the Schumann-Runge system and shows 
at the bottom the resulting spectrum for low values of NV, the rotational quantum 
number corresponding to the total angular momentum apart from spin!. In a 
3 state, for each value of NV, except N = 0, there are three sublevels designated 
F,, F2, Fs; depending on whether the quantum number J of the total angular 
momentum including spin is V + 1, N, or N — 1, respectively. The selection 
rules applicable to a *> — *2 transition are 

AJ =0,2+1, AN=+1, ser-a, +<¢—-. 


All of these rules, except AN = + 1, apply rigorously. The two symmetry 
rules are automatically fulfilled since, on account of the zero nuclear spin of 
O'®, only symmetric levels (s) occur and since, therefore, in the upper state only 
‘negative’, in the lower state only “positive” rotational levels occur. On the 
basis of the selection rules for J and N, the twelve branches P, P», P3, Ri, Ro, 
R; and ? Qi, ’Qes, ’Ris, “Qei, “Qs2, “P31 result which are indicated by full and 
broken vertical lines in Fig. 5. According to well-known intensity rules (see 
Herzberg (10)) the branches with AJ # AN have rapidly decreasing intensities 


'Previously this quantum number was designated K. We follow here the recent recommendation 
of the Joint Commission for Spectroscopy (J. Opt. Soc. Amer. 43: 410. 1953). 
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Fic, 5. Energy level diagram of a *27 — 839 band. 

The figure is drawn to scale for the 15-0 band of the Schumann-Runge system. The vertical 
lines representing the transitions are drawn in such a way that their projections below give a 
representation of the spectrum. The main branches are indicated by full lines, the satellite 
branches by broken lines, and the ‘‘forbidden’’ branches (with AN = + 3) by dotted lines. 
The (nonexisting) N = 0, J = 0 levels and, correspondingly, the band origin vo are marked 
by dot-dash lines. Estimated intensities are indicated by different heights of the ‘‘lines’’ in 
the schematic spectrum. The lines 7Ru(3) and 7Ry(5) originating from higher rotational 
levels of the upper state could not be shown in the energy level diagram. Their positions have, 
however, been marked in the schematic spectrum below because they fall into the region 
of the band head. The lines ?Qi2(1) and ?Ri3(3) as well as the satellite lines for the R(3) tran- 
sition (not drawn) have not been observed in any of the bands under investigation. 
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with increasing J; they are satellite branches, while the branches with AJ = AN 
are the main branches with the regular intensity distribution characteristic of 
normal branches. In the present case in addition to the main branches, lines of 
all six satellite branches have been observed at low N values. 

The selection rule AN = + 1, while not rigorous, is a fairly strong rule for 
>y — & transitions. Exceptions to this rule have up to now been observed only 
in forbidden transition in which the normal branches are rigorously forbidden 
by the symmetry rules, e.g., in the *2{ — *2> bands (11) of Oo. Such transitions 
violating the AN = + 1 rule (but in accord with the other rigorous selection 
rules) may be produced by spin-orbit interaction. The present case appears 


‘to be the first one in which such transitions have been observed in an allowed 


band system. The transitions with AN = + 3 are of this type. They are marked 
by dotted lines in Fig. 5 and designated 7R3, and *P;; (or, for short, 7R and 
YP). It is easily seen from Fig. 5 that they are the only branches of this type 
possible in a *2; — *27 transition. They have been observed to occur weakly 
at low N values in all the Schumann-Runge bands investigated here and can 
be identified particularly well on the low temperature spectra (see Figs. 3 and 
4). The 7R branch gives rise to lines on the shortward side of the band head 
formed by the main R branches. The satellite and forbidden branches are 
mainly responsible for the crowding of lines near the band heads. 

The main P and R branches of the 12-0, 13-0, 14-0, 15-0 bands can be readily 
picked out by inspection. The assignments agree with those of Knauss and 
Ballard (12). The branches can be followed to high N values even when they 
pass through several bands of lower v’ values which obscure certain regions of 
the branches under consideration. In Fig. 3 some of the lines of high N (up to 
N” = 33) are marked. 

The identification of the satellite and forbidden branches is not quite as 
straightforward. However, as mentioned before, it was greatly aided by the 
low temperature spectrograms. On the other hand, it can be seen from Fig. 5 
for the 15-0 band how the exact wave numbers of the lines of the forbidden 
branches *P3 and 7R,3 allow one to ascertain which components of the main 
branches are P;, Ri, which P2, R2, and which P3;, R;. It turns out that in all 
unperturbed bands the order is with increasing wave number P;, P2, P; and 
Ri, Ro, R3. (Only R2(7), R3(7) and P2(9), P3(9) of the 19-0 band were found 
to be interchanged owing to a perturbation.) 

The triplet splitting is very small for the 12-0 band, and is only partially 
resolved for N < 17. It increases rapidly with increasing vibrational and 
rotational quantum numbers, e.g., for the R(29) “‘line’’ of the 14-0 band (see 
Fig. 3c) the total splitting has reached 9 cm™', and in the 19-0 band a total 
splitting of 18 cm~! occurs for the R(11) “‘line’’ (see Fig. 4a). 

For the 16-0 and higher bands, perturbations (cf. Fig. 4) involving extra 
lines and intensity anomalies occur. These perturbations in addition to the 
large fine structure splittings made the analysis somewhat difficult. The known 
energy levels of the ground state and a comparison of line intensities at the two 
different temperatures were used as the main clues. Once a strong perturbation 
was encountered it was usually not possible to follow the respective branch to 
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TABLE 16 


WAVE NUMBERS OF THE OBSERVED LINES CLOSE TO THE CONVERGENCE LIMIT 


This table lists the observed lines from the limit to the head of the 17-0 band (corresponding 
to the upper spectrogram in Fig. 4) except for the lines definitely assigned to the 18-0 and 19-0 
bands. For a complete list of wave numbers in this region, the lines of the 18-0 and 19-0 
bands have to be included beyond 57029.88 cm. Intensities are visual estimates from the 
plate reproduced in Fig. 4. The designations "P;,"R;, "PX 4, and "®X ,, will be explained in 
section D(1). 























Foagtem*) | Intensity | Assignment | Yyaccem 4) Intensity Assignment 
Sete eseRcag id ba onesaoseee aadelienstiianaaiment hominis ci assets aesigceceleaeiscnaeeiipssialemaace aptamer aa ae 
a 
57130.73 1 | PX;3(1) | 57051.94 3 

129.50 | 1 "Ri() 051.09 Ot 

128.65 2 mPX 3(1) 050.44 3 20R3(5) 
117.70 Ot I 048. 26 4 2Ro(5) 
116.82 2 "R;(3) 1] 047.12 | “4 21R3(7) 
116.72"! 2 046.84 it 

115.30 2 pee | ane 3 an 
114.36 | 3s se Patl) 045. 2 21R2 
112.98 | 2u, b | R21) 043.80 3 20 P2(5) 
112.25 | 1 | = *RQ32(1) | 042.67 | 3 1P3(7) 
108.73 0b 042.14 a 

103.68 | Ot 040.78 2 | 21P2(7) 
102.36 | 5 21R3(3) 038.83 2bb,d? 

101.58 1 036 .53 lu 

100.55 | 4 | R23) 035.21\ | 1b 

100.01 3 | 1 P3(3) 034.78 f 

099.00 | 1 | 034.14 2 

098.41 | 2 21P3(3) Hi 032.85 3bb,d? * 

095.69 | Ot,b? 030.91 3 * 

093. Ot 029.88 2 * 

093.30 | 2 | ™R5) | 026.48 1 

092.58 | 1 024.69 i 2 

091 2 022.26 1t,.b 

091. 2 "Pi(5) | 020.90 3 

089. | 3 019.10 ee 20R3(7) 

087. Ib 016.83 3 20R2(7) 

085.6 2 013.01 4 20P3(7) 

084. 1 010.89 4 20P2(7) 

083 | 3 009. 62 1 | 

082.92 | 1t,b 20ROo1(1) 009.14 1 | 

082.22 | 1 20R Ps(1) | 004.15 1 21R3(9) 

| : ( 

081 2 20Ro(1) | - 1 

080. 4 21R3(5). +2°RO32(1)|| 002.90 2 

079.59 1 | 001.62 fa 1R2(9) 

078.17 | 3 21R2(5) || 56998.90 1 

076.46 | 3 1P3(5) | 998.27 Y a3 |  ™Pa(9) 

076.16 | 0 996.77 Ag 

074.80 | 2 21P2(5) 996.30 1 21P2(9) 

073.00 | 4 HI 991.33 0 

072.00 3 |} 20R3(3) 990.65 | 0 

069.68 s |  0Ra(3) | 989.85 a 

067 .69 3b,d? 20 P3(3) 985.92 icey | 

066.82 | 3 20P2(3) 984. 23 0 

065.94 | 1 982.77 0 

063.91 3 i 979 .03 | Ot | 

058.94 | Obb® | 977 .96 a: |  20R3(9)° 

057.95 1 975.48 L233 |  20Ra(9) 
2 9 | : | | 

056.78 2 ies | 974.03 Ot 

055.93 1 P;(7) 973.45 ot 

055.24 1 972.23 | 5 

054.50 2 971.39 1 94@ 

053.42 2 970.26 | 3 20P5(9) 

052.56 0 H 967.82 | 10 20P2(9)* 
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TABLE Ib (Concluded) 
WAVE NUMBERS OF THE OBSERVED LINES CLOSE TO THE CONVERGENCE LIMIT 


























Yyae(cm™!) Intensity Assignment | Pyac(cm™?) Intensity Assignment 
ee x ees 

56965.80 2 || 56922.47 0 

964.52 1 917.61 10 20P,(11)* 

961.13 2 | 915.11 2 20P2(11) 

959.71 1 || 902.54 0 

958.09 1b? * | 901.31 1 

957.16 1 | 900.80 2 

953.97 2 898.93 2 * 

951.40 1d? 891.08 1b 

950.13 3b,d? * 889.74 0 

948.88 2b 888.64 2 

947.63 4 1R2(11)* | 885.02 1 

945.51 2 * 875.54 1 2 P2(13) 

943,22 4 1P3(11)* 872.37 2 

941.32 2 2P2(11)° 865.02 5 * 

938.77 1 | 863.35 1 

937.50 2 862.32 1 

936.32 1 861,17 0 

932.96 | oO 860.34 3 

931.08 0 858.74 6s? 

926.68 2 | OR) 854.85 2 | 0 3(13) 

924.22 8 20R2(11)* 852.22 2 | 29413) * 














| 
| 
| 
| 
| 
| 
| 


“Possibly several lines. 

Possibly one broad line. 

“Unresolved group of lines between the two adjacent lines. 
4Probably two lines of intensity 4 and 2. 


Abbreviations: b: Broad; bb: very broad. 
d: Double. 
s: Shaded to shorter wave lengths. 
u: Unsymmetrical. 
1: Difficult to measure, wave number has a greater uncertainty. 
*: Mainly or partly assigned to 18-0 or 19-0 band. 
°: Probably blended with weak line of 18-0 or 19-0 band. 


higher N values without ambiguity. The assignment of the lines in the 16-0, 
17-0, and 19-0 bands is partly shown in Fig. 4. 

Tables Ia and Id contain all the lines that have been assigned to branches. 
Between 1803 A and 1759 A approximately 100 of the weakest lines, out of a 
total of 500 measured ones, remain unclassified. The fraction of unclassified 
lines is much larger between 1759 A and the end of the discrete absorption 
spectrum at 1750 A. We have, therefore, listed all wave numbers measured in 
this region (Table Id). A discussion of the assignments "’X , and "*X ,, will 
be given in section D(1). 


(2) Term Values and Rotationai Constants 

For the lowest vibrational level of the X*2>7 ground state (v’’ = 0) the 
rotational levels are very well known (2,3, 1). We have not attempted to deter- 
mine them independently from the Schumann-Runge bands. Rather, the 
knowledge of the energy levels of the ground state was used to derive the energy 
levels of the upper state from the observed lines. The usual procedure of using 
combination differences or combination sums was found to be less practicable 
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in the present case since many of the branches are incomplete. However, 
sufficient unblended lines from different branches were usually available to 
derive the levels of the upper state from the known levels of the ground state. 

The following term values of the ground state were adopted: for the F2 
levels (J = N, cf. Fig. 5), the A:F./’ data given by Babcock and Herzberg 
were added up, using the mean observed values (Table IV, column 5 of Reference 
(2)) up to N = 25, and the computed ones (column 6) for N > 25. The dif- 
ferences F, — F, and F, — F; (measured by Burkhalter, Anderson, Smith, and 
Gordy (3) and Anderson, Johnson, and Gordy (1) up to NV = 25) were taken 
from (11), and for N > 25 from (2). Finally 2.87 cm~! were added to refer all 
term values to the (nonexisting) NV’ = 0, J’ = 0 level, which incidentally is 
1.09 cm—! above the lowest actual level V’”’ = 1, J’ = O (see Fig. 5). 

According to Schlapp (see subsection 3) the term values 7» of the F, compo- 
nents of the B *2z state should follow the simple formula: 


[1] T. — v) = F,(N) = .B,N(N + 1) — D,N2(N + 1)? 


where vp is the wave number of the band origin and B, and D, have their usual 
meanings. The constants vs, B,, and D, were derived from the ‘‘observed”’ 
term values by graphical methods. Only for the 13-0 band was equation [1] 
found to hold within the experimental error (up to N’ = 26). The data are 
incomplete for the 12-0 band. For v’ = 14 and v’ = 15 slight deviations from 
equation [1] occur, but inclusion of a term H,N*(N + 1) did not seem to be 
justified. For v’ > 15 large perturbations limit the accuracy of the rotational 
constants obtained in this way. 

For low v’ and N’, the Ps and P3; as well as the Rz and R; lines are not resolved. 

TABLE II 
ENERGY LEVELS FOR v’ = 21 

(Term values in cm~) 

















| | 
| | From Mean, ‘ 
+ ow J’ | R(N —1) P(N’ + 1) observed Calc. 
| : Auten: 
ae rs] - 57115.18 57115.18 | 
m3 2 | 57115.85t 115.66 115.66 57115.62 
® |. #1 =e 117.58 117.55 
a 4 | 117.80t 117.93 117.89 117.91 
F; 6 5 | 121. 26+ 121.20 121.20 
m: | © “oat aa 121.28 121.29 121.26 
ae oe 125.68 125.67 
Fy : 8 | 125.64 125.66 125.65 | 125.67 
Be lig 2 | ee (130.99) 131.56 | 
F, 10 130.98 131.01. | 130.99 | 130.99 
| 
Fz | 49 11 | 
AR. 22) 137.04 137.04 | 137.02 


From ®P3(1) : 57115.35t, from ®Q30(1) : 57115.12. 

oT, = [57114.77 + 0.159 NN’ + 1) — 0.000105 NN’ + 1)*]cem™. 

tA dagger indicates a larger uncertainty. 

()Term values in parentheses have been derived from blends in which the line in question is prob- 
ably completely hidden by another stronger line. 
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In order to obtain even in such cases approximate term values for the F; levels 
the following assumptions were made: (1) the calculated F2 levels are good 
approximations to the real F, levels; (2) both components (P2 and P; or Re 
and R;) have the same intensity. Let [F.] and [F3] be the term values derived 
from the wave numbers of the blends using the correct ground state levels. 
Then obviously 


[Fo] = F,—x, [Fs] = F;+x, 


and therefore 


[2] F; = [Fs] + ((F2] — F2) 


where for F, we can use the calculated value (i.e., calculated from the B’ and 
D’ values obtained from unblended lines). 

The final data are collected in Tables II, III, and IV. Table II compares, as 
an example, the term values obtained from different lines (listed in Table Id) 
for the important v’ = 21 vibrational level. The calculated F; levels are given 
in the last column. The band origins vo and the rotational constants B,’ and D,’ 
are collected in Table III. In this table, for the sake of completeness, the data 
given by Knauss and Ballard (12) as well as Curry and Herzberg (4) are 


TABLE III 


BAND ORIGINS, VIBRATIONAL QUANTA, AND ROTATIONAL CONSTANTS OF THE ScCHUMANN-RUNGE 
ABSORPTION BANDS 
































vo (cm™!) | AG'(v' +43) (cm) | B,’ (cm™) 
v | KB(12) | CH(4) | KB(12) CH(4) KB(12) CH(4) 10 en 
| Present Present Present | 
work work work 
| 

0 ‘ere 688.05" | 0.813% 

1 | 50045.65 | 665.1 0.798 

2 | 50710.7; 641.4; 0.785 

3 | 51352.2 617.64 0.770 4.38° 
4 | 51969.6 | 51969.8, | 593.0 591.8,° 0.754 

5 | 52562.6 | 52561.6,°| 561.7 561.15° | 0.735° 

6 | 53124.3 | 53122.8°| 534.2 534.0° 0.719 0.719° 

7 53658.5 | 538656.8° | 500.4 0.703 0.702° 

8 | 54158.9 | | 465.5 0.671 

9 | 54624.4 | | 428.9 0.651 

10 | 55053.3 | | 388.2 | 0.633 

11 | 55441.5 | | 345.4 | 0.593 | 

12 | 55786.9 | 55784.59| 300.7 300.92 0.554 | 0.562; 13 
13 | 56087.6 | 56085.51| 255.0 | 254.96 0.523 0.5247 16.8 
14 | 56342.6 ae 208.4 210.07 0.492 0.4836 21.2 
15 | 56551.0 | 56550.54 | 168.96 0.431 0.4399 25.7 
16 | 56719. 50 132.91 0.395; 34.3 
17 56852.41 | 102.13 0.347 | 45 
ig | | 56954. 54 75.64 | 0.296 152 
19 | 57030. 18 | 52.65. | 0.258 49 
20 | 57082. 83 | 31.94 0.207 76 
21 | 57114.77 | | 0.159 | 105 


"Obtained from the observed origin of the 1-0 band and the observed AG'(4). This is not an 
extrapolated value as erroneously stated by (12). 

>From Lochte-Holtgreven and Dieke (13). 

©From bands with v"’ = 1. 
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TABLE IV 
TERM VALUES OF THE TRIPLET COMPONENTS IN THE B*Z7 STATE 

The term values are given in the form F, — F2°*!*, etc., where F2°*!® are the calculated 
F; levels, using the constants vo, By, and Dy given in Table III. All data are in cm. The values 
for v’ = 21 are contained in Table II. 

Whenever the P2 and R: lines are blended with the corresponding P3 and R; lines, the F2 
levels derived from the wave numbers of these blends are given in square brackets []. The F; 
levels derived from these blends, corrected by the amount F2 — F2°*!°, as described in the 
text, are marked by a § sign. Values in parentheses () are interpolated ones. 

Extra levels connected with perturbations are given as F,’, etc. 
{Values marked by a dagger have a greater uncertainty. 
?A question mark indicates that the level is doubtful. 


ie 


ppp F.—Fo*™ p,-F 


-04 
3.73 


-1.6T | 
—1.84§ 
1. 
1. 
1. 
-1, 
-1, 
1. 
-1. 
1. 


[+0.04] 
[+0.07] 
[+0.10] 
[+0.15] 
{[+0.17] 
(+0. 
[+0.30] 
[+0.36] 
0.00 


2.56 
94 
2.75 
91 
05 
3.21 
.34 
.54 


22 





16 





= scale calc F; _ fale Fi —peale _ Peale wR 


F.—F2 F: 


5.57 ‘ E —9.02 





.17) 
[+0.14] 
[+0.26] 
{+0. 

+0. 
0. 
-0. 
—0. 
—0. 
+0 
(+0. 
+0. 
—0.06 
(—0.22) 
—0.38 


06 
00 
03 
04 
02 
03 
10) 
11 


©_1.71§. 
f40.15¢. 
9_4 188. 


37 


=~], 


—1.63" 


-1, 
-1, 


—1.2 
—0. 98 
—0.6: 
—0. 2: 





+0.12 
+0.6 


+1. 


+1. 


+2. 
+2. 
+2. 
+3. 


_ 


| —12. 


74 
94 | 
67 | 


4.23 
5.34 
-60 
.08 
.53 
.08 
.59T 
16 
9.71 
—10.31 
—10.80 
—11.30 
73 
21 
79 
47 


-9 


—12. 


—13. 


| —14.41 


| 


h 


Fe ~ 
: Fi’ -F 


j Ff- Fcale 


-1 
-1 
—1 
-1 


0.00 

[+0.14] 

[+0.13] 

—0.52h 

+0.15' 
—0.01 
—0.02 
—0.01 
0.00 
+0.02 
14 


—0 


+1 
+2 
+2 


56 +3 


37 


= +0.45f. 
= —1.29. 
= +6.81. 


. cale 
"2 


—1.% 


ped 
+0.3% 
+0. 


.30T 
.82§ 





—8.49 
—8.56 
—8.85 
—9.38 
—9.90 
—10.45 
—10.91 
—11.35 


0.00 
[+0.30]t 
[+0.18} 

—0.01 
+0.08t 

0.00 

+0.08 
+0.01 


-1. 


-1 
| 
-1} 
—-0 
-0 


+0. 


+1 


22 
85 
638" 
.30 
.78 
.10 
72 


64 
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TABLE IV (Concluded) 








calc , cale 
‘2 


|Fi-F F.—F2 Fo—Po 


—14.80t = ip 
—14.32 +0.01 —1.02 
—14.41 0.00? 1.712 
—14.75 [(-0.02] —2.01" 
—15.36 +0.05 —2.43° 
—16.01 —0.05 -+1.66° | 
—16.74 —0.02 +1.52 
—17.81 +0.01 

—18.60 


11 40.34). | Pr 40.16). 


 _1.33§. | 9% —1.69§. 
np? FS . ~16. 232. r 


Ope Cale 
Fi —F:* = 





—2.00§. 
nate 
Fi —F*"" = +3.48t. 


Fi? —F,** = —3.204. 


—10.41f?7, | 8 


t 








included. In Fig. 6, a plot of the B,’ values against v’ is shown. Instead of present- 
ing the actual term values we give in Table IV, in order to save space, the 
deviations of the observed term values from the calculated F» series, calculated 
from equation [1] with the final B, and D, values. The signs attached to some 
of the numbers of this table are explained in the heading of the table. 
By 
(cm-') 
0.8 


° 10 20 


Fic. 6. By,’ curve of the B%dzZ state of Ox. 


In Fig. 7 the observed energy levels above 56835 cm. are plotted against J. 
This figure gives an idea of the relative magnitude of the vibrational quanta, 
the rotational energies, and the triplet splittings near the convergence limit. 
It also shows the perturbations which occur in various places. 


(3) Triplet Splitting 
The observed triplet splittings F; — F,and F; — Fe: are readily obtained from 
Table LV by taking the differences of columns 1 and 2, and 3 and 2, respectively. 
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0 J 


Fic. 7. Observed energy levels of the B%ZZ state of O2 close to the convergence limit as a 
function of J. 

The F; components are represented by circles, the F, and F; components by crosses. Where 
perturbations and, associated with them, extra levels occur the lines connecting the levels of a 
given v’ are dashed. The small arrows at J = 0 point to the band origins vo (for the F2 levels); 
the larger arrow at the top shows the convergence limit (see section D). 


A graphical representation of the observed splittings in their dependence on 
N is given in Fig. 8 for v’ = 18, 15, 17, 19. It shows the rapid and accelerated 
increase of the splitting with increasing v’. The broken line curves give, for 
comparison, the triplet splitting of the ground state. Extrapolating the splittings 
of the upper state to smaller v’ values one sees readily that for v’ < 11 the 
splittings will be very similar in magnitude to that of the ground state and 
therefore, since the main branches correspond to AJ = AN, the triplet splitting 
of the lines in these bands is very small and has not been resolved in spite of 
the high resolution employed here. 

According to Schlapp (18) the three triplet components of a * state, neglect- 
ing centrifugal stretching terms, should follow the formulae 


[3] Fi(N) = B.N(N + 1) + (2N 4+ 3)B,— > : 


— V/(2N + 3)? B,? +d? — 20B, + y(N+)), 
[4] F.(N) = B,N(N + 1), 


[5] F(N) =B.N(N4+1)-(QN-1B,-v\ 
++/(2N — 1)?B2 +d? — 20B, — yN, 


where \ and y are constants: \ is a measure of the spin-spin interaction of the 
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TRIPLET SPLITTING 


(cm-) 
5 


visl9,F5 


D v'=I7, Fs 


aaa eeenecenenee = Fs 


"seas eeeneacan «= = F, 


v'=0 


OD v'si9,F, 


-20 
oO 10 20 30 N 


Fic. 8. Triplet splitting in the B*2Z state of O2 as a function of N for various v’ values. 

For comparison the triplet splittings in the ground state X #29 (v’’ = 0) are shown by broken 
lines. The splittings for v’ = 12, 14, 16, 18 have been omitted in order to avoid overcrowding 
of the diagram. These splittings are entirely similar to, and intermediate between, those shown. 
A strong perturbation is evident for the F; component of the v’ = 19 level. 


unpaired electrons, y is a measure of the interaction of the resultant electron 
spin with the magnetic moment produced by the rotation of the molecule. 
On the basis of Schlapp’s formulae the following sum rule should hold: 


[6] [F2(N) — Fi(N)] + [F2(N + 2) — Fs(N + 2)] = 244+ 7, 


that is, the sum at the left formed from the observed splittings should be inde- 
pendent of N. This is found to be the case for v’ = 13, 14, 15 within the experi- 
mental error (compare Fig. 9, see below). For v’ = 16, 17, 18, 19 the sums are 
no longer independent of NV. However, perturbations are present in all these 
bands and they may well be the only cause for the apparent breakdown of 
equation [6]. In Table V the average observed values of the sums, i.e., of 2A + vy, 


are listed. For v’ > 15 the values given refer to low N values only; for v’ = 12 
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TABLE V 
SPLITTING CONSTANTS IN THE B8ZZ STATE 





13 14 15 16 17 18 19 
2aA+y7 4.5 5.02 5.62 6.48 (8.0;) (10.2) (12.8) (15.8) cm™ 
- : . 
Values in parentheses refer to low N values only and are very uncertain. 


the partially resolved triplet splitting of the 12-0 band has been analyzed under 
the assumption that the sum rule is valid. In spite of these limitations the table 
shows clearly the rapid and smooth increase of 2A + y near the convergence 
limit, very similar to that found for the A*Z¢ state (11). Considerable effort 
was spent to find assignments for the components F; of the 20-0 and 21-0 
bands and thus determine 2\ + y for the v’ = 20 and 21 levels. These F, 
levels are certainly present near the energies expected by extrapolation. How- 
ever, owing to the presence of strong perturbations no definite assignments 
could be made. 

A further test of the validity of Schlapp’s formulae in the present case is 
obtained by ascertaining whether they represent the observed variation of 
the splitting with NV. From equations [3], [4], and [5], it follows immediately 
that 
[7] F,(N) — Fi(N) — 3(2\ + y) = Fi(N + 2) — Fo(W + 2) + 3(24 + 7) 


= — (2N + 3)B, +V/(2N + 3)°B,? + \*—20B, — y(N + 9), 


that is both F, — F; and F; — Fy» follow the same course if properly combined 
with (2A + y). In Fig. 9 the observed values of the left-hand side of [7] for 
the v’ = 15 level with (2A + y) = 6.48 cm™ are plotted against N. The 
validity of the sum rule is equivalent to the fact that the points for F; and F; 
follow the same curve. Since the quantity 2\ + y has already been determined 
and used (left-hand side of [7]) only one parameter is left to fit the observed 
points. The three broken line curves of Fig. 9 represent the right-hand side of 
equation [7] with y = 0, A = 3.24; y = — 0.18, A = 3.83; and y = 0.23, 
A = 3.35; cm! (i.e., in each case, 2A + y = 6.48). It is evident that the general 
character of the observed splittings can be represented by [7] with y = — 0.18 
cm, but no quantitative agreement can be brought about by varying y and A 
while keeping 2A + y at 6.48 cm~'. Similar discrepancies even though somewhat 
smaller are found for v’ = 14 and 13. Because of the perturbations no attempt 
was made to carry out similar comparisons for v’ > 15. 

Schlapp’s formulae have been found to hold within the accuracy of optical 
measurements for the X *Z7 ground state (2) and the A*Z} state (11) of Os. 
With the very high accuracy possible in the microwave region a slight deviation 
has been detected for the ground state by Gordy and his collaborators (3). 
Several attempts (14, 15) have been reported recently to improve the formulae 
in order to obtain an accurate representation of the microwave data. The 
modification of Schlapp’s formulae published by Miller and Townes (14) gives 
a slightly better fit for our data than the Schlapp formula (Fig. 9), but none of 
the formulae available at present gives a satisfactory agreement. This failure 
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y=-0.23 cm"! 7 
@ 


-1 
(cm~') yf 
7 O F2(N)-F,(N)-3.24cm"! ue 
“Oo 
X F3(N#2)-Fp(N+2) +3.24 cm~! Pog a 





Y=-0.18 cm~! 


° 10 20 30N 


Fic. 9. Test of Schlapp's formulae for the triplet splitting in the v’ = 15 level of the 
Bq state of Os. The observed values of the left-hand side of equation [7] have been plotted 
as circles (F:) and crosses ( F;). The three broken lines correspond to Schlapp’s formulae for 
various values of y. The full line has been calculated with Schlapp’s formulae as modified by 
Miller and Townes (y = — 0.18 cm™'). The experimental value 2A + y = 6.48 cm™ has been 
used for all curves. 


is probably connected with the fact that the coupling parameter y is very 
large in the present case, of the order of —0.18 cm~! for v’ = 15 of B*Zz as 
compared with only —0.00837 cm~'! forv’’ = 0 of X *Z7 and with —y < 0.016 for 
A* Xt. With the theoretical discussions still in progress we have not attempted 


to represent our data by empirical relations, e.g., by adding in Schlapp’s formulae 
correction terms to y which are proportional to N°. 


D. VIBRATIONAL ANALYSIS AND DISSOCIATION ENERGY 


(1) Convergence Limit 

By the rotational analysis described in the preceding section, six new vibration- 
al levels, with v’ = 16, 17,...,21, have been established beyond those pre- 
viously known. The last four of these, v’ = 18, 19, 20, 21, are higher than the 
-convergence limit that has until now been generally accepted.’ 

2Watanabe, Inn, and Zelikoff (20) have recently measured the absorption coefficients of oxygen 
in the Schumann region under low resolution. They noticed already that the discrete absorption 
extended beyond the previously accepted dissociation limit. However, the feature at 57140 + 60 cm 
is probably not the 20-0 band as they believed since, as may be noted from Fig. 4, none of the heads 


of the last four bands are conspicuous and there are much stronger accidental groups of lines which 
will give rise to peaks at low resolution. 
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In Table III the band origins v9 of all the absorption bands with v”’ = 0 
are given including the older measurements of Knauss and Ballard (12) and 
Curry and Herzberg (4). The fourth and fifth columns give the vibrational 
quanta AG'(v + 4) derived from them. These as well as the second differences 
A°’G’ are plotted in Fig. 10. 


4G’ 


0 
° 10 20 v' 


Fic. 10. AG’ and A*G’ curves for the B%Zz state of Ox. 

The observed AG’(v + 3) values are plotted with the abscissae v’ + }, the negative values 
of A?G'(v) = AG'(v + 3) — AG’(v — 4) with the abscissae v’. A magnified portion of the AG’ 
curve with the most probable extrapolation is shown by the light line in the upper part of the 
diagram (ordinate scale at right). 


Since the last vibrational quantum is as small as 32 cm and the last A?G’ 
values follow accurately a smooth curve (Fig. 10), only a very short extrapola- 
tion is necessary to obtain the convergence limit. For this purpose the maximum 
of the extrapolated G(v’) curve was determined where G, as usual, denotes the 
total vibrational energy. The continuous function G(v’) which we used was a 
polynomial of degree m exactly representing the last m + 1 vibrational levels 
(v’ = 21, 20, 19...,21 — m). In order to obtain the coefficients of the poly- 
nomial a difference table AG, A°G,..., A"G was constructed and extended to 
higher v’ values by keeping constant the last value in the A"G column. The 
maximum of G(v’) and the corresponding v'max, were found by interpolation. 
Table VI shows the results for various values of m > 1 and, at the same time, 
gives corresponding extrapolations for the upper state (A*Z{) of the near 
ultraviolet bands using the recent measurements by one of us (11). In both 
electronic states these extrapolations indicate the existence of one additional 
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TABLE VI 
CONVERGENCE LIMITS OF THE B82; AND A®Z{ STATES EXTRAPOLATED BY USING POLYNOMIALS 
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vibrational level. The ‘‘last AG’”’ and the difference 6G between the extrapolated 
limit and the extrapolated last level have been listed separately in Table VI. 

As long as a smooth extrapolation is assumed to be reasonable at all, the 
approximation should become better with increasing degree m of the polynomial. 
However the influence of the experimental errors of the input data on the 
differences A"G also increases with m. For this reason from m = 6 on these 
differences show an irregular variation rather than a definite trend with v’ 
for both the A*Z+ and B*>; states. The extrapolations with m = 6 were, 
therefore, made with a linear (least squares) representation of the A5G values. 
The limits so obtained are considered to be the most reliable ones of those given 
in Table VI. The tail of the AG’(v’ + 4) curve of the B*dz state in Fig. 10 is 
drawn according to the best numerical extrapolation just discussed. 

With increasing m the extrapolated limits of Table VI appear to approach a 
value which in the case of the B*2z state may be given to +1 cm~". This does 
not take account of the possible errors of the band origins (0.10 cm~') which 
have not been smoothed in the above-described extrapolation method and which 
may lead to a cumulative error of certainly less than +4 cm~'. We, therefore, 
conclude that the extrapolated convergence limit of the Schumann-Runge 
absorption bands can be safely given as 57127.5 +5 cm™!. 

According to the extrapolation of the B*2; state the v’ = 22 level is just below 
the limit and may have a few stable rotational levels. The origin vo of the 22-0 
band (F: component) is predicted at about 57127 cm~'. There is indeed at the 
very end of the discrete absorption spectrum near 57129.5 cm! a group of three 
absorption lines (see Table Ib and Fig. 4). Since all bands are strongly degraded 
to the red, these lines must be assigned to a vibrational quantum number 
v’ > 21. Assuming that no further lines are hidden by the adjoining continuum 
the lines of highest wave number must necessarily arise from combinations 
with NV” = 1 of the ground state. Combinations of the same levels with NV” = 3 
should then lie to the red by about 14 cm! (see Fig. 5). A gap of just this 
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magnitude actually exists between the last three lines and the next group at 
longer wave lengths. Furthermore, the lines at 57130.73 cm~! and 57128.65 
cm! show a separation of 2.08 cm~! which agrees exactly with the difference 
between the VN” = 1, J” = 0 and N” = 1, J” = 2 levels of the ground state. 
If this agreement is real, it would fix the higher one of the two levels with 
J’ = 1, that is, Fi(0) or F3(2) at 57129.64 cm—!. The line intensities favor the 
F, assignment. A P and an R branch can be picked out including all three lines 
mentioned. These branches can be followed up to V = 7, and they are designated 
"P, and "R; in Table Ib (or ""X x, and "*X x for the satellite branches) since 
neither their assignment to a specific vibrational quantum number nor to a 
specific triplet component is definitely established with the present data. 

Irrespective of the detailed identification of the lines of shortest wave length, 
it can be seen from Fig. 5 that an upper limit for the energy (v9) of the Fy» 
component of the highest stable vibrational level can be obtained by adding 
1 cm— (=F,'(1) — F./’(0)) to the wave number of the last discrete line. 
Hence, assuming again that no absorption line has been missed, the last vibra- 
tional level must lie below 57131.7 cm™, that is, the last vibrational quantum 
AG’ must be less than 16.9 cm~!. The energy difference between the convergence 
limit and the last vibrational level must be less than the last vibrational quantum 
and therefore an upper limit of 57148 cm~! can be given for the dissociation 
limit. This upper limit is independent of any extrapolation. 
(2) Dissociation Energy 

The. extrapolated convergence limits (see Table VI, m = 6) for the A* Xt 
state and the B*Z;7 state differ by 15868., cm~!. This difference agrees surpris- 
ingly well with the term difference 'D2 — *P, in atomic oxygen which is 15867.7 
cm~'(16). Thus there can be no doubt that the convergence limits under discussion 
represent real dissociation limits (see also below) and that upon dissociation 
the states of the resulting atoms are *P, + *P; for the A*Z+ state and 'D,. + *P; 
for the B* 7 state, with *P; the same for both states (the closest alternatives 
8P, + ®P, and 'Dz + *Po, or *P2 + *Po and 'D,. + *P), give an energy difference 
that deviates from the observed by +68 cm~! respectively). According to the 
non-crossing rule (see the discussion by Herzberg (11)) *P; can only be *P2. 
We shall from now on make this assumption. The dissociation energy Dy’ of 
the oxygen molecule into normal (*P:) atoms is then directly given by the 
dissociation limit of the A*2{ state or by the dissociation limit of the B*D7 
state less 15867.7 cm™. 

In Table VII the various data concerning D»® are summarized. The first three 
lines give the following information derived from the A* Z* state (11): 

(a) The last observed vibrational level at zero rotation gives a lower limit 
for Do. 


(b) The extrapolated convergence limit (Table VI) gives the dissociation 


energy directly. 

(c) There is evidence for predissociation by rotation in the v’ = 10, N’ = 17 
level and therefore an energy between that of the N’ = 15 and N’ = 17 levels 
must be an upper limit to the dissociation energy. 
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TABLE VII 
DATA REFERRING TO THE DISSOCIATION ENERGY OF THE O2 MOLECULE 





| : Dissociation limit (cm~') 
Obtained from —_——, 











3P2+'D. | 5P. + 5P, 
se ene Re Ua re eee s ~ = 
(a) | At vf = 10, N’=0 | >41153 
(6) Extrapolation = 41259 
(c) Predissociation (v’ = 10, N’ = 16) | <41284 
(dq) | Bs v' = 21, N’=0 | >57114.s | > 41247 
(e) | Highest level with J’ = 1 | 2 57129. | 241262 
(f) Extrapolation 57127.,45 | 41260 + 5 
(g) Last line and last AG’ <57148 | <41280 


Accepted value: 41260 + 5 cm= 
The following data arise from the present investigation of the B*27 state: 
(d) The energy of the well-established v’ = 21 level at zero rotation must be 

lower than the dissociation limit. 

(e) The energy of the highest rotational level with J’ = 1 as obtained above 
from the lines of shortest wave length represents a closer lower limit than (d). 
Conceivably it might even be very slightly above the limit if the possibility of a 
rotational potential maximum is considered. 

(f) The extrapolated convergence limit (Table VI), because of the short 
distance of extrapolation, may be expected to agree with the true dissociation 
limit within the inaccuracy of the extrapolation. 

(g) An upper limit for the dissociation limit was derived in subsection (1) 
from the observed line of shortest wave length. 

No evidence of predissociation by rotation was found in the Schumann-Runge 
bands. This is probably due to the presence of perturbations which made it 
impossible to predict the accurate positions of the high rotational levels for 
which such a predissociation could be expected. 

The data in the last column of Table VII (obtained for B*27 by subtracting 
15867.7 cm~! from those in the second last column) support each other very 
well, and they provide narrow limits for Do°®. However, the limits (a), (6); 
(d), (e), (f) are based upon the assumption that the potential curves of the two 
electronic states A*2Y*} and B*zz7 have no potential maxima. If there were 
potential maxima the extrapolations (6) and (f) would merely give upper 
limits while the data (a), (d), (e) would no longer give lower limits. 

It cannot be proved conclusively that there are no potential maxima present. 
However, the excellent agreement of the limits obtained from the A*Z}{ and the 
B*>z states can only be understood if either both states have no potential 


‘maxima or both states have potential maxima of the same height. It is difficult 


to believe that the latter coincidence would occur, particularly since the AG 
curves and, therefore, the potential curves of the two states are quite different 
(compare our Fig. 10 with Fig. 8 of Reference 11). In addition it appears that 
for both states the B,’ values extrapolate to zero at the convergence limits 
rather than to some finite value as would be the case if potential maxima existed. 
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This may be seen from a plot (not shown) of (B,’)" versus G(v’) where 7 is chosen 
in such a way that a straight line results (n = 7 for the A*Z} state, n = 4 
for the highest vibrational levels of the B*27 state). 

The preceding arguments do not exclude potential maxima that are only a 
few cm~' high. In order to take into account this uncertainty we shall assign a 
conservative error to our result. As the final value we shall adopt for the dissocia- 
tion energy of O2 (referred to N” = 0, J’ = 0; see section C(2)): 

D,°(Oz) = 41260 + 15 cm-!. 


This value does not depend on any extrapolation since the limit of error given 
includes a lower limit and almost includes a (conservative) upper limit. Using 
the conversion factors of DuMond and Cohen (5) the above value corresponds 
to 5.1143 + 0.002 ev. and 117.96 + 0.04 kcal./moleg.em (or with the older 
conversion factors used by Herzberg (10), 5.1133 ev. and 117.94 kcal./mole). 

The change in D(O.) affects of course all those dissociation energies which 
have been determined from thermochemical data involving oxygen. A partic- 
ularly important one is the heat of dissociation of NO. One finds from the above 
D(Qz) value, the dissociation energy of Ne (Gaydon (7), Douglas (4a)) and the 
heat of formation of NO at 0°K. (Rossini (19)) 

D,°(NO) = 52477 + 50 cm™! 


or, again using DuMond and Cohen’s conversion factors, D o°(NO) = 
6.5052 ev. = 150.03 kcal./mole. 

Note added in proof: Recently DuMond and Cohen (5a) have published a 
revised adjustment of atomic constants. Using the corresponding conversion 
factors increases the above values in ev. and kcal. by approximately 1 part in 
10,000. 


E. CONCLUSION 


The old value of the dissociation energy of Os was between 40990 and 41009 
cm~! depending on the extrapolation chosen. This extrapolation was considered 
to be accurate within the stated limits and in fact this case has often been quoted 
to illustrate safe and precise spectroscopic methods of determining dissociation 
energies. The new value for Do°(Oz) is higher by 260 cm~!. While this is not a 
large amount it is far outside the error considered possible previously. The 
discrepancy is of course due to the fact that the inflection of the AG’ curve and 
the “‘toe’’ of positive curvature (see Fig. 10) could not be detected as long as 
the vibrational levels were known only tov’ = 15. This example again emphasizes 
the need for caution in stating the results of such extrapolations. 

Like the old value the new Do°(O2) value is based on the assumptions (a) that 
on dissociation at the limit of the B*Sz state one O atom is in the *Pz, state 
(rather than *P,; or *Po) and (0) that no potential maximum occurs. The first 
assumption is in agreement with the non-crossing rule while any alternative 
assumption would contradict it. Arguments against the second assumption 
have been presented above. We therefore feel that the value for Do°(O2) here 
derived must be considered as definitive. 
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The analysis of the Schumann-Runge bands at the very end of the discrete 
absorption is still incomplete and a number of lines remain unclassified. As 
mentioned before our failure to assign these lines is due to the presence of per- 
turbations. Incidentally, the perturbations establish the presence of another 
stable electronic state presumably going to the same limit (*P,+ 'D2) as 
B*>;. However, with only a few “‘extra’’ levels available, no identification of 
the new electronic state is possible at present. 
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PHOTONEUTRONS IN A HEAVY WATER PILE! 


By M. W. Jouns? AND B. W. SARGENT? 


ABSTRACT 


Decay curves of the neutron intensity in the low power heavy water pile 
(ZEEP) at Chalk River as a function of time from shutdown with cadmium rods 
were obtained with three proportional counters of different sensitivities. The 
range of counting rates spanned, including operation at 50 w., was 8 X 10°. 
The durations of steady operation before shutdown were 10, 50, and 230 min. 
The mean lifetimes and absolute yields of delayed and photoneutron groups 
were found by analyzing each decay curve into a sum of decreasing exponentials 
and using the theory of pile kinetics. The yields were corrected to infinite irradia- 
tion time. These results are compared with data obtained by others. The compari- 
son shows that the energetic y-rays of fission products which produce photo- 
neutrons in the heavy water are absorbed to the extent of about 26% in the 
heavy water and 74% in the uranium rods. A table of recommended mean life- 
times and absolute yields of delayed and photoneutrons is given for use in the 
kinetics of heavy water piles such as ZEEP. Some quantitative data on the 
effects of the control plates and shutoff rods are included. 


INTRODUCTION 


The threshold for the photodisintegration of the deuteron occurs at 2.230 Mev. 
(2). During steady operation of a heavy water pile, photoneutrons are produced 
by the prompt y-rays of fission, by the y-rays from the radiative capture of 
neutrons in the various materials in the pile, and by the y-rays from the accumu- 
lated fission products if the energies of the y-rays are above this threshold. The 
last-mentioned source depends on the previous history of the pile. On the whole, 


these contributions are small in comparison with the neutrons from fission. 
However, after shutdown the neutron intensity decreases with time in a com- 
plicated manner owing to the delayed neutrons from fission and the photo- 
neutrons. These sources of neutrons are then dominant in the pile behavior. 
Information on the delayed and photoneutrons is therefore of practical interest. 
The delayed neutrons are evaporated by certain fission fragments after they 
have undergone §-decay. The intensities and periods have been thoroughly 
studied with the aid of cyclotrons and later of graphite piles (4, 14, 19, 22). We 
shall adopt the results given in a paper by Hughes, Dabbs, Cahn, and Hall (14). 
The photoneutrons produced by the energetic y-rays of fission products should 
possess the half-periods of these products, which may amount to minutes, hours, 
etc. In 1944 Elliott (5) showed that 75-day and 10-day old irradiated uranium 
produced photoneutrons in heavy water. The y-emitting fission product was 
identified as La'*° by chemistry and the half-period of the parent—300 hr. Ba!*”. 
The neutron yield was measured in these experiments with the integration 
apparatus described by Elliott, Hincks, and May (7). Later, Elliott (6) showed 
directly with a 8-ray spectrometer that La'*® emits about 0.05 quanta of energy 
1 Manuscript received October 5, 1953. ' 
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2.55 Mev. per disintegration. The cross section for the photodisintegration of 
deuterium at 2.55 Mev. was found to be (1.4 + 0.1) 10-*5 cm.? per deuterium 
atom. It is not practicable to isolate short-period fission products and make 
similar studies of photoneutron yields. 

A decay curve for the photoneutron intensity in the Argonne heavy water pile 
after shutdown was analyzed by May (17). The mean lifetimes reported were 
17.3, 78, and 238 min. Fenning (8) measured the photoneutron intensity in 
ZEEP with a boron trifluoride chamber from 30 min. to 24 hr. after shutdown of 
a 225-min. operation at 50 watts. The derived mean lifetimes were 30, 128, and 
275 min. Hincks (13) has examined these decay curves and found that the 
disagreement can be removed if certain justifiable changes (allowances for long 
period background and counting losses) are made. The mean lifetimes obtained 
by Fenning are then supported by the work at the Argonne Laboratory. 

The purpose of the experiments described in this paper was to make careful 
measurements of the intensities and mean lifetimes of some of the photoneutron 
groups, especially those that just escape being masked by the delayed neutrons. 
It was realized that the photoneutron groups might be numerous and that an 
analysis of a composite decay curve might be empirical only. The experiments 
were carried out on ZEEP in December 1946 and January 1947. Shortly after- 
wards, similar work was reported by Bernstein, Preston, Wolfe, and Slattery (3) 
and by Spatz, Hughes, and Cahn (24). The results of the three experiments will 
be discussed together in a later section. 

DETAILS OF ZEEP CONSTRUCTION 

The ZEEP (Zero Energy Experimental Pile) was designed and built at Chalk 
River in 1945 by a group* headed by L. Kowarski. Since only an incomplete 
account of its construction has been given (21), further details are given in the 
present paper. 

The aluminum tank that holds the heavy water is cylindrical in shape, has a 
“dished” bottom, and a flange at the top from which it is suspended. The 
thickness of the tank wall is } in. and of the bottom } in. The outside diameter is 
6 ft. 91 in., and the radius of curvature of the bottom is 13 ft. 6 in. The length of 
the cylinder to the joint with the dished bottom is 8 ft. 5} in., and the over-all 
length is 8 ft. 112 in. ; 

The uranium metal is in the form of slugs, 1.285 in. in diameter and 6 in. long. 
These are jacketed in aluminum ‘“‘stockings’’ of inside diameter 1.295 + 0.005 
in., wall thickness 0.040 + 0.003 in., and length 9 ft. 63 in. Each “‘stocking”’ 
contains nine slugs. A total of 148 rods, made up in this way, form the lattice, 
which is square with a spacing of 6 in. between rod centers. The rods hang 
vertically by a gimbal arrangement, supported by the steel beams that bridge 
the top of the tank. The lower ends of the rods lie in the horizontal plane con- 
taining the joint of the cylindrical tank with the dished bottom. 

The graphite reflector is 24 ft. thick under the tank and 3 ft. thick around it. 
The latter portion extends upwards to within about 1 ft. of the top. There are 


*L. Kowarski, C. W. Gilbert, C. N. Watson-Munro, G. J. Klein, E. P. Hincks, F. W. Fenning, 
A. H. Allan, G. Fergusson, and D. B. Nazzer. 
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two horizontal tangential channels through the graphite. The control plates are 
cadmium-plated steel sheets, suspended by wires in the annular gap, about 
13 in. wide, between the tank and the graphite. The plates are curved to conform 
to the curvature of the tank, and are spaced around the tank at approximately 
90°-intervals. Three are linked to move up and down together and to provide 
the coarse control, while the fourth can be moved independently and provides 
the fine control. The dimensions of the coarse control plates are 20 in. X 12 in., 
and those of the fine control plate 20 in. X 43 in. The two'sets of plates are 
moved by manually turning handwheels at the main control board, while selsyn 
dials indicate their positions. 

The shutoff rods are cadmium-plated steel tubes completely enclosed between 
aluminum tubes, and are suspended on wires inside and near the axis of the 
tank. The approximate length of the shutoff rods is 1 ft. These are arranged in 
two rows (A + B) of four each, spaced 8 in. apart within the row, while the 
two rows are also 8 in. apart. The rods can be raised or lowered manually by 
handwheels at the control board, and are normally in the ‘‘up”’ position when 
the pile is operating and in the ‘‘down’’ position when the pile is shut down. 
Intermediate positions are not used. Either set can be released to descend under 
gravity into the heavy water in about one second, either by push button or by a 
“trip” of the amplifier attached to an ionization chamber. 

In order to introduce neutron detectors, etc. into the tank for activation, a 
central aluminum “‘thimble’’, 24 in. O.D. and 2.17 in. I.D., and 21 aluminum 
“thimbles”, 1 in. O.D. and wall 0.022 in. thick, are provided. These are spaced 
at 6 in. intervals along two diameters at right angles, and pass through the mid- 


points of lattice cells. These tubes are usually filled with heavy water of inferior 
purity. The heavy water in the tank is 99.76 atom % D. 


EXPERIMENTAL DETAILS 


(a) Neutron Counting and Recording Equipment 

!t was desired to measure the neutron intensity in ZEEP during steady opera- 
tion at as high a power as permitted by the radiation hazard, and after shutdown 
(with all cadmium rods) over a duration of hours. The operating power level was 
chosen to be 50 watts, as indicated by one of the ion chambers. Since the neutron 
flux at 50 w. was of the order of 8 & 10° times the flux after a long shutdown, 
three neutron counters of very different sensitivities were arranged at suitably 
chosen points in the pile to follow the flux changes occurring in these experi- 
ments. To obtain the flux as a function of time over the periods of operation and 
shutdown, an intercalibration of the relative sensitivities of the three counters 
in their chosen positions was carried out and all counting rates (i.e. relative 
neutron fluxes) were expressed in terms of that of the most sensitive neutron 
counter. 

The counters used were designed and built by B. H. Flowers (9) and N. Veall 
for the measurement of slow neutron fluxes in a pile. The construction of a coun- 
ter is shown in Fig. 1. All the counters had the same physical dimensions, and 
were operated in the proportional region with a positive potential of above 1400 
volts on the central wire. Further details are given in Table I. 
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6 9 4 8 7 13: 12 0 
PROPORTIONAL NEUTRON COUNTER 


Fic. 1. Neutron counter. 

Over-all dimensions: length 7-3/8 in., diameter 1 in. Parts are of aluminum, with the 
exception of the end plugs and screws 2, which are of brass. 1—Threaded ring for applying 
pressure to rubber gaskets; 2—end plugs; 3—section of inner sleeve for clamping platinum 
strip supporting uranium or boron; 4—sections of inner sleeve; 5—outer tube; 6 and 7—end 
caps; 8—rubber gasket; 9—Kovar seals; 10—Amphenol plug; 11—brass contact plug; 12 
glass sleeve; 13—piano wire spring. The central wire is 8 mil tungsten. 


TABLE I 
NE UTRON COUNTERS 


Lab. No. | Neutron detector wh Gas filling | Operating voltage (v.) 





c 12 0. 14 mgm. u | 15 cm. -” 1320 
| | 5 cm. / | 
| | 

C8 1 mgm. U 15 cm. CH, 1560 
enriched in 5 cm. A 
| [236 | | 

c2 | B-coating 7.5 cm. CH, 1320 
2.5cm. A | 


The least sensitive counter (No. C 12), containing a foil of natural uranium 
and operating as a fission counter, was located in the central thimble of ZEEP 
somewhat below the position of maximum flux. It was used during the period of 
operation at 50w. and for about two minutes after shutdown, at which time the 
counting rate became too low to yield significant data. A similar counter (No. 
C 8), containing a foil of uranium enriched in U?**, was placed about two feet 
above the natural uranium counter. It was used chiefly during the first 15 min. 
after shutdown. The most sensitive counter (No. C 2), containing a boron coat- 
ing as neutron detector, was placed in one of the tangential channels in the 
graphite reflector about one foot from the wall of the tank. At this point the 
neutron flux was roughly one-fifth of that in the central thimble. This counter 

- provided a record of neutron flux from four minutes after shutdown to the end 
of each run. 

The electronic equipment was as shown in the block diagram, Fig. 2. This 
equipment was duplicated to allow any two of the three counters to be Se 
simultaneously. The voltage characteristics and bias curves of the counters were 
essentially the same as those obtained by Flowers (9). 
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Fic. 2. Block diagram of electronic equipment. 


(b) Method of Recording Data 

Two weeks of ZEEP time were allocated to these experiments, one in Decem- 
ber 1946 and the other in January 1947. A total of 12 runs were carried out. In 
addition, counting losses were investigated after each of the two series of runs. 

Before the experiments were begun all substances for the production of radio- 
nuclides for tracer work were removed from the pile. A shutdown period of at 
least eight hours preceded each series of runs in order to reduce the effect of 
previous operations. Activation periods of 10, 50, and 230 min. were chosen as 
suitable for bringing up the activities to various fractions of their saturation 
values. 

In most of the runs the counts in units of 128 were recorded on moving paper 
by two Esterline-Angus pen recorders. The speeds chosen were 5, 80, and 300 sec. 
per in. 

The power during the steady operation of ZEEP was monitored independently 
with a quartz-fiber electroscope* (Type TQQ) located some thirty feet from the 
side of the pile. Although this measured ionization by y-rays, it provided a more 
sensitive indication of flux changes than the power-indicating ion chamber in the 
graphite reflector. The rates of deflection, when corrected for changes in density 
of the air in the electroscope, were related to the counting rates in the uranium 
fission counter. 

The fine control plate of the pile was usually kept at a fixed and well-elevated 
position during operation and shutdown. The power was kept steady with the 
coarse control plates. A few seconds before shutdown, the coarse control plates 
were moved slightly to a predetermined position so that, after shutdown, the 
multiplication constant of the reactor would have the same value in all runs. At 
zero time for recording the decay of the neutron sources both sets (A + B) of 
the cadmium shutoff rods were dropped into the pile. The end of recording 
occurred from 3 to 43 hr. later. Each run provided the background for the next. 
(c) Counting Losses 

The observed counting rate m’ (in units of 128 per minute) of a counter and 
its associated equipment was related to the rate of drift P of a TQQ electroscope 
placed between the water shield and the graphite reflector of the pile. The power 
was varied from | to 50 watts. The rate of drift was assumed to be directly pro- 


*Designed by Dr. H. Carmichael, Chalk River Laboratory. 
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portional to the neutron flux and hence to the true counting rate of the neutron 
counter. This assumption appeared to be valid down to a power of about one 
watt, below which the effect of long-lived y-activity from previous runs became 


noticeable. 
The values of n’/P were plotted against P to give a curve that could be fitted 


by an expression of the form 
n'/P =a — bP + cP?, 


where a, 5, and ¢ are constants. Since aP represented the true counting rate n, 
the relation could be rewritten as 


a a 


or n’ =n — Bn? + Cn’, 


from which a graph of (x — n’)/n’ asa function of n’ was drawn and used to make 
the corrections for counting losses. Table II contains the values of the two 
constants B and C so determined. 
TABLE II 
CORRECTIONS FOR COUNTING LOSSES 


December 1946 January 1947 


Cx 107 


Counter B X 108 
Natural U 1.14 { Not repeated 
Enriched U | 0.54 . 

Boron 1.07 


The surprisingly large counting losses were probably caused by an unneces- 
sarily long dead time in the discriminators. For the counting data on which 
most of the decay analysis depends, the losses were 5% or less. 


(d) Treatment of Recorded Data 

The following procedure was adopted to extract the recorded data from the 
paper charts. The counting rate of the natural uranium counter, prior to shut- 
down, was established by tabulating the time on the chart corresponding to 
every hundredth count of 128. For a few seconds before and for each of the dozen 
or so counts after shutdown the time for each count (of 128) was determined. 
The instant of shutdown was thereby established to one second. In the case of 
the enriched uranium counter, the time corresponding to every tenth count 
(of 128) was recorded for the first minute following shutdown, the time for every 
fifth count for the next two minutes, and for every count for an hour thereafter. 
Beyond this time the counting rate was too low to be useful. In the case of the 
boron counter, the time corresponding to every hundredth count (of 128) was 
recorded for the first 90 min., every fiftieth count for the following day, and every 
tenth count thereafter. 

Wherever possible the counting rate was determined from the time for 100 
counts (12,800 particles), and the rate correlated with the mean time for the 
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counting interval. Where the counting rates were low or where the decay period 
was changing rapidly this procedure was not possible, and the rate was deter- 
mined from 50, 10, 5, or even 1 count with the concomitant loss in statistical 
accuracy. For the natural uranium counter prior to shutdown and for the boron 
counter throughout the run, from 50 to 100 counts were used for each deter- 
mination of counting rate, which should be accurate to about 1%. Less accuracy 
was obtained in determining the rates with the enriched uranium counter in the 
interval from two to six minutes subsequent to shutdown, when the effective 
half-period of decay was short and changing rapidly and the rates were too small 
to use time intervals short compared to the half-period. Under these adverse 
circumstances the rates were based on 5 counts (640 particles) and a decay curve 
was drawn. This curve was then corrected for errors introduced by measuring 
the rates over intervals comparable to the half-period. 

Counting loss corrections were finally applied to all the rates, and the corrected 
counting rates were plotted against the times. Fig. 3 shows part of a decay 
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Fic. 3. Observed counting rates as a function of time after a typical 10-min. run (No. 11 
Shutdown of the pile occurred at time zero. 


curve* for a typical run (No. 11). By intercalibration of the counters in the 
manner to be described in the following paragraph, it was then possible to 
express all counting rates in terms of the rates registered by the boron counter. 
Fig. 4, drawn from the data of runs No. 8, 9, and 12, shows in a convincing 
manner that the y-emitting fission products are built up to various degrees 
during pile operations of different duration. 

*The slight rise in the power just before shutdown was caused by the small shift in the position 


of the coarse control plates already referred to. This rise in power has been omitted in drawing 
Fig. 4 
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O 24 48 72 96 (20 144 68 2 216 240 264 288 32 3 
MINUTES AFTER SHUTDOWN 


Fic. 4. Counting rates (referred to the boron counter) as a function of time after steady 
operation of the pile for times of 10, 50, and 232 min. (runs No. 8, 9, and 12). 


(e) Intercalibration of Counters 

From the graph of counting rate against time for a run, such as Fig. 3, the 
ratios of the rates of the boron and enriched uranium counters and of the 
enriched and natural uranium counters were directly determined. In addition to 
the counter sensitivities, these ratios include the normal spatial variation of 
neutron flux and shadow effects of the control plates and shutoff rods. The ratio 
of the counting rates of the boron and enriched uranium counters remained 
constant to 1% while the boron counting rates fell from 600 to 60 counts (of 128) 
per minute. This result checks the counting loss corrections for the boron 
counter. No similar check was possible with the uranium counters. Their ratio, 
obtained from graphs such as Fig. 3 and from some direct comparisons at low 
power operation, is regarded as known to 3 or 4%. 

Table III contains the intercalibration ratios and other information for the 
various runs. The values of the counting rate m» of the boron counter at 50 w. 
were calculated from the rates of the natural uranium counter using the inter- 
calibration data. These values of m» are required in finding the absolute inten- 
sities of the delayed and photoneutron groups. We have no satisfactory explana- 
tion of the discrepancy in the values of m» for the December and January runs. 
The counters were removed between the two series of runs and then replaced but 
not exactly in the same positions. Moreover, the control plates were not at the 
same positions in the two series. Shadows cast either by the shutoff rods on the 
counters or by the control plates on the ion chamber may have been responsible 
for the discrepancy. The absolute yields of the delayed and photoneutron groups 
are therefore rather uncertain; the relative yields should be more accurate. 

(f) Measurement of the Change in k, at Shutdown 

The effective multiplication constant, k,, of a pile is unity under conditions of 

steady operation and is less than unity in all shutdown conditions. As will be 
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TABLE III 
INTERCALIBRATION OF COUNTERS 








_ Rum No.1 — 7 Runs No. 8 — 12 
Ratio B/enrichad U 
Ratio enriched U/natural U 
Rate* natural U at 50 w. 
Ratef TQQ at 50 w. 
Calculated rate* B at 50 w. (no) 
Ratio B/TQQ at 50 w. 
(¢ + 61)/no | 
*In units of 128 particles per minute. 
{In divisions per minute. 
The probable errors shown have been derived from the residuals. 


| 
December 1946 January 1947 
| 
| 





shown in the section on pile theory, the calculation of the fraction of the neutron 
population represented by delayed or photoneutrons requires a knowledge of 
6 = 1 — k, for the subcritical condition under which the data following shut- 
down were recorded. This value of 6 will be designated 6,. Two independent 
methods of measuring 6; were used. 

The first method requires a measurement of the ratio of the sudden initial 
drop in power accompanying the release of the shutoff rods after a long period of 
steady operation to the residual power immediately after shutdown. For this 
purpose, the power-indicating ion chamber was used to record automatically on 
an Esterline-Angus recorder. From equation [16] it will be seen that this ratio 
is 6,/c, where c is the fraction of neutrons of all types that are delayed. Using the 
value c = 0.0080 (to be justified later), the results of three trials are presented 
in Table IV. 

TABLE IV 
DETERMINATION OF 6; 
Method 1 


Trial No. Initial drop Residual 
in power _Power 
40.6 
1.47 
1.48 


Method 2 


| | 6,(° A ) 
Ah 3 2 ‘8; [en 
(cm.) >m.) | From 5o | From 8s 


2.62 : | >. 36 3. 2.54, | 2.34 
2.04 3.6: 3.37 4.4 2.62 2.46 


iuecaiee: 2. 49 


In run No. 4 the control plates were inadverte ntly “moved after shutdown. For these shutdown 
conditions 6 was found to be 3.60%. 
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The second method of obtaining 6, involves a comparison of the neutron 
counting rate in the normal shutdown condition with those in other subcritical 
conditions described by values 62 and 63; (smaller than 6:). Subsequently, 5: and 
6; were obtained from the increases in the depth of the heavy water necessary 
to bring the reactor to the critical condition (6 = 0). 

In the subcritical pile the neutron population is directly proportional to the 
strength of sources (photoneutrons and spontaneous fission), which are nearly 
constant several hours after shutdown. The neutron counting rate at a point in 
the pile is inversely proportional to 6 for constant sources. Hence, if the counting 
rates in the two subcritical conditions described by 6; and 62 are m and nz 
respectively, the experimental ratio 11/mz is 62/61. 

The subcritical condition described by 6; is the normal one with shutoff rods 
“down” and the control plates in a preselected position very near that for 
normal operation, while 62 and 6; refer to slightly subcritical conditions where the 
shutoff rods are ‘“‘up’’ and the control plates are moved to more effective posi- 
tions. With the pile in the conditions described by 62 and 63, the increments 
Ah» and Ah; in the depth of heavy water necessary to restore the critical condi- 
tion were measured. Fig. 5 is included to illustrate the relation between the 
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SELSYN READING FOR COARSE CONTROL PLATES 


Fic. 5. Manometer height in centimeters of heavy water as a function of the height of the 
coarse control plates for a self-sustaining chain reaction. Curve C: fine control plates set at 
190. Curve D: fine control plate set at 114. The lower ends of the uranium rods are at approxi- 
mately 17.2 cm. on the manometer scale. 


positions of the control plates (as read from the selsyn dials) and the depth of 
heavy water for a self-sustaining chain reaction in ZEEP. The values of 6, and 6; 
_were calculated, using equation [17], from Ah2 and Ah; and the known pile 
constants. Hence, two values of 6; were derived. From the results in Table IV, 
the mean value of 6, for the shutoff rods (A + B) is 2.53%. 
PILE THEORY 

The theory of the kinetics or time behavior of the chain reaction was developed 

and put on record by Nordheim (18) and others (10, 23). 
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For each neutron absorbed in the pile, &, neutrons are made available to carry 
on the chain reaction. Of these, k,(1 — c) are immediately available and k,cy; 
are stored in nuclei of type 7. The photoneutrons may be included with the 
delayed neutrons in spite of their different origin. The total number of free 
neutrons, 7;, in the pile and the stored neutrons, r;, associated with various 
delays change with time ¢ according to the following differential equations: 


dr; ke(1 — ¢) 


oa eee Be\t — C) F 
[1] a - te ts - +h Ni lis 


[2] os os Reals wh Be 
where 1 = 2, 3, 4,..., 7. 
The following symbols are being used: 
k, = effective multiplication constant of the pile, 
c = total fraction of neutrons delayed, 
u; = fraction of delayed neutrons of type 7, 
\; = transformation constant associated with delayed neutrons of type 2, 
and + = mean lifetime of a neutron in the pile. 
The following changes in notation lead to shorter differential equations: 


[3] 5=1-k,, 
[4] Ar = (6+ — 68)/7, 
[5] and A; = (1 — 8)cp,/r. 
Equations [1] and [2] become 

dri 


[6] dl ~My ti-+ Lire 


z dr; 
[7] a a tet Bee 

During the period of activation of duration 7 the value of r; is kept constant. 
It may be set equal to m») as measured with the neutron counter. The delayed and 
photoneutrons, r;, are built up from negligible values to fractions of their 
saturation values, F;, given by 


[8] F,:= 1 —e™*". 
The solutions of equations [6] and [7] under these conditions are 
[9] ri = No, fe (A; F;/d;)no, 


which are the initial conditions for the solutions of equations [6] and [7] following 
shutdown with the cadmium rods. The value of 6 is changed to 6; = 2.58% 
within one or two seconds of steady operation. During this interval 6 and ), 
change rapidly and come to their steady values defined above. An exact solution 
is not possible if the production of further stored neutrons is included. However, 
Xi is much larger than any of the ,’s, and for this reason a good approximation 
to the solution can be found for this perturbation problem. 
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The solution for 7; may be written 


[10] n= > Ae 
i=1 


The values of A; and a; may be expressed in terms of the y;’s and \,’s through 
the following equations: 

(11] a; = A,(1 + €), 

where e; < 1, 


i=2 Ai — ri ’ 
Ai Ai 


[13] ae Fing = At = ri f 


x a 
[14] ‘ee ( + Fino ’ 


and 


[12] 


[15] cu; = (c + 6) Fe. 

The experimental decay curve may be analyzed into a sum of exponentials 
of the form of equation [10]; and the values of A; and a; so obtained may be 
related to the desired neutron yields cu; and transformation constants \; through 
the use of equations [11] to [15]. 

Equation [15] may be directly used to determine the sudden relative drop in 
power that occurs in the first few seconds of shutdown. This ratio is evidently 


[16] |. - ¥ al / sent 


Rumsey (20) has considered the modifications to pile theory that are necessary 
for the interpretation of experiments involving the kinetics of piles with reflec- 
tors. The following formula pertaining to ZEEP in a subcritical condition was 
derived: 


i a (2h + Ah)Ah rckoceers 

[17] 6 = Hh > Ah)? (M* + 2K°L'L,’). 

For ZEEP the square of the slowing-down length, L,?, and the square of the 
diffusion length, L*, for neutrons in the lattice are both approximately 118.5 cm.? 
The migration area M? = L,? + L? = 237 cm.? From the spatial distribution 
of slow neutrons in ZEEP, we have found that the Laplacian K? is 794 & 10-® 
cm.~? and the vertical distance h between the two horizontal planes where the 
neutron density vanishes is 148.5 cm. (15). The value of M? + 2K°L°L,? is 
therefore approximately 260 cm.? Equation [17] determines the values of 6 for 
the various subcritical conditions discussed in the previous section. 


ANALYSIS OF DECAY CURVES 


Each neutron decay curve was analyzed as the sum of a number of decreasing 
exponentials (equation 10) by the well-known “‘peeling’’ process. Figs. 6-8 show 
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Analysis of Run 4 showing 

(a) Determination of Background 
(b) 6.0 hour period 

(c) 2.40 hour period 


T=2.40 hr. 


15 20 25 
Time in hours from shutdown 


Fic. 6. Part of the analysis of a 50-min. run. 


Experimental 
Curve 
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3. 7. Part of the analysis of a 50-min. run. 
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Fic. 8. Part of the analysis of agp0-min. run. 


the application of the method. When, as in these cases, there are many com- 
ponents, the method is subject to serious cumulative errors. Attempts to mini- 
mize such errors have been made by varying the activation time and by having 
the curves analyzed independently by two persons. 

The limited time allotted to these experiments made it impossible to follow 
each run until the counting rate reached the apparently steady value determined 
by long-period sources such as 12.8-day Ba—La'*® and by spontaneous fission. 
The backgrounds used were based on the values obtained in the two longest 
runs, Nos. 4 and 12, in which the counting rate had become apparently constant. 
These values of 0.37 and 0.69 units of 128 per minute (measured at different 
5's) are approximately three times those estimated for spontaneous fission alone. 

THE TOTAL FRACTION OF NEUTRONS DELAYED 

Since the analyses could not be extended to give directly the intensities cu; of 
the very short lived groups of neutrons, it was necessary to choose ¢c (= )ocyu;) 
in order to obtain the absolute yields of the various photoneutron groups ob- 
served. The value c = 0.0080 was arrived at in the following manner. The total 
fraction of fission neutrons that are delayed was found to be 0.00755 by Hughes 
et al. (14). Bernstein et a/. (3) found that the total yield of photoneutrons was 
16% of the saturated delayed neutron activity when their measurements were 
‘corrected to an infinite volume of heavy water around an irradiated sample of 
uranium having negligible y-ray absorption. Evidence to be presented later 
shows that about 26% of the energetic y-rays escape from the ZEEP rods. The 
photoneutron contribution to c in ZEEP is therefore 26% XK 16% X 0.00755 
= 0.00031. Hence c = 0.00786. If we use a photoneutron yield estimated at 


28% of the delayed neutrons by Spatz ef a/. (24) in an experiment similar to that 





150 CANADIAN JOURNAL OF PHYSICS. VOL. 32 


of Bernstein, the photoneutron contribution to c is 26% X 28% X 0.00755 
= 0.00055, and hence c = 0.00810. Gilbert and Fergusson (11) found c = 0.0083 
in an experiment in which the ZEEP power was modulated with a rotating cad- 
mium shutter. We have adopted c = 0.0080. 


EXPERIMENTAL RESULTS 
The final values of the mean lifetimes and yields cu; of the various groups of 
delayed and photoneutrons are shown in Table V. The probable errors shown 
have been calculated from the residuals. 


TABLE V 
SUMMARY OF MEAN LIFETIMES AND YIELDS 
OF DELAYED AND PHOTONEUTRONS 


Mean lifetime Yield X 105 


5.8 +0.3sec. | 
30.0 + 0.5 sec. 
81.9 + 0.5 sec. 

5.45 + 0.08 mina. | 
22.8 +0.3 min. 
72.4 +1.5 min. 

2.48 + 0,02 hr. 

6.15 + 0°06 hr. 


DISCUSSION 

The first three groups of neutrons listed in Table V are mostly delayed neu- 
trons from fission. The mean lifetimes 5.8, 30.0, and 81.9 sec. agree with 6.5, 
31.8, and 80.2 sec. found by Hughes et a/. (14) in their study of the delayed 
neutrons. The yields 257, 172, and 33 (relative to 10° fission neutrons) are also 
in approximate agreement with 213, 166, and 25 measured by Hughes. The 
photoneutrons produced by the y-rays of short-period fission products in ZEEP 
are therefore strongly masked by the delayed neutron groups. 

For the purpose of comparison, Table VI presents a summary of the most 
detailed investigations of delayed and photoneutrons. The yields of the photo- 
neutron groups have been normalized to 166 X 10-5 for the 31.8-sec. delayed 
neutron group (14), which is sufficiently prominent for the purpose. The mean 
lifetimes 6.15 hr. and 2.48 hr. found in the present work agree with 6.3 hr. and 
2.4 hr. found by Bernstein ef a/. (3). Marked disagreement in the mean life- 
times occurs at the weak group of photoneutrons with a mean lifetime of 72.4 
min. observed in the present work. Since errors in analysis are cumulative, the 
disagreement in lifetimes persists into the stronger 22.8 min. and 5.45 min. 
groups. The delayed neutron groups at 80.2, 31.8, and 6.5 sec. are sufficiently 
strong that even small admixtures of photoneutrons and cumulative errors do 
not prevent our obtaining agreement with Hughes ef al. (14). 

The photoneutron yields in ZEEP are very much less than those found by 
either Bernstein or Spatz in their investigations using thin y-ray sources of 
irradiated uranium. Since the short-lived groups are strongly masked by the 
delayed neutrons in our work, an over-all comparison is less likely to be signi- 
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TABLE VI 
A COMPARISON OF VALUES OF MEAN LIFETIMES AND YIELDS OF DELAYED AND PHOTONEUTRONS 


Mean lifetimes Yields relative to 10° fission neutrons 





Delayed 


Spatz | Bernstein 
(24) (3) 


Present 
| work 
| Delayed 


Hughes 
(14) 


Spatz 
(24) 


Bernstein 
(3) 


Present 
work 


Photoneutrons only | and photo- 


only | neutrons 


0.07 sec. 
0.62 sec. | 
2.19 sec. 
112 
| 5.8 sec. | 
130.0 sec. | 5 


79 


6.50 sec. 
31.8 sec. | 


80.2 sec. \81.9 sec. 
3.4 min. | 5.45 min. 
min. |22.8 min. 
min. |72.4 min. 


2.4hr. | 2.48 hr. 





12.4 
3.6 
2.0 
1.6 





0.5 


>.3 hr. 6.15 hr. 


0.123 


120.5 
16% 


Total yields 
Total yields 





47 
6% 


211 
28% 


755 
100% | 











* Normalized to 166 to agree with Hughes (14). 
{This sum represents only photoneutron contributions; the 5.8 sec. and 81.9 sec. yields are 
assumed to include photoneutron contributions of 248 — 213 = 35 and 32— 25 = 7 respectively. 


ficant than one involving only the longer lifetimes. The total yield for the five 
longest lifetimes in ZEEP is 4.67 X 10-5 of fission neutrons, which is 26% of 
Bernstein’s value for comparable periods. This reduction of 74% is in agreement 
with the following calculated absorption of energetic y-rays in the uranium rods 
of ZEEP. Since the absorption coefficient for y-rays in the energy range 2.23 to 
3 Mev. changes very slowly with energy, it will be sufficient to select the value 
4.45 X 10-? cm.?/gm. obtained by Hatton (12) of this laboratory for the 
absorption of the 2.62 Mev. thorium C” y-rays in uranium. Using this coefficient 
and Kushneriuk’s graph (16) of the absorption of y-rays from uniformly dis- 
tributed sources in a rod, the absorption is found to be 67%. The neighboring 
rods in the lattice absorb about 12% of the y-rays that escape from the source 
rod. The total absorption of y-rays in the uranium rods of ZEEP is therefore 
71%, which is in agreement with the photoneutron reduction of 74%. 

For calculations in pile kinetics, we have adopted the mean lifetimes and 
yields of delayed and photoneutrons measured by Hughes (14) and Bernstein 
(3). The photoneutron yields were reduced to 26% of the measured values to 
take account of y-ray absorption in the ZEEP rods. Table VII contains the rec- 
ommended information, which has been used at Chalk River. Some confirma- 
tion of the essential correctness of these data has been found by Bayly (1), who 
obtained good agreement between measured and calculated periods of ZEEP 


when divergent. 
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TABLE VII 


DELAYED AND PHOTONEUTRONS IN ZEEP. RECOMMENDED DATA 


Mean lifetime 
Ti (sec.) 


Delayed 0.07 
neutrons .62 
.19 
».50 

8 


9 


- 


Photo- 3.6 
neutrons 59.0 
204 
660 
2340 
8640 


ai X 108 


31.7 
107. 
306. 


5.50 


23 
07 


659 
743 
. 102 


22680 | 
273600 | “0327 


c = Seu; = 0.8044%, 


It has not been possible to identify unambiguously the hard y-ray-emitting 
fission products responsible for the production of photoneutrons. 
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DISINTEGRATION OF Au'®’! 


By L. G. Etitrotr, M. A. PREsSTON,? AND J. L. WOLFSON 


ABSTRACT 


In addition to the 411.77 kev. y-ray transition, two y-ray transitions of energies 
676.5 + 0.8 kev. and 1088.9 + 0.9 kev. have been observed in Hg'®* following 
the disintegration of Au'®*. The intensities of these transitions were found to be 
(8.20 + 0.56) X 107% and (1.63 + 0.12) X 10-* per disintegration respectively. 
The K conversion coefficient of the 676.5 kev. transition was measured as 0.0224 
+ 0.0019, identifying the radiation as a mixture of magnetic dipole and electric 
quadrupole. The K conversion coefficient of the 1088.9 kev. transition was 
measured as 0.00450 + 0.00034, identifying this radiation as electric quadrupole. 

A 8-transition to the ground state of Hg'®’, of energy 1371 kev. and intensity 
(2.5 + 0.5) X 10-4 per disintegration, has been observed and identified as a 
third order forbidden transition. The shape of this spectrum has been studied 
above the end point of the intense 959 kev. 8-transition. Of the four combinations 
considered for the 6-decay interaction form, viz. SA, VA, ST, VT, it was found 
that the experimental shape can be fitted by any one of VA, ST, or VT, but not by 
SA. The log fot value for the transition is 11.7. 

A disintegration scheme is proposed in which the excited states of Hg'®* at 412 
and 1089 kev. are each assigned spin 2 and even parity and the ground state of 
Au!’ spin 3 and odd parity. 


I. INTRODUCTION 


Au'®* (2.69 day) decays by 8- emission, the main branch proceeding to the 
first excited state in Hg'®® at 412 kev., followed by the emission of a y-ray to 
the ground state which is partially internally converted. The end point of the 
main 6-spectrum has been measured by many observers (4, 15, 16, 35, 39), the 
most recent determination being that of Fan (6) who reported a value of 
961 kev. Langer and Price (15) have shown the spectrum to be of the ‘‘allowed”’ 
shape, consistent with the results of Ridgway (27) and Stump and Frankel 
(40) who find no @-y angular correlation. Renard (26) has shown that K 
capture if present occurs in <0.4% of the disintegrations, while Mims and 
Halban (20) have found that 8+ emission occurs in <3.3 X 10-°% of the 
disintegrations. 

The most precise determination of the energy of the 412 kev. y-ray in 
Hg'** has been made by Muller et a/. (22) who reported a value of 411.770 
+0.036 kev. In addition to this y-ray, two others of energies about 680 kev. 
and 1090 kev. have been reported by Cavanagh et al. (3), Pringle and Standil 
(24), Elliott and Wolfson (5), and Hubert (13). Cavanagh (2) has shown that 
the 412 kev. y-ray is in cascade with the 680 kev. transition. All workers agree 
that the 1090 kev. y-ray is a crossover transition. Elliott and Wolfson (5, 42) 
have observed a very weak 8~ transition to the ground state of Hg'®*, and 
Brosi et a/. (1) have observed a third weak 8~ transition of energy 290 kev. by 
counting 8-rays in coincidence with the 680 kev. y-ray. They conclude that 
this low energy transition proceeds to the 1090 kev. state in Hg!®*, 

The A, Z, and M conversion coefficients of the 412 kev. transition have been 
measured by many observers (6, 12, 13, 32, 35, 36, 39), the most recent mea- 
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surement being that of Simons (37), who has found the value for the K conver- 
sion coefficient to be in excellent agreement with the calculated value of Rose 
(28) for electric quadrupole radiation. Measurements of the half life of this 
transition have also been found to be consistent with assignment of electric 
quadrupole character. Values of <3 X 107" sec., (2.1 + 0.5) X 10-™ sec., 
and 5.5 X 10~" sec. have been obtained by Graham and Bell (9), Moon and 
Davey (21), and Malmfors* (19) respectively. 

In the present paper the results of a detailed study of the disintegration of 
Au!®’ are presented. The y-rays of energies 676.5 and 1088.9 kev. have been 
observed both in internal and external conversion and their intensities and 
conversion coefficients measured. The energy and intensity of the weak 8- 
transition to the ground state of Hg'®* have been measured, and the shape of 
the spectrum above 975 kev. has been studied. From the data obtained the 
spins and parities of the states in Hg!** and the ground state in Au'®* have been 
deduced. A disintegration scheme incorporating these results is proposed. 


II. APPARATUS 
Two magnetic lens spectrometers have been used in the studies. For the 
examination of 8-ray spectra and internal conversion lines, a spectrometer 
using two coaxial magnetic lenses was employed. The instrument is shown 
schematically in Fig. 1. The source is mounted on a thin backing and placed 
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Fic. 1. Double lens 6-ray spectrometer used in the study of 8-ray spectra and internal 
conversion lines. In these experiments only the outer three windings on each magnet were used. 








TO PUMP 


at the center of one of the coils. The detector, an anthracene disk 3 in. in di- 
ameter and about | mm. thick, is situated at the center of the other magnet 


*We have multiplied Malmfors’ value (1.6 X 10~"' sec. for the mean life) by a factor of 5 to take 
into account the statistical weight of the 412 kev. state. 
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coil. A lucite light pipe 30 in. long conducts the light flashes from the anthra- 
cene to a magnetically shielded 5819 photomultiplier. Ring focusing is em- 
ployed, the focus occurring at the annulus in the baffle nearest the detector. 
The distance from source to detector is 50 cm. A set of spiral baffles made of 
aluminum coated copper is installed, so that for a given direction of magnetic 
field the instrument selects electrons of one sign only. However, in the present 
application these baffles serve the purpose of reducing the transmission of 
electrons scattered by the annular baffles and walls of the vacuum chamber. 
Additional quantities of lead, not shown in the figure, are disposed in the 
chamber so as to reduce the background due to primary and secondary quanta. 
As adjusted in these studies, the instrument has a measured transmission of 
about 0.5% and the line width at half maximum intensity is 0.5% in momen- 
tum. The momentum measurements are made in terms of the voltage drop 
across a single manganin resistor, placed in series with the magnet current, and 
kept at a constant temperature in an oil bath. The magnet current is stabilized 
to better than one part in 5000. The A interval conversion line of the 411.77 
+0.036 kev. (22) transition was used to calibrate the spectrometer. 

Because of the low intensity of the high energy 8- transition it was felt 
necessary to investigate carefully the background due to scattered electrons 
above the end point of a simple 8-spectrum, since an appreciable contribution 
of scattered electrons from the intense 959 kev. spectrum could easily mask 
the true shape of the higher energy spectrum. The investigation was carried 
out by examining the region near and beyond the end point of the P® spectrum 
which is known to be simple. It was found that the counting rate in excess of 
counter background extended for three to four line widths beyond the end 
point of the spectrum, after allowing for the finite resolution of the spec- 
trometer, and was just detectable above the counter background. The correction 
applied to the measured high energy Au'®* spectrum on account of this effect 
amounted to less than 1% of the counting rate for the experimental point at 
975 kev. For the remainder of the data, the corrections required were com- 
pletely negligible. 

The background due to y-rays was investigated using a 35 mc. source of 
Au'®’ covered with sufficient beryllium absorber to remove all primary §-rays 
from the beam. The background counting rate which was measured at the 
same momentum values as used for the high energy Au'®* 8-spectrum was 6 
counts per minute, of which 2 counts per minute was the counter background. 
There was no variation of background over the momentum range of interest. 
With the 7-mc. sources used in the experiments the y-ray background was 
therefore 0.8 counts per minute. 

Another subsidiary experiment was performed to determine the amount of 
magnetic shielding required by the photomultiplier. The front face of the 
anthracene crystal was covered with thick aluminum and a Co® source was 
placed outside the vacuum chamber. The counting rate was then measured 
over the same range of magnetic fields as used in the study of the high energy 
8-spectrum. Sufficient shielding was then added until the counting rate became 
independent of the magnetic field. 
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For y-ray studies a single lens spectrometer was employed. The line width of 
this instrument is 1.5% in momentum at half maximum intensity, and the 
transmission 0.16%. The detector consisted of an anthracene disk 1 mm. thick 
and 4 mm. in diameter. A lucite light pipe conducted the light to a magnetically 
shielded 5819 photomultiplier outside the intense magnetic field. The studies 
were concerned exclusively with measurements of y-ray intensities by com- 
parison of the intensities of the corresponding K shell photoelectron lines 
externally converted in an 8 mgm./cm.? thorium radiator, 4 mm. in diameter. 
The end plate of the spectrometer was provided with an aluminum rod which 
projected about 13 in. into the vacuum chamber to the end of which the 
thorium radiator was affixed. A hole 4 mm. in diameter was drilled into the 
rod to within 1 mm. of the end, and sources enclosed in aluminum containers 
were inserted into this hole from outside the vacuum chamber so as to occupy a 
standard position relative to the thorium radiator. 

III. THE B-SPECTRUM 

The gold sources used were all 1 mm. in diameter mounted on 2.7 mgm. /cm.? 
thick aluminum foil. The aluminum foil itself was mounted on a dural ring 
with an inner hole 1? in. in diameter. After centering, the source was placed in 
the spectrometer as shown in Fig. 1. Sources of Au'®* of thickness 0.2, 1.0, 2.5, 
and 5.0 mgm./cm.? were used in different trials. 

Sources of thickness 0.2 mgm./cm.? were used in the study of the 959 kev. 
spectrum and in observing the internal conversion lines. For the detailed 
examination of the high energy spectrum above 975 kev., and for measuring 
the K/L ratio of the 1089 kev. y-transition, sources of thickness 2.5 mgm. /cm.? 
were used. An experiment was performed to demonstrate that for these pur- 
poses the 2.5 mgm./cm.? sources could be considered “‘thin’’. Counting rates 
were measured at 17 energy values from 440 kev. to 1350 kev. using sources of 
thickness 1, 2.5, and 5 mgm./cm.? of known relative activities. From these 
data it was found that above 850 kev., for a given source strength, the variation 
of counting rate with source thickness was negligible over the range of thick- 
nesses used. 

The spectrum above 2560 gauss-cm. taken with a source of thickness 0.2 
mgm./cm.” is shown in Fig. 2. Evidence was presented in Section II to show 
that the continuum above the end point of the 959 kev. B-spectrum was not 
due to electrons from the lower energy component scattered into the detector. 
In addition, however, the electrons were demonstrated to be of high energy 
by a simple absorption experiment in which the anthracene detector was 
covered with 140 mgm./cm.* of aluminum foil, and the counting rates were 
observed at spectrometer settings of 600 kev., 1005 kev., and 1025 kev. The 
second point lies at the peak of the K internal conversion line of the 1089 kev. 
transition, so that the electrons counted are known to be of high energy. The 
point at 1025 kev. lies between the K and ZL internal conversion lines of the 
1089 kev. transition, and at this energy only those electrons belonging to the 
high energy continuum are counted. At the 600 kev. setting it was found that 
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_ Fic. 2. Plot of the 6-spectrum of Au'®*® above 2560 gauss-cm. The Z internal conversion 
line of the 676.5 kev. transition, not shown, is about one sixth of the K line. 


the absorber caused the counting rate to be reduced by a factor of 5.3 + 0.5, 
while at the 1005 kev. and 1025 kev. settings the reduction factors were 2.0 
+ 0.1 and 2.0 + 0.2 respectively. Finally, using a more intense source, the 
high energy continuum was examined over a period of eight days, and it was 
found to decay at the same rate as the K internal conversion line of the 412 kev. 
y-ray transition. From the results of these experiments it is concluded that the 
high energy continuum constitutes a genuine B-ray spectrum resulting from 
the disintegration of Au'®*. 

The main component of the 8-ray distribution, at energies above 412 kev., 
was found to have the ‘‘allowed’’ shape, in agreement with the results of 
Langer and Price (15). A Fermi plot of this component yielded an extrapolated 
end point of 959 + 2.5 kev. An allowed Fermi plot of the high energy com- 
ponent yielded an extrapolated end point of 1370 + 4 kev., which agrees with 
the sum of the energies of the 411.77 + 0.036 kev. y-transition and the 959 
+ 2.5 kev. main 6-ray transition. This 8-ray transition is accordingly believed 
to proceed to the ground state of Hg'®*. A discussion of the shape of this spec- 
trum is given in Section VII below. An estimate of the intensity of the 1371 kev. 
component relative to the 959 kev. component was made by constructing the 
‘two spectra from their Fermi plots extrapolated to zero energy, and measuring 
the ratio of their areas. The intensity of the high energy component was found 
in this way to be (2.5 + 0.5) X 10-4 of the 959 kev. component. The corre- 
sponding log fof value is 11.7. 
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IV. INTERNAL CONVERSION SPECTRUM 


Internal conversion lines corresponding to three y-ray transitions have been 
observed. In addition to the K, Z, and M lines of the well-known 411.77 kev. 
transition, K and L lines of a 676.5 + 0.8 kev. transition and K, L, and M@ 
lines of a 1088.9 + 0.9 kev. transition have been observed. These are shown 
in detail in Fig. 3. The intensities of the K lines of each of the transitions have 
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Fic. 3. Internal conversion lines of the 676.5 and 1088.9 kev. transitions in Hg'®’. The 
vertical bars represent standard deviations. The intensities of the lines are proportional to the 
line areas divided by their corresponding momenta. 


been followed in the spectrometer over a period of three half-lives. The inten- 
sities of these lines, relative to the 412 kev. K internal conversion line, were 
found to be constant within the statistical errors, indicating that the transi- 
tions were associated with Au'*8 decay. The K—Z differences of both sets of 
lines were also measured. The values of 67.9 + 0.6 kev. and 68.0 + 0.5 kev., 
found for the 676 and 1089 kev. transitions respectively, constitute evidence 
that conversion takes place in the Hg atom and, in particular, that the tran- 
sitions do not follow K-capture. 

The ratios of conversion coefficients ag/az were found for each transition 
by comparing the intensities of the K lines with the intensities of the corre- 
sponding ZL lines. 

The values obtained were: 


y-Ray energy aK /az 


676.5 kev. 5.7+0.5 
1088.9 kev. 6.340.5 
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For the measurements on the lines of the 677 kev. transition a source 
0.2 mgm./cm.? thick was used, while for measurements on the lines of the 
1089 kev. transition, a source 2.5 mgm./cm.* thick was used in order to obtain 
better statistics in the measurement of the LZ line. The results quoted are the 
means of two separate determinations. 

The intensities of the K conversion lines of the 677 and 1089 kev. transitions 
relative to the intensity of the K conversion line of the 412 kev. transition were 
also measured, using a 0.2 mgm./cm.? thick source for all measurements. The 
values found were (5.91 + 0.24) X 10-* and (2.41 + 0.03) X 10-‘ for the 
676.5 kev. and 1088.9 kev. transitions respectively. These quantities, which 
are the means of two determinations made with different sources, are used in 
the determination of internal conversion coefficients in Section VI. 

It is to be noted that the difference of the energies of the 1088.9 kev. and 
676.5 kev. transitions, i.e. 412.4 + 1.2 kev., is in close agreement with the 
energy of the 411.77 kev. transition. It is accordingly believed that the latter 
transition is a crossover transition from a state in Hg'®* at 1088.9 kev. to the 
ground state in Hg!®, 


V. y-RAY INTENSITIES 


The intensities of the 676.5 kev. and 1089 kev. y-rays were measured by 
observing the K photoelectron lines from a thorium radiator, 8 mgm./cm.’? 
thick, using the single lens spectrometer described in Section II. The source 
was made from three disks of pure gold, each 10 mgm./cm.* thick, 3 mm. in 
diameter. The gold foils were separated from each other by two disks of 


“‘superpure’’ aluminum of the same diameter and 0.009 in. thick, and the 
assembly mounted in a small cup at the end of a ‘‘superpure’’ aluminum rod, 
and covered with an aluminum disk 0.009 in. thick. The source and mounting 
were then irradiated in the Chalk River reactor to an activity of about 350 mc. 
A similar rod of ‘‘superpure’’ aluminum was irradiated along with the source, 
and examined afterwards to verify that activities were not formed in the rod 
material which would interfere with measurements made on the Au!* source. 
A disk of gold 0.2 mgm./cm.? thick and 4 mm. in diameter was also irradiated 
along with the intense source and an absolute determination of its activity 
made by coincidence counting. The activity of the intense source*was then 
found by comparison with this source of known activity by y-ray counting 
with a scintillation counter. 

In order to determine the y-ray intensities accurately from the measured 
intensities of the K photoelectron lines it was necessary to calibrate the spec- 
trometer arrangement for intensities because of the variation of the photo- 
electron absorption cross section with energy. This was done at y-ray energies 
of 412, 602, 794, 1172, 1332, and 1360 kev. by using sources of Au!®’, Cs!4, and 
Co® of known activity. The Au'® source was the intense one described above. 
A Co® source of the same geometrical size was calibrated against a 
Co® standard. The y-ray flux of energy 1360 kev. from a Cs'* source, also 
of the same geometrical size, was found by comparing the intensity of its K 
photoelectron line with the intensity of the K line due to the 1332 kev. y-ray 
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from the Co® source. The relative photoelectric absorption cross section for 
the two quanta was taken to be 1.06, and this value was later checked against 
the complete curve of effective relative photoelectric absorption cross section 
vs. quantum energy. The intensities of the 602 and 794 kev. y-rays from Cs", 
relative to the 1360 kev. y-ray, were taken to be 33.4 + 1.7 and 32.4 + 1.7 
respectively. The latter value was obtained by direct comparison of the 794 
and 1360 kev. quantum intensities using a NalI(T1) scintillation spectrometer 
calibrated for intensity. The former value was then deduced assuming the 
order of emission of the 602 and 794 kev. y-rays proposed by Waggoner e¢ al. 
(41). These values* are, however, not sensitive to the details of the disintegra- 
tion scheme of Cs", 

The K photoelectron lines of the 676.5 kev. and 1088.9 kev. y-rays from 
Au!*’ are shown in Fig. 4. From the intensities of these lines, the y-ray inten- 
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Fic. 4. K photoelectron lines of the 676.5 and 1088.9 kev. y-ray transitions. The ordinates 
are the observed counting rates. The vertical bars represent standard deviations; where no 
vertical ‘bars are shown, the standard deviations are less than the radii of the circles. The 
intensities of the photoelectron lines are proportional to the line areas divided by their corre- 
sponding momenta. The dashed curves give the shapes of the background under the lines as 
found by examining the spectrum with the thorium radiator removed. 


sities relative to the 412 kev. y-ray were found to be (8.42 + 0.56) X 107% 
and (1.70 + 0.12) X 10-* for the 676.5 kev. and 1088.9 kev. y-rays respectively. 
When allowance is made for the internal conversion of the 412 kev. transition 
(a = 0.045), the values found for the intensities of the transitions per dis- 
integration, J, are: 

*In the earlier work (5, 42) the values used were 24.6 and 23.6 respectively. These were deduced 
from the ratio given by Waggoner et al. (41) for the intensity of the 1360 kev. y-ray to the B-ray 
intensity, using a B branching ratio of 30/70. Accordingly the quantum intensities of the 676.5 and 


1089 kev. transitions in Hg'*® were underestimated by about 25% and 10% respeciirely, with 
corresponding overestimates of the conversion coefficients. 
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Energy of transition I 
676.5 kev. (8.20 + 0.56) X 10°% 
1088.9 kev. (1.63 + 0.12) X 10% 


In deriving these values allowance is also made for internal conversion of 
the transitions as found in Sections IV and VI. The results quoted are the 
weighted means of three separate determinations using different Au'®* sources 
for each determination. 

VI. CONVERSION COEFFICIENTS 

The A conversion coefficients, ax, have been determined for the 676.5 kev. 
and 1088.9 kev. transitions, utilizing the data of Sections IV and V along with 
the theoretical value (28) of ag, 0.0317, for the 412 kev. transition which is 
taken to be £2 in character. There is no doubt as to this assignment for the 
412 kev. transition since all experimenters have found this to be the case on 
the basis of the measured value of ax (6, 12, 13, 32, 35, 36, 39), the ag/az ratio 
(6, 11, 12, 32, 35, 36, 37, 39), and the lifetime (9, 19, 21) of the transition. 
The conversion coefficients were found to be: 


Energy of transition aK 


676.5 kev. 0.0224 + 0.0019 
1088.9 kev. 0.00450 + 0.00034 


The calculations of Rose (28) give the following values of ag for the different 
orders of multipolarity : 


y-Ray energy E\ M1 E2 M2 E3 


676.5 kev. 0.0041 0.047 0.0111 0.105 0.0263 
1088.9 kev. 0.00173 0.0148 0.00453 0.030 0.0097 


It is seen that the measured value of ag for the 1088.9 kev. transition is 
consistent only with #2 radiation, and since the ground state of Hg! is 
almost certainly of even parity and zero spin, it follows that the state in 
Hg'®* at 1088.9 kev. has even parity and spin 2. Because the state at 411.77 
kev. also has even parity and spin 2, the 676.5 kev. transition involves no 
change in parity. It is therefore concluded from the measured value of ax 
that the radiation from the 676.5 kev. transition is a mixture of M1 and £2, 
with relative intensities of (32 + 6)% M1 and (68 + 6)% E2. This finding 
‘ will be discussed further in Section VIII. 

VII. SHAPE OF THE 1371 KEV. 6-SPECTRUM 

It has been shown in Section VI that the states in Hg!** at 412 and 1089 kev. 
are both of even parity and spin 2. The log fof of the 959 kev. 6-transition is 
7.4 and the spectrum has the allowed shape. Since the nuclear charge is high, 
these facts are consistent with a first forbidden 8-ray transition (18). The 
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log fot value of the 283 kev. 8-transition is 7.6 as determined from the measured 
y-ray intensities (Section V), so that it too is expected to be a first forbidden 
transition. The parity of the ground state of Au!® is accordingly expected to be 
odd. By assigning the odd proton to ad state and the odd neutron to a # state, 
the same conclusion is reached from nuclear shell structure considerations. 
The spin of the ground state of Au'®* must be 1, 2, or 3, all other possibilities 
being inconsistent with the spectral shapes and log fot values of the 959 kev. 
and 283 kev. 8-transitions. A spin of 1 would require that the 1371 kev. transi- 
tion be first forbidden, which is inconsistent with the forbidden shape found 
for the spectrum and with the log fo¢ value of 11.7. A spin of 2 would require 
that the transition be first forbidden, which is inconsistent with the log fot 
value, and it has also been found that the shape of the spectrum does not 
conform to the a-type (43) required by this assignment. A spin of 3 remains 
the only possibility consistent with all the experimental findings, and leads to 
the conclusion that the high energy 8-spectrum is third forbidden. 

The shape of the spectrum might therefore be expected to yield information 
as to the interaction form governing 8-decay, and it has been studied with this 
end in view. The interaction Hamiltonian is written as \sHs5 + AyHy + ArHr 
+ \4H,4 + ApHp, where the subscripts S, V, T, A, and P stand for scalar, 
vector, tensor, axial vector, and pseudoscalar. The \’s are coupling constants 
and the #7’s are as defined by Konopinski (14). The pseudoscalar term makes 
no appreciable contribution to a third order transition, and assuming the 
Fierz conditions (7), only the four interaction combinations VA, ST, SA, and 
VT are considered. The shape factors are taken from Greuling (10) and 
Pursey (25). The effect of the finite radius of the nucleus has been taken into 
account by an extension of the work of Rose and Holmes (29, 30) for Z = 83. 

The nuclear matrix elements are of the form Q;;,(a) defined by Greuling 
(10). The matrix elements which appear are those with a = r, 6 Xf, a, 6r, 
86 X r, Ba. If a and b are any two of these vectors, it is known (17) that the 
ratio Q;;x(a)/Q;;.(b) is independent of 7, 7, k and is real or pure imaginary. In 


fact 
Qijx(@)/Qinx(¥) = pt, QinlO X v)/Qin(t) = ot, Qij-(Ba)/Qij(a) = 7 
where p, og, and 7 are real. Also 


Qijx(6r) = —Qin(t) and Qijx(86 Xr) = —Qinx(é X r) 


when small components of the wave functions are neglected in comparison 
with the large components. 

Denoting the shape factor by C, which is proportional to (N/7?F)-(E,)—E)-?, 
the quantity C'” is found from the experimental data at eight energy values. 
These are plotted in Fig. 5. The calculated values of C'”? for each of the com- 
binations VA, ST, SA, and VT were then plotted, normalizing all curves to 
the experimental value of C!” at 1200 kev. This procedure required estimates 
of the pertinent matrix elements and these were made using a special nuclear 
mode! based on nuclear shell theories. Since the ground state of Au!%* has a 
spin of 3, it is indicated that the odd proton is in a d3,. state and the odd neu- 
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Fic. 5. The square root of the shape factor vs. electron kinetic energy. The circles represent 
the experimental data. The vertical bars represent standard deviations. Only the best fit for 
each of the combinations SA, VA, VT, ST is shown. For the various curves, the values of the 
parameters as defined in the text are: 

(a) o = —4/3; (c) o = —4/3, p= 4.94, r= 4/3; 
(b) o = —4/3, p = 4.5; (d) o = —4/3, p=19, += 4/3. 


tron in a p3,2 state. A reasonable configuration is (d3/2)* for the protons and 

(713/2)'2p3;2 for the neutrons. It is only necessary to consider the state with 

M; = 3 in calculating the matrix elements (38). This requires m values of 

3/2, 1/2, —1/2 for the protons and 3/2 for the neutron. The one-particle 

transition is therefore from a p3,2 (m = 3/2) state to a d3;2 (m = —3/2) state. 
The following values of the matrix element ratios were found: 


p ~ 9/2, o¢ = —4/3, rt = 4/3. 


The value for p is highly approximate and may well be in error by a factor of 
.two or more. The values of « and 7 do not depend on the radial parts of the 
one-particle wave functions. 

Using the above values of ¢ and 7 and various values of p in the neighborhood 
of 9/2, attempts were made to fit the theoretical results to the experimental 
data. The other parameter in each case, \7/As, Aa/Av, Aa/As, OF Ar/Av, Was 
adjusted to give the best possible fit. From the families of curves obtained only 
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those giving the best agreement with experiment are shown in Fig. 5, one for 
each of the interaction combinations VA, ST, SA, and VT. A brief discussion 
of these results is given in Section VIII. 
VIII. DISCUSSION 

The measured value of ax for the 1089 kev. y-ray transition is in excellent 
agreement with the theoretical value for E2 radiation, and since this transition 
leads to the ground state, the radiation is identified as pure £2 with certainty. 
It follows that the 1089 kev. state in Hg'®* has spin 2 and even parity. The 
value of ag for the 677 kev. y-ray, identifying this radiation as a mixture of 
M1 and £2, is also consistent with the assignment of spin 2 and even parity to 
the 1089 kev. state. Recently, Schrader et a/. (34) have investigated the angular 
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(0.98 0.05)% 


B- 
959.0+ 2.5kev i088.9 kev 


5£0.5)x10*% 1088.9£0.9 kev 
7 (0.16 £0.01) % 


411.77 kev 


0 5 i 0 


Fic. 6. Proposed disintegration scheme of Au'®%’. All energy measurements are based on 
the energy of the 411.77 + 0.036 kev. y-ray measured by Muller et a/. (22). 


correlation between the 677 kev. and +12 kev. y-rays and have also found the 
677 kev. y-ray to be a mixture of 1/1 and £2. They have found the mixture to 
be (40 + 5)% M1 and (60 + 5)% E2 as compared with the results (32 + 6)% 
and (68 + 6)°% obtained in this investigation. Schiff and Metzger (33) have 
also studied the angular correlation and have estimated the mixture to be 
(30-50)% M1 and (70-50)°% E2 radiation. The spin and parity assignments 
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found from these angular correlation experiments are in agreement with those 
given above. 

The ratios ax/a, found in Section IV provide a further means of checking 
the multipolarity of the y-rays. The result of 6.3 + 0.5 for the 1089 kev. y-ray 
is in fair agreement with the extrapolated empirical curve of Goldhaber and 
Sunyar (8) for E2 radiation which predicts a value of about 6.8. The result 
of 5.7 + 0.5 for the 677 kev. y-ray is however not in agreement with the 
extrapolated curves for 1 and E2 which predict a value of ~7.1 for 32% M1 
and 68% E2 mixture. A mixture of 40% M1 and 60% E2 as found in the angu- 
lar correlation experiments (33, 34) leads to a value of about 7.2. It is to be 
noted, however, that the extrapolation of the curves is rather uncertain. 

The investigation of the shape of the 1371 kev. 6-transition suffers from the 
severe limitation that only that portion of the spectrum above 959 kev. is 
available for study, owing to the presence of the intense 959 kev. spectrum. It 
can only be said that the combination SA is in definite disagreement with 
experiment. For the other three combinations, VA, V7, and ST, fairly good 
fits to the experimental data can be obtained with reasonable values of the 
parameters. As indicated by other experiments (23, 31), the combination VA 
may be rejected since it contains the axial vector rather than the tensor form 
for the Gamow-Teller part of the interaction, but the present data nevertheless 
do not provide a clear-cut preference for either one of the remaining two inter- 
action combinations. 
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ON THE NUSSELT MODULUS AS A FUNCTION OF FRICTION 
FACTOR, REYNOLDS AND PRANDTL MODULI, FOR HEAT 
TRANSFER TO A FLUID FLOWING TURBULENTLY 
THROUGH A TUBE OF CIRCULAR SECTION! 


By A. W. Marris? 


ABSTRACT 


Theoretical expressions for the Nusselt modulus and radial temperature 
distribution in terms of Reynolds number, Prandtl number, the ratio & of the 
eddy diffusivities for heat and momentum, and the friction factor are developed 
for the case of heat transfer from a tébe of circular section to a fluid flowing 
turbulently through it. The results, appropriate to fluids with Prandtl numbers 
greater than 1/2k, are derived on the assumption of the constancy of fluid 
properties and require a linear longitudinal temperature gradient in the tube 
itself. They enable the effect of the insulating ‘“‘laminar’’ and “‘buffer’’ layers 
adjacent to the tube wall to be demonstrated at various Prandtl numbers and 
are appropriate for correlation with experiment for the investigation of the 
constancy or otherwise of the ratio of the eddy diffusivities for heat and mo- 
mentum implicit in the momentum transfer analogy. 

A similar result is worked out on the vorticity transfer analogy for the par- 
ticular case of Prandtl number equal to 1/k’ and is found to be invalid for such 
a configuration. Here k’ is the ratio of the eddy diffusivities for heat and vorticity. 


1. SUMMARY OF NOMENCLATURE 


= mass density of fluid at point P in fluid. 
mean pressure at P. 


axial distance of P. 
radial distance of P. 
radius of tube. 
velocity of fluid passing over P. 
viscosity of fluid at P. 
n/p = kinematic viscosity of fluid at P. 
temperature at point P. 
specific heat of fluid at constant pressure. 
molecular conductivity of fluid. 
s/cp = thermal diffusivity of fluid. 
éy = eddy diffusivity for momentum defined by 
€,ou/dr = — U;u,, u, and u, being the turbulent velocity fluctuations, or by 
ty = ly°du/dr, ly being the “mixing length’’ for momentum transfer. ¢~ may 
be regarded as the effective increase in kinematic viscosity due to turbulence. 
é, = eddy diffusivity for heat defined by e,dt/dr = — u,t’, U being 
the turbulent temperature fluctuation, or by ey = /,’du/dr, ly being the 
“mixing length” for heat transfer. Again ey may be regarded as the effective 
‘increase in thermal diffusivity due to turbulence. 
Rk = én/€y. 
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ey = eddy diffusivity for vorticity, the virtual viscosity due to turbu- 


1a Ou = 
lence on vorticity transfer theory, and may be defined by ss (- a = u,t’ 
r r 


where ’ is the turbulent vorticity fluctuation. 
k' = en/ey. 
h_ = heat transfer coefficient for inside surface of tube wall. 
@w = heat flow per unit time normal to wall surface: 
A, = area of wall surface across which g » flows. 
A = 2arz = area of cylindrical surface in fluid across which heat flows. 
tw = temperature of wall surface. 


‘“‘mixed”’ temperature of fluid defined by tm =f ‘unadr [ f “rdr. 
0 0 


. e 
mean velocity of flow, defined by um = f urdr | f rdr. 
0 0 


lw aP 


Fanning friction factor = — ; 
B , PUm? Oz 


10P 2r, 


shear between wall and fluid defined by — es 
Pr w 


Reynolds’ modulus = 27, 1m/v. 


; ; ; Tw u 
nondimensional velocity = u ac rae 
p Um 


= r/r» = nondimensional distance from axis. 
RevV/(f/8) . 
Q(1 — x) = nondimensional distance from tube wall. 
a value of y stipulating extremity of laminar layer. 
B value of y stipulating extremity of buffer layer. 
¥ a[l — log (a/B)]. 
= value of y stipulating extremity of laminar layer when buffer 
layer is omitted. 
K, L, M, N = constants in equations for velocity distributions. 
Nu = Nusselt modulus = 2hr,,/s. 
Pr Prandtl modulus = cpv/s = cu/s. 
P = kPr. 
P’ = k'Pr. 
2. INTRODUCTION 


Much of the recent theoretical work on the problems of heat transfer by a 
turbulently flowing stream has been based on the assumption of a formal 
analogy between heat transfer and momentum transfer. In Reynolds’ modifi- 
cation for turbulent flow of the Navier-Stokes equations there arises a quantity 
ey, the eddy diffusivity for momentum, an additional “turbulent” viscosity; 
likewise in the equations for heat transfer (equations of similar form), there 
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arises ¢,, the eddy diffusivity for heat, an addition to the molecular thermal 
diffusivity occasioned by turbulence. The solution of the heat transfer problem 
on the basis of analogy necessarily requires a relationship to be assumed be- 
tween the quantities ey and ey. If one sets ey/ey = & then the constancy of 
the quantity k is a measurement of the validity of the analogy. It is re- 
membered that momentum is a vector quantity whereas heat is a scalar. 

Previous authors employing the momentum transfer analogy, von Kar- 
man (6), Boelter (1), Martinelli (1, 6), Lyon (5), have each set k = 1, however 
the theoretical investigation of H. Reichardt (8) of turbulent heat transfer 
from a flat plate and in a tube indicated that if k was computed from the 
measured temperature profiles it would have some value between 1.4 and 1.5. 
The temperature measurements of Sheppard (9) indicated a value 1.6 for this 
quantity; meanwhile measurements of heat transfer in wakes yield a value 2. 
However, it is known that the momentum transfer theory of turbulence does 
not predict the correct velocity distribution in wakes; but that the vorticity 
transfer theory propounded by G. |. Taylor (11, 12) yields satisfactory results. 
So the question arises as to whether vorticity transfer analogy considerations 
would yield better consistency. It is known however that the vorticity transfer 
theory breaks down in the neighborhood of solid boundaries. 

In this paper, concerned with the heat transfer from a tube to a fluid flowing 
turbulently through the tube under the prescribed condition that at any point 
of the fluid the longitudinal temperature gradient is independent of the radial 
position of the point, i.e. d¢/dz = dt,~/dz (see nomenclature), the Nusselt 
modulus and temperature distribution are expressed in terms of the friction 
factor f, a quantity Q equal to Rev/(f/8), Re being the Reynolds number, 
and P, the product of the Prandtl number and the ratio k. Expressions are 
obtained first on the momentum transfer analogy for all P greater than 3, 
and then on the vorticity transfer analogy for the particular case P’ = 1. 
Here P’ = k’Pr where ’ is the ratio of the eddy diffusivities for heat and 
vorticity. 

It is found that the requisite integrals admit exact evaluation over the 
complete velocity distribution consisting of a “laminar” adjacent-to-wall 
region, a semiturbulent ‘‘buffer-layer’’, and a fully turbulent core. Martinelli’s 
device (6) of taking the velocity at a given point as equal to the mean flow 
velocity as a first approximation is avoided as is the series approximation to 
the logarithmic velocity distribution employed by Reichardt (8). 

It appears that the vorticity transfer analogy cannot give a reasonable 
result for the particular configuration considered so that attention is refocused 
on the momentum transfer analogy. 

Assuming constancy of fluid property values throughout the tube, the 
results are appropriate for correlation with experimental results for infinitely 
low heat transfer rates, for example with the work of Eagle and Ferguson for 
water (2), so that the constancy or otherwise of the ratio e¢y/e¢y may be further 
studied. 

The results yield quantitative information concerning the heat “barrier” 
effect of von Karman’s “buffer-layer’’ between the laminar surface layer 
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and the turbulent core, at various values of the Reynolds and Prandtl moduli, 
and indicate the magnitudes of errors likely to be introduced by uncertainties 
in the ‘‘near-to-wall’”’ velocity distribution. 

The temperature distribution formula predicts the zero temperature gradient 
at the tube axis as found experimentally by Sheppard (9). 

3. GENERAL INTEGRALS FOR (tw — ts)/(tw — tm) AND Nuf DERIVED 
FROM THE MOMENTUM TRANSFER ANALOGY 

On the statistical theory of turbulence of Reynolds for which momentum 
and temperature are transferable properties the respective hydrodynamical 
equations representing flow of fluid through a pipe of circular section at a large 
distance from the end, and for heat transfer through the turbulent fluid can 
be written 


fio Ou 

7. ar Ee = v) ae ’ 
1a 

rT 


ot 
2 ee A) ie 
2] 8 | we +0) # 
Equations [1] and [2] indicate the nature of the eddy diffusivities ey and é, 
as ‘‘turbulent additions’”’ to the molecular thermal diffusivity and kinematic 
viscosity respectively. 

Heat is flowing radially into the fluid from the internal wall of the tube 
and is taken away as the heated fluid moves on down the tube. The heat 
transfer coefficient h is defined by the equation 

Vw 


8) "= Alle — tm) 


so that 


Zhrw _ ww 

s  AwS(tw — tm)” 

Following Lyon (5) the ratio (t~ — t)/(tw — tm) and the Nusselt modulus 
are expressed as iterated integrals. The suffixes 1, 2, 3 refer to current positions 
in the fluid. 

A heat balance on an annulus of fluid of radius 7 gives 


ot 0q1 
5 ee: 
[5] u; Ai pc a dr; 


[4] Nu = 


so that using [2] 


he ig A ae | se) 
Aicp on. On; ri(ex + a) or; 
which, since gq; = 0 at 7; = 0, integrates to 


A SM Ne ot 
[6] - Aicp Aywcpx (cx + a) Or, 


Again by [5] 


at (°” 
= ¢p af. u, A, dr, 
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ot - ae” 
= 69 — u, Ai dr, = Com > A,dn 
Iw p 02 0 on Oz 0 
so that 


Ts RTo els 

= Uy ci Uy 

[7] q2 = aw f — Ay, an | A, dr; = a0 | — x; dx}. 
0 Um 0 0 Um 


By [6] and [7] 


ot qu Gane 
aaa. sinte _—eomome <7 Ss dx, 
Ore colex ae a) Xs se 0 Um 


so that 


or q vl a Za u 
alo , 1 
[8] by — ts = “ J eee aS -J —— wn dx, dx, 
SAy Zs (ex + a) 0 Um 
al z 
4r - Us 1 * uy 
ty — bn = =o Je —~ X3 te ee te oer dx, dx» dx3, 
A 7) vi Um Vr, cp (€x 4 7) Xo 0 Um 


: olf E Us nf — f° ee | 
[10] Nu = 9 bes X3 eres Bes x1 dx; dx. dx3 |. 


Expressing the equations for steady motion of fluid as 


[11] T= XTw, 
[12] V(t w/p) = Um (f/2) , 


nondimensional velocity and distance from tube wall are defined thus: 


+ u 2 
[13] y= #4/?, 


[14] y = Qi — x) 
[15] where Q = Rev(f/8) . 


Equation [1], which can be written 


[16] Se, 1 OT eu + v) ie 


Pl wv r or 


vields, on integration and transformation, 


~ eu _ AS : - dU 
[17] (0 + dx - 


v dx 
Writing 
=f —1 
aac: €x# -_ €m 
[18] ta ta (« + :) (p v + ) 


where Pr = cu/s = v/a, the Prandtl modulus, and 
[19] P = kPr, 


<a 
20) ey ta a 
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Using [13] and [20] the basic equations [8], [9], and [10] now become 


tw — bs _ Lt a im sal 
[21] ok = nu4/f es Xe (1 P) dU /dx; ; U; x; dx; dx 


and 
8 . | _PQxs r rf 
[22] Nu f = if U3 x3 Ea Xe2 (1 — P) dU /dx; U,; X1 dx, dx dx3. 
The integrals [21] and [22] have to be evaluated for the appropriate radial 
velocity distribution U = U(x). 
Following the work of Nikuradse (7), Stanton and Pannell (10, 4), and 
von Karman, the velocity distribution will be considered in the general form 
U = U'(x) = Q(1 — x) for 1 > x > 1 — (e/Q), 
the ‘‘laminar’’, adjacent-to-wall region; 
[23] }U = U(x) = K+ L log Q(1 — x), for 1 — (a@/Q) >x > 1 — (6/Q), 
the ‘‘buffer-layer”’; 
U = U'"(x) = M+ N log Q(1 — x) for 1 — (8/Q) >x > 0, 
the turbulent core; 


where between the several constants a, 8, K, L, M, and N the boundary 
relations 
[24] L=a; K =a(1 —loga);} M+N log B = afl — log (a/8)] = v 


hold. They express continuity of U and dU/dy at y = a, the boundary be- 
tween laminar and buffer layers, and continuity of U only at y = 8, the 
boundary between the buffer layer and the turbulent core. The relations [24] 
used in conjunction with experimental results of Nikuradse etc. yield the 
numerical values a = 5, 8 = 30, K = — 3.05, L=5, M = 5.5, N = 2.5, 
y¥ = 14. 


4. EVALUATION OF INTEGRALS 


Equations [21] and [22] may be written 


‘ ly — ts ae) f 
[25] ae nu4/ Is, 
[26] 1/Nuf = I 


v here 


1 —1 z 
so ee ca Pon | r 
[27] J X2 E P) dU /dx; ; U, x; dx; dx 


and 


1 
[28] I f U3 X3 Tz dx3. 
0 


Before proceeding to the general case the particular case of laminar flow 
throughout the tube section is considered. 
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u _ /2 6 
In this case eg = ey = O and U = 44/2 == "e (1 — x?) throughout the 


rangeO <x <1. 
Thus 


eo 2 ok © ae affine 
[29] n= V2 ff (x1 — x1) dx, dx2 = f Sener mae as 


and 


Pa 2 


1-8 f (xa — xs°)(3 — 4x3? + xs") 
7 f J: 16 48f 


[30] 


so that 


[31] ans * = 4.364 


and 


(32) fo 8 — 4c* + 2'). 
L- 1 

The general results valid for P > 1/2{1 — (8/Q)] > 1/2[1 — (@/Q)] > 3 
will now be derived. 

The results will involve convergent infinite series of the form >> F,(PQ)/Q” 
and the question arises as to the lowest power of 1/Q that can be neglected 
without causing appreciable error. For the simpler case of buffer layer omitted, 
a result valid over the restricted range } « P < 1 was worked to O(1/Q*). 
It is found that in this range, that is for P « O(Q), the error involved in 
neglecting terms O(1/Q?) was about 1% at Q = 300 corresponding to Reynolds’ 
number 104 and about 0.2% at Reynolds’ number 5 X 104. 

For the sake of avoiding undue arithmetical complication, the derivation 
for all three layers is restricted to O(1/Q). Inaccuracies may exist in the region 
P = O(Q) owing to the fact that certain terms of the form P/Q are neglected, 
but it is considered that uncertainties in the near-to-wall velocity distribution 
very effective at the higher Prandtl numbers will have vitiated the results 
anyhow before these values are reached. 

The problem is that of evaluating the integrals [27] and [28] for the velocity 
distribution given by the equations [23] and [24]. 


Zs 
Writing I2(x2) =f Ux,dx,, integrating, and using the relations [24], the 
0 


following results are obtained. 
(a) For0 < x2 <1 -— (8/Q) 


Io(x2) = Io!" (x2) = J [M + N log Q + N log (1 — x1)] x1 dx 


2y-N WN N N 
= (ae — — log 8) a 2 (1 — x2") log (1—x2). 


4 2% 0 2 





174 CANADIAN JOURNAL OF PHYSICS. VOL. 32 
(b) For 1 — (8/Q) <x. <1 — (a/Q) 
1—(8/Q) 
I2(x2) = Ie!" (x2) = J [M + N log Q + N log (1 — x1)] x1 dxy 
0 


i‘ J "IK + Llog Q + Llog (1 — =:)] x1 dx: 


1—(8 /Q) 
[34] = (fa — 4a) x2” — Jax2 — Sa(1 — x2") log (1 — x2) + E(Q) 


35) where BAR al denim v4 log ; +" a |. 


(c)- For 1 — (a/Q) < x2 <1 
»1—(8/Q) 


Filey 0 Eee J IM + log 0+ N log (1 — »,)) 21d, 


0 


Ts 


21—(a/Q) a 
+ | [K of i log Q + i log (1 = x1)] xy dx + | Q(1 <= x1) xX dx 
e v1 


1—(8/Q) —(a/Q) 
[36] = I Q(3x2” — 2x2") + F(Q) 


N,.., £ *—2(a—N 
[37] where F(Q) = — = log 2 = ee 
hg 


The function 
I; 


may now be evaluated. 
Again the range is divided into three intervals. 
(d) For 1 — (a/Q) <x3 <1 where dU/dxz = —Q 
*! 5Q(3x2" — 2x2") + F(Q) 


hin) =e) = f Pea pe 


9 she 2 
= Qa _ «(2 —- ws) a a Pet) log ws| 


20 -—83 
+ fees E> ae + GQ) |toe {1 — P(1 — x3)] — log x3] 
} 1—P 6P 
[38] 4 
where G(Q) = F(Q) + e 
. te 8 aN Ee -8e N (3) 
= 9 log 0 + 1 + 20 +0 oO): 
If P becomes greater than Q/a, then the term log [1 — P(1 — x3)] becomes 
undefined over a certain subinterval of the range 1 > x3; > 1 — (a@/Q) for 
which it is required. It will be seen below that to O(1/Q) the term involving 
T;'(x3) contributes nothing to the value of the Nusselt modulus. So if P > Q/a 
but < Q/@say where 6 < a, one can integrate over the subregion 1 — (@/Q) 
< x3 < 1 of the laminar layer instead of over the whole layer and to O(1/Q) 
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the result will be unaltered. This is clear also from the physical fact of relatively 
low heat transfer over regions of laminar flow. 


(e) For 1 — (8/Q) < x3 < 1 — (a@/Q) 
I3(x3) = T3'" (x3) = f x2" E “= ae - Ts dx 
po] =f" BlO)-+al 02 og (a/O) ae — bee— Y=) oe =H] 


xo[(1 — P) — (PQ/a) x2(1 — x2)] 
f° iQeare = 2s) + FO gy, 
e 1—(a/Q) Pxo"[1 + (1 — P)/Px:] 
by [34] and [36] and since dU /dx, = — Qxein the interval 1 — (a/Q) < x3 < 1 
and dU/dx, = — L/(1 — x2) = — a/(1 — xe) in the interval 1 — (8/Q) < x; 
<1 — (a/Q). 

Putting x3 = 1 — (a/Q) mn [38] it is found that 


*  4Q(8x." — 2x") + FQ) 
aa I wi Peell + 1 ~ P a: 


a N ] 
5 Q |-4 > log Q ns ae aN | oft - 
To evaluate 


: “/2 F(Q)+a[}(1—2 log (a/Q))x2"— 4x2—4(1—x2" ) log (1=x2)] 5 
J, xo{(1 — P) + (PQ/a) x2(1 — x2)] 


the denominator should be factorized, thus 


3 


X29 


r 
v2 (1 a P) + FO x, (1 —_ «| = Fe (A on Xo) (mu + X2) 


where \ = 1+ aa ne +O (2) ; 


-sigP+o(Q) 
be a pai O ay 
and the integral resolved into partial fractions. 

However, since P > 1/Q, it is clear that a close approximation to the 
answer will be obtained by writing the denominator x.[(1 — P) + (PQ/a)x2 
X (1 — x2)] as (PQ/a)x2?(1 — x2), a device which simplifies the analysis consider- 
ably. The outline of the more exact reduction is given in the appendix, and it 
is found that the result obtained by. the above device is correct to O(1/Q) 
except for two inconsequential terms. 

Writing the integral then in this form, one has 


a 1—(a@/Q) E(Q)+al[i a oni tom (a Os -3x2— 4 (1—x2" *) log. (1=*2)] 4, 
PQ. + = i Xo) 


‘ 1—(a@/Q) 3 ae ef a) 1—(@/Q) dx. 
ane ee = —2 log — ane 
po £) J hl +i dx2+ 1P0 1 log 0 i; is 
2 1—(a/Q) 1—(a/Q) 
a 1 I a 1 1 : 
~ 2P0 J, & *T- ‘J om — 386 f E ? 1] me Cane 


[41] = — it log ; a 4 x 
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x|(4- ~) + tog (1 — 2) — tog #1 + lo (1 #5) — oe 3 
oa IO . Q ee ’ Q 


¢ — 2 log 2)| tog (1 — - — log 


_a® [log (a/Q) _ log on (1 — 2) | 
* 5P0 i a” ee * 
where, since 


= const. = (1)=— 


oo 


[42] (x) = logxlog (1 — x) + p i » : z 


1 


x) 


one has 


= [1 — (a/Q)J" — xs" _ ee ( se =| 
43] » rm = log x; log (1 — x3) log & 0 log \ 1 0 
4 Gan)" = @/Q)"* 


n° 


(f) For0 <x3 <1 — (8/Q) 
al =i 
I3(x3) = T3'"" (x3) = § x2 la =p oe fe I, dx2 


a 1—(8/Q) 
a [4 Qy—N) — 3Nlog (8, ‘Q) x2" — 3Nx2— (1 =x" )log (1 — xe) 
[44] = J xd — P) + (PQ/N) m(l — x)] os 


- ie E(Q)+a[3[1 —2log (a/Q) x2” — }x2—}3(1—x2 “)log(1=x2)] 4 
1—(8/Q) xe[(1 cet od Me (PQ/a) ost. — x2)] 


1-(a/q) Pxe [1 + (1 — P)/Px2] 
by [33], [34], and [36] and since dU/dx2 = — Qxz in the interval 1 — (a/Q) 
< x3 < 1,dU/dx. = — L/(1 — x2) = — a/(1 — x2) in the interval 1 — (6/Q) 
< x3 <1 — (a/Q), and dU/dx, = — N/(1 — x2) in the interval 0 < x; 
<1 — (6/Q). 
Putting x3; = 1 — (8/Q) in [41] and reducing to O(1/Q), and adding the 
integral [40], the last two integrals of [44] have the value 


[45] H(PQ) = (a < eS ., x a £) i 4 o( * ) 


where y = a(l — log (@/8)). 


ia sQ(3x2" — 2 2x2") ee Xo 


If the exact analysis is carried out, one obtains the small additional term 
[4(2y — 3N) — EN log (8/Q)](a/PQ) log [P + (a/8)(1 — P)], see Appendix, 
giving the true value to O(1/Q) as 


os -|( 26 +5 log | P+2c i eqs me 5)o+o(g): 
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We thus have, approximating as before, 
I," (x3) = H(PQ) 


& Ss [4 (2y—N) — $Nlog (8/Q)|—3Nx2- — 4N(1—x2”) Jlog(1=a2)_ 
PQ x2"(1 — x2) 


eee oh 1-(8/Q) 
= H(PQ) + 556 nar 4 8) f i 


N? a1—(8/Q) f 1 ! 2 2 1—(8/Q) | 1 1 | 
2s 2PO J re +- 1 penvee - dx. — 2P0 4 a bt oa log (1—x2) dx» 


[47] = H(PQ) + ea Gis — J log ai (1 — xs) — log 8] 
N 18 (1=xs) __ log (8/ O14 4 U=G/or = x 
~ 2PQ X3 1—(8/ ‘Q) 2PO0 n° 

which, using the approximate value [45] for H(PQ) and oo gives 

N° | tog BY’ me “* (Qa — 4y + NON |g 8 a(2y — 3N) 


heak" aaa G Rg ae 


—y)(2y—-3N Nix’ 2y- 
—_ | eon@r-am) _ -- “|. Q a 2P0 2 ee —wN log s| log (1 = X3) 


- | iy 
[48] — 3P0 log. en xs) 


If the unapproximated form of H(PQ) is used, and the denominator of the 
first integral of [44] is factorized according to 


v.| ( ~ P)'+ 7e (I = x) | = re (§ — x2)(n + x2) 


and the integral evaluated by resolving into partial fractions, then, as in the 
appendix, the right-hand side of [48]: will contain the additional small term 


(27-6 -- * tog aa log | P oh atl - P)| — N log E + ey 7 : 


The substitution of equation [48] in 


BY ly — ts ae f 
[25] to = ug/L q; 


now gives the requisite temperature distribution in the turbulent core once 


Nu is determined. 
Proceeding now to the Nusselt modulus one has 


1 1—(8/Q) 
Nuf > [= J [M + N log Q + N log (1 — x3)] x3 Is’ (x3) dx 


2 i—(a/Q) 
+ J [K + Llog Q + L log (1 — x3)] x3 Is’’ (x3) dxs 


1—(8/Q) 


op of (1 - x3) X3 T;' (x3) dx 
1—(@/Q) 
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where = _7;'’’(x3) is the value of J3(x3) in the interval 0 < x3 < 1 — (8/Q), 
T;/'(x3) its value in the interval 1 — (8/Q) < x3 < 1 — (a@/Q), 
and I3'(x3) its value in the interval 1 — (a/Q) < x3 <1. 


By equation [47] 


+1—(8/Q) 
J [M + N log Q ++ N log (1 aa X3)] X3 T3'"" (x3) dx3 
0 

a 1—(8/Q) 


[50] = H(PQ) J [M +N log 0+ Nlog (1 — x:)) xa des 
N =. i. B 
* 2P0 (2 : 3 ’) 
1-(8/@) B 
x< J [M+ N log Q + N log (1 — x3)] x3 E (1 — x3) — log 0 dx3 


yy? rte) 


+ 359 


2PO [M+ N log Q + N log (1 — x3)] x3 


log (8/Q)_ _ nett ai] 
i — (6/0) xe m* 


N? re 1—(8/Q) 
“fb 2P0 , [M+ N log Q + N log (1 — x3)] x3 


Xx E x3 log (1 — x3) — log oe ¢ _ 


the last term following from the identity [43]. 
Substituting for H(PQ) from [45], integrating, and simplifying, one obtains 


»1—(8/Q) 
[51] J [M+ N log Q + N log (1 — x3)] x3 Is’ (xs) dx 


0 


Flees) +i [aps R= agate) 


—aN(2y — 3N) , Nila — y)(2y — 3N) — (y° — GN +9N 2) tog 2 
. | 4Q * 4PQ log 6 


ee 
} 0 

7N*y _ 71N" . zs 2y - any pee (4, ) 
+ [2 oes 4 ee 


If the more exact expression for J;’’’(x3) is used, then a term 


2y — 3N Nv NG — P 
= roe 9 log 8) log (p+2 +2 - p)) — N log (1 + — =e BP ) | 
2 1— ae 
oo: [M + N log Q + N log (1 — x3)] x3 dx 
ae ae 
x | log (p + < 


must be added. 
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Using equations [40] and [41] for J;’’(x3), and equations [38] for J;’(x3), it 
is found that neither of the last two integrals of [49] yields terms to O(1/Q); 
that is, to this order, in the third integration, the laminar and buffer layers 
contribute nothing to the value of the Nusselt modulus. 

It follows finally that 
[53] 

1 i ( 8) 1 [ aN? 
=a ge oa =e eas 
iw” “mero T4la" 2PO 

—aN(2y —3N) , N[(a — y)(2v — 3N) — (7° — 6yN 
+ —* ee Vila — y)(y — 3N) — (vy — 6yN 


2y ee 
‘ee 


7N? 71N 3 N 3 
- k ae teen 


+ €E a O (2) ’ 
where 


N 2y—3N |’ Ni-P I 
.-|} log 6 =f 1 | | log] P+ —P) | log] 1+ oP) PO 


and is small. 
5. ATTEMPT TO OBTAIN SIMILAR RESULTS ON THE BASIS OF 
THE VORTICITY TRANSFER ANALOGY 

Before analyzing the above results obtained on the basis of the momentum 
transfer analogy, an attempt will be made to obtain corresponding results on 
the ‘‘vorticity transfer’ theory. 

The vorticity transfer theory was propounded by G. I. Taylor (11, 12) who 
pointed out that a statistical theory of turbulence with momentum conserved 
as the transferable property, as originally proposed by Prandtl, required that 
the fluctuating turbulent pressures do not influence the mean transport of 
momentum, and who further proved that this condition only holds when 
momentum in a given direction z say is transferred in a plane xy by turbulent 
motion in which lines of particles once parallel to the motion remain so 
throughout the motion. He suggested, having shown that the momentum 
transfer theory cannot hold for purely two-dimensional motion, that if the 
whole motion is confined to a plane, xy say, then by Kelvin’s circulation 
theorem, vorticity in the direction z is conserved and is the transferable 
property. The theory was subsequently generalized to three dimensions, and 
a modified vorticity transfer theory set up, one of the main assumptions of 
which is. that components of vorticity are transferable in the sense that heat 
is transferable (3). 

The modified vorticity transfer theory yields a turbulent velocity distri- 
bution for pipe flow in as good agreement with experiment as that derived 
from the momentum transfer theory. The theory seems to hold for wakes 





180 CANADIAN JOURNAL OF PHYSICS. VOL. 32 


(where the turbulence might well be approximately two-dimensional as in the 
case of a von Karman vortex street) but it appears to break down in the 


neighborhood of solid boundaries. 
The equation representing fluid flow through the pipe is now written 


oS P= r+ ot 2 (r%) 


so that in analogy with equation [17] one now has 
ey _ au) / 4. a 
vo {20 + 4 + oe | «dt I} dx L* dx 


. Sages he oer 
[56] até Pini) +1 
= a _ 2P’ Pie | ; 
or " la ) ae at hale 
so that one now has 
ty ex T i 
[58] os = Nu 4/4 J3, 
[59] 1/Nuf = J, 
where 
1 e 2P’Qx , ¢ 
pe ee ok IN ce ae ee T 
[60] J3 i Xo a P dx» dU /dxs}/ /dx. 5 Ll 1%1 dx, dx» 
and 
al 
[61] J J U3 X3 J3 dx3. 
0 


Owing to the complicated nature of the second integration only the case 
P’ = 1 will be considered. Also the buffer layer will be omitted. 
The integrals J; and J have thus to be evaluated for 
(U = U(x) = QU — x) for 1 > x > 1 — (8/Q), 


[62] lu = U(x) = M+ N log QUI — x) for 1 — (8/Q) > x > 0 
where from continuity of velocity at y = Q(1 — x) = 6, 
[63] 6=M+WNlogé. 
Writing Jo(x2) = t Ux,dx,; , one has: 
(a) For0 < x. < 1 — (8/Q) 
Paley ac Teo . ‘(aN 0+ Vie ~ ated 


—N N N 
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(b) For 1 — (6/Q) <x. <1 
1—(6/Q) 
Jo(x2) = Je'(x2) = J [M+ N log Q + N log (1 — x1)] x1 dx; 


+ of (x1 — x1”) dxy 
1—(6/Q) 


[65] = §Q(3x." — 2x2*°) + F(Q) 










where 





3 __N, 5 Q, 2%-3N , #& —2(6-N)s (2) 
[66] F(Q) = > log 0 6 = i +> 20 +O oO) 
[67] For P’ = 1; Js (x3) == 20 Lf x (25 + X2 ae ) Je 2(X2 ») dx» 






so that: 








5 dU dU 
ere i = eas L where © , + xs os —-@Q 
al ’ 
, * Jo'dxy 
J3(x3) = Js’ (x3) = a 
al ~l nee 
= Sse — 2x2) dx_ + a4 ore by [65] 
. ; 1 F \ 
[68] - e (x3| — 3x3) + 2 (2 ae F(0)) + me 





and 













, oe ee 
(d) forO <x, <1 - 0 where dx: + aR (1 — x2)’ 
e Wr ee dx: a (x2) 
[69] Js(xs) = Js" (xs) = ae x2(1 — x2)" J. (/q) 2x2" 7 






_N 2 -N_N, #) "gee dx, N* Pie 
=20\ 4 "ol, G-ar 4. niil-—ar 


N* ¢-/® (1 + x2) log (L — x2) dx (a -3N N =) A ) 
~ 40, xe'(1 — x2) TX” 4 2 oe Ce 






- 






by [64] and by putting x; = 1 — (6/Q) in [68] 


1/22-3N NWN, 6\f6 .N\.N* 6 \’ 
. i( ls ale . 5) +H (lo 2 + (og 0) | 







[70] N ee 2) 2 
= 20 4 2 8 Q 1 — as 
- | be X3 ' flog (1 x3)] + 22 n- 3 + O @ é 





Equation [70] substituted in equation [58] gives the temperature distri- 
bution throughout the turbulent core once Nu is determined. 
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The Nusselt modulus is now given by 
1 


1—(8/Q) 
J= f U3! x3J3"dx3 + U3'x3J3'dx3 
0 


1-—(8/Q) 


Sears 
Nuf os 


2 1—(6/Q) 


Us!'x3J3"dx3 + O (2) 


re a)(5+ 5) + [Ce 6) + (me 5) 75] 


1—(8/@) 
x f [M+ N log Q + N log (1 — x3)] x3 dx 
0 


5 +1—(6/Q) x: 

: J [M + Nlog Q +N log(1 — x3)];—**—dvs 
0 — wa 

a 1—(6/Q) 


[M+ N log Q + N log (1 — x3)] log (1 — x3) dx; 
3)? (8 v) 
- 5) (3 ae 


P oD 73 
N og ste-®) iad (: _< 1 ha 
2 Q 40 Q 4 r2n/2Q 
In equations [70] and [71] the results are vitiated by the respective terms 
[4(26 — 3N) — 3N log (6/Q)]N/5 and [}(26 —3N) — 3N log (6/Q)/?N/6 
which are not O(1/Q). For any sort of a result at all these must vanish. 
N being fixed at about 2.5 from the experimentally well defined velocity 
distribution for fully developed turbulence, the only possibility is to alter 6. 
This could reasonably be done since the velocity distribution adjacent to the 
wall is uncertain. 
Of the two possibilities however, it is seen 
(a) that 3(26 — 3N) — 3N log (6/Q) = 0, giving the impossible result Nuf 
= Q/11.34, requires that logio(Q/6) = 0.6514 — 0.17376, which is equally 
impossible, 
(b) that (6/Q) + (N/6) = 0 requires 6 to be imaginary, which is also im- 
possible. 
It thus seems that vorticity transfer will not give a reasonable answer to 
this problem, at least until definite modifications are made. 


6. FORMULAE OBTAINED ON MOMENTUM TRANSFER ANALOGY: 


NUMERICAL RESULTS 
Inserting the numerical values a = 5, B = 30, y = 14, N = 2.5, equation 


[53] for the Nusselt modulus with all three layers becomes 


1 _ _ 47.69 0) a2 sae) a 
Nuf PQ (1ogi T\"9 + pO logie 6 


146.9 ah) 30: , 180.3. 1783 
aa ~eaoneees ——— } ] — an mm ole 
ogee wo ge ae ee 
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TABLE I 





a a az 


.67152 . 5366 . 12436 16 . 50672 —5.681818 

- 38037 ; a | . }. 27492 —4.80769 
.05109 , } 22, -77146 —4,16667 
.14096 | 6.91162 | 19. 23 | .03376 —3.67647 
5.39111 | 50 | 3. 56836 3.57865 —3.12500 
5.59984 | -646 » 5 24 38580 —2.08333 
- 96282 .f . 28933 —1.56250 
3.43836 | —32. | 142 9.14466 —0.78125 
.13861 | —27.6 76 5.09644 —0.52083 
38.13233 | —24.812 2. 5b .57234 —0.39063 
3.17014 | —23. 65 | 3.65787 —0.31250 
} —s 2.61276 —0. 22321 


Thee OOOO yy 


OER ae eS 


SSssssusensy 


moo 


TABLE II 
JES OF Q, RE, AND f, AS TAKEN FROM VON KARMAN’S LAW: 
1/f = (4.0 logio Q + 1.4062)? 


9,600.73 | 127.939 0.007816 
44,997 .8 | 186.295 | 0.005368 
129,948: 234.540 | 0.004264 
178,918 | 250.095 0.003998 
229,130 | 262.507 | 0.003809 
280,331 272.858 | 0.003665 
332,355 | 281.790 0.003549 
385,082 289.633 | 0.003453 
438,426 | 296.6837 | 0.003371 
492,317 | 302.976 0.003301 
1,052,749 | =: 346.343 0.002887 


Q( =Rev/(f/8)) Re \/f 





TABLE III 
NUSSELT MODULUS 


0.65 0.75 0.85 | ‘ 1.05 | 1.10 1.50 


Re 


| 
So 


: 


~ 9,600.73 i ee 5} 29.2176) 31.1608) 32.9671) 35.4647) 36.2482) 37.0078) 42. 4263) 
44,997 .8 1.9849) 90.1662) 97.3867) 104.279 | 113.833 | 116.838 | 119.755 | 140.667 | 162. 

| 

| 

| 


8 


129,948 .375 | 205.378 | 223.958 | 241.341 | 265.553 | 273.192 | 280.623 | 334.183 | 389. 
178,918 | 237.606 | 263.869 | 288.301 | 311.195 | 343.125 | 
288.344 | 320.768 | 350.969 | 379.297 | 418. | 431.278 | 443.532 | 531.563 
376.430 | 412.331 | 446.034 | 493 
431. 143 | 472.690 | :725 | 566. 583.601 | 600.455 | 726.194 | 850.36 
532.197 | 576.517 | 638.54: 
590.970 | 640.545 | 709.998 | 
| 528. | 590. 649.104 | 703.901 | 780.667 | 805.012 | 828.729 |1000.97 | 1182.; 
032749 5.523 |1094.95 |1206.99 [1312.77 [1461.46 |1508. |1554.82 |1891.00 | 2247. 


— 


623.5 


OOM Gre Co 


| ui 


Re ce hate 
4.00 6.00 8.00 | 10. 2.00 | 14. 100 
9,600.73 | 62.7278] 71.7416, 77.9940) 82.6479| 86.2 27 565) 
44.997.8 | 220.740 | 257.125 | 282.629 | 301.740 | 316.703 | 328.793 | 
129,948 543.662 | 640.785 | 709.478 | 761.254 | 801.961 | 834.947 |1041.5667|1105.2780| 1170. 5929 
178,918 | 
229,130 | 879.892 [1043.19 |1159.28 [1247.07 [1316.25 |1372 |1725.8842|1834.5586| 1946.9480 
280,331 
| 332,355 1208.82 1438.49 |1602.33 [1726.49 /1824.49 |1904.15 | 2406.8158]2561.9601] 2722.0826 
| 385,082 
| 438,426 | | 
1693 . 2¢ 22.81 |2258.67 [2437.85 |2579.51 2694. 3425.005 |3650.797 | 3883.859 
| 
' 


89 i 2234) 112 7315| 118. 4401 


| 492,317 | 
|1,052,749 13262. 25.99 .: 38.46 [5057.71 |5293. 





CANADIAN JOURNAL OF PHYSICS. VOL. 32 


TABLE IV 


VALUES OF ANU = NUSSELT MODULUS CALCULATED WITHOUT BUFFER LAYER MINUS 
N USSELT MODUL US CALCULATED WITH ALL THREE LAYERS 


kPr 


0.75 0.85 1.00 1.05 


2.0345 | 1.6468 | 1.1865 | 0.3964 | 0.1154 
"8440 | 4.193 3.097 0.296 : 
037 8.522 3387 “627 ‘043 | —16675 
2/410 | 10.600 ‘970 
‘690 | 12.565 ‘473 
3860 | 14.456 "921 
‘979 | 16.309 2'346 
.999 18.071 3.696 
2.976 19.936 5.015 
24.923 | 21.505 3 
2°78 37.19 





io 


9, 600. 73 | 

44,997 .8 
129,948 
178,918 


eo 


| —25.446 
—37.428 


ON oe 


—45. 406 
—81.46 


Com RO Whe 
to y bo 





6.00 | 8.00 


~ 93.1329|— 28,4183] — 
79.703 |— 99.931 |— 
2 |— 245.471 | — 286825 | — 320.: 


| 


.518 | — 309.762 |— 396.526 |— 465.504 |— 521. 
— 640.48 « 719.2 


492,317 ' 371.61 |— 588.86 |— 760.72 |— 898.57 |-1011.33 | —1105. 26 
052,749 46 | —696.97 | -1120.66 | ~ 1460. 62 |-1735.63 |-1961.91 | —2151.2: 


543.94 


| 
| 

| . | 

| 

| | 
| 


39.063 | — 


where « is small and given by 


_ [5.1 — 2.9 logio (30/Q)]* 
ia PQ 


x {is logio [P + 0.17(1 — P)] — 5.8 logio E + 0.083 ( =2)0 


Numerical values for Nu are given in Table III. 
The corresponding equation for buffer layer omitted, obtained by substi- 
tuting a = B = 11.5, y = 0, N = 2.5 in equation [53], becomes 


1 _ _ 47.69 ay a= was) nay 
Nuf * PO (10g: 0 + 0 + PO logio 0 


256.5 25.2 11.5 172.7 23.48 
[530] a > " eg! o +t = _ “PO. + € 


where e. is small and given by 


(7.75 — 5.8 logiy (11.5/Q)7° ; | uP ¢ —P 
PO 1.4 logio} 1 + 0.22 P 


Equations [53a] and [536] enable ANu, the effect of the buffer layer and in 
fact the order of the error caused through uncertainty in the adjacent-to-wall 
velocity distribution, to be investigated. In particular it is found that ANu = 0 
for P = 1.07. See Table IV. 

When P and Q are large such that P/Q tends to zero, equations [53a] and 
[536] take the respective forms 
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1 aot aay _ 146.9 1,30 __ 130.3 
re 


Ny @ V8" 


1 _ 95.27 (10g nay _ 256.5) U5, 172.7 
Nuf Q Q Q Q Q 
Putting Q = 2 X 10‘ corresponding to Re ~ 10%, it is found that at this value 
the heat barrier introduced by omitting the buffer layer more than halves 
the Nusselt modulus. 
For the temperature distribution equations [25] and [48] give 
Q tw — ts 


‘ 30 30\° 
NuV(/2) te— ty neg Fe (tom = 


n 


[72] + (« + a5 logio - logio (I — X3) + a6 logie (1 — xs) + az), ~ 
Q X3 1 ” 


where a), .. . , @7 are functions of P whose values are given in Table I. 
Writing 
@ n ez 
% log, (1 — x) dx 
a) « Bea J log. (1 — x) dx 
1 ” 0 x 


it may be shown that 


e(1) = =, 


9 rT = - . os € - =) | 
o(2 sin a | og. (2 sin 10 : 

pe 3 : = x a! | (: : 2yy 
of (2 sin =) ] 5 log, \ 2 sin io) |° 


7 2 
6(3) = 75 — allog. 2], 


0(0) = 0. 
Values of (tw — ts)/(tw — tm) are given in Table V. 
Equations [25] and [48] prophesy a zero temperature gradient on the tube 
axis. 
Thus in the region 0 < x3; < 1 — (8/Q) 


A=). we ae V5 | ee eae | 
OXaLty — ted a. 2Lx.[(1—P) + (PQ/N) (1—x2)x2] 


= ~ nug/d 


x{ er — N) — 3N log(B/Q) x2 + 3N x2 log(1—x2)—3N[1 + log(1 <2) 2] 
(1 — P) + (PQ/N) x2(1 — x2) 
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TABLE V 

C tw — ts 

VALUES OF ———;— 

tw —tm 
tw = wall temperature, tm = mean temperature of fluid, ¢; = local temperature at x3 = 13/ro. 
bP is aPr.) 








1.50 5 2.00 
| 


Al 


z wm rirem 0 





| | } | 
300 1.31072 2 | 1.22409 | 1.17910 .15211 | 1.10014 | 1.07681 | 1.062% .05367 
1,150 | 1. 24342 . 22101 1.18563 | | | | 
5,000 | 1.19467 | 18215 | 1.15479 | 1.13039 | 1.12069 | 1.07953 1.06208 .05112 | 1.04348 
10,000 1.17739 .16781 | 1.14321 | 1.12179 | 10734 | 1.07555 | 1.05923 .04154 
20,000 1.16257 .15527 | 1.13297 | 1.11396 . 10099 1.07180 1.05650 1.04667 .03969 


Il. x3 =r/ro = [2 sin (r/10)]? = 0.38197 


.04783 | 1.03935 | 1.03365 
03801 03106 | 1.02622 
| 


300 1.19129 17033 | 1.13890 | 1.11057 .09402 06215 | | 
1,150 | 1.15178 14117 | 1.11564 | | 

5,000 | 1.12210 | 1.11737 | 1.09681 | 1.08124 .07112 | 1.04908 | 
1 02946 | 1.02481 
1 


10,000 11401 10839 1.08961 | 1.07597 | 06675 | .04648 | 1.03608 
.02790 .02343 


20,000 10210 | 10042 | 1.08316 | 1.07093 | .06268 | .04399 .03419 


| 
| 
| 


III. x3 =r ro = 0.50 


paeeeaion aiensagal aeaepaaeaiaiad = . 


300 | 1.14010 | 12719 | 1.10141 | 08120 .06912 | 1.04586 | 1.03542 .02923 | 1.02506 
1,150 | 1.11234 | 1.10688 | 1.08553 | | | | 

5,000 | 1.09076 | 1.08940 | 1.07178 | 1.06016 05248 03594 .02762 | 1.02240 | 1.01877 
10,000 | 1.08287 .08270 1.06643 1.05617 | 1.04920 03391 02607 | 1.02108 01757 


| 
20,000 | 1.07594 .07667 | 1.06160 | 1.05238 04611 03196 | 02454 .01978 | 1.01640 
| | | | 


IV. xs =r/ro = 2 sin (7/10) = 0.61804 


07618 | 1.07331 | 1.05535 | 04452 .03803 | 1.02552 | 1.01991 | 1.01658 | 1.01434 
06296 | 1.06381 | 04781 | | | } 

05145 | 1.05431 .04037 | 1.03359 | 02910 | 1.01944 | 01458 | 1.01154 | 1.00942 
.04705 | 1.05044 .03734 | 03130 | 1.02717 | 1.01814 | 1.01350 | 01056 | 1.00849 
-04306 04684 | 1.03451 .02902 1.02528 .01683 | 1.01241 | 00957 | 1.00756 


; ‘ log (1 — x2) 
which vanishes at x. = 0 since lim — er =-1. 
r2350 A2 


APPENDIX 
The integral 


J @/®) FQ) + af }[1 — Blog (a/Q)]xs* — fre — HCL — x2") log (= x2) ] 


x[(1 — P) — (PO/a) xo(1 — x) 
of equation [39], required for the range 1 — @ Q) < x3 < 1 — (a@/Q) is ex- 
pressible in partial fractions as 


a&(Q)1 (“dx ok(Q) 1” 
PQ Xwd:, x2 PQ X(A+2n) 


F aE(Q) 1 tn i@ 


~ PQ w(dQ+u) 


a'{l — 2 log (a/Q)) ere eRe. 
4PQ(A + u) zs 2 ptm 


ee 1—(a/Q) 
oan +—— | des 


2PQ(A + x) b+ X2 


2 1—(a/Q) 2 2 
a 1d 1—A 1 be 1 
a. es . r = amen | ine fh ng.) a 
aa. 1 Mie aie ae jana) te 


4. 
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where PQ(A — x2)(u + x2) = a(1 — P) + PQx2 — PQx?? 


so that A= 


A+ nu) 
oa bey et 
wA+hn) 
{—>» 


(A + us) 


pint PO (2) 
He =1--~*5+0|-). 
u(r + w) a(t 2) Q 
Evaluating the integral and using equation [40], one obtains the result 


“7. \\_ & N B 27 —3N 
I; (x3) = ail got oat 


+ —S i log 5+4 + 3 ae 8] | tog ¢ —_— s) _ log «| 


a |2y—-3N_ N a) 
+H). 1 = 


x E (1 — x3) + bg (1 +5 


+ a he + 4-8] [1 (= 
6 ve 


a — 3N 


+ oq — Pp) =| toe (1 — x3) log x; — log ( 3 log 6 
+e ae = ee («/Q) | 


1 
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. ae ~ 30- nnd iS " —_ oe 


and provided P < ae < Q/a : the G are defined throughout the range. 
If P becomes higher than this value, an ‘‘incision”’ in the range of integration 
would leave the result to O(1/Q) unaltered. 

The additional O(1/Q) term in H(PQ), the value of J3(x3) for x; = 1 
— (8/Q), arises from the terms containing log [1 + u/(1 — x3)] and log [1 
+ (uQ/a)] in the second line of the above expression. Thus while log [1 — (a/AQ)] 
and log [1 — (8/AQ)] are both O(1/Q), log [1 + (uQ/a)] = — log P + O(1/Q) 
and log [1 + (uQ/8)] = log [1 + a(1 — P)/8P] + O(1/Q). 

Again the integral 
, ae [3 (2y — N) — 3N log (8/Q)] x2'— 3N x2 — 3N (1 — x2") log (1 — x2) 5. 

x[(1 — P) + (PQ/N) x2(1 — x2)] 


Zs 


of equation [44] required for the range 0 < x; < 1 — (8/Q) is expressed as 


N 2y-N_ N ,4) | ; . 7 
POE + n) (2 4 2 18 6 i: aaa 


1— oe 1 | 
- SPOR +n) = a sta 


N? co he 1 ee 1 ae ; 
~ BP ies — Tt _ | hog (1 — x2) ds 
2PQ Re ie oio nae ae 


where PQ(E — x2)(9 + x2) = M(1 — P) + PQx2 — POQx-?, 
so that similar relations hold between é and n, a now being replaced by NV. 
In place of equation [47] one now has 


Is!" (x3) = H(PQ) + + BI ae 


cn ste 35)-m gt) 
+oloeli-!54) mel 8) 
+5 son Y | oe «1 — x3) log x — log (1 -£) log 6 
+= ee 


J eed 
2PQ-2(1 — 


| — xs) /é]" — (8/80)" ] J) 
n° +? (4 
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and the term log[{l1+ N(1 — P)/@P] arises from the fact that 
log [1 + (nQ/8)] = log [1 + N(1 — P)/BP] + O(1/Q) . 
It is noticed that if P > 1, the terms log [1 — (1 — x3)/é] although required 






N/P -1 
for the range 0 < x3 < 1 — (8/Q) cease to exist for x3 < ce) a See 





eo [(1 — ane 









becomes divergent over this subinterval of the range. 


that os 








It follows then that, strictly speaking, equation [47a] cannot be applied for 
x3 = 0, that is to give the important ratio (twar — faxis)/(twait — tmean) - 
It can be shown however that the ‘‘device’’ of putting x; = 0 in the ex- 
pression [48] does in fact give this ratio to O(1/Q) . 
Calling 7,(0) the value of J3(x3) at x3 = 0, then, since 


fie a) a eth =) 49(. 7s 
X3 rs3=N(P—1)/PQ r330 X3 Q 


it is seen that the result valid for the point x; = N(P — 1)/PQ is obtained to 
O(1/Q) by substituting 7,(0) in equation [25]. 

It remains to show that ele tr = 0)/ /(twait — tmean) differs from 
(twant — te, = n(P- »/p@) /(twait — tmean) by O(1/ Q?) 

From equation [44], for x3; < 1 — (8 Q), one ioe 


_ OTs — [4(2y — N) — 4N log (8/Q)] = - 3N x3 — 3N(1 — Xs ) log (1 — xs) 














Ox; {(P ~-— 1) — ?@, IN) x3(1 - a 


ols 
so that [ a] = (}). 
so tha ms 0X3 J2,=Nn(P-1)/P@ Q 


The application of Taylor’s theorem to the function /; about the point 

= N(P — 1)/PQ proves the result. 

Again, in the final integration for the Nusselt modulus, the result to O(1/Q) 
is clearly unaffected if the subinterval 0 < x; < N(P — 1)/PQ of length 


O(1/Q) is omitted from the range of integration. 
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AN EXPERIMENTAL INVESTIGATION OF THE HEAT TRANSFER 
TO TURBULENTLY FLOWING PRESSURIZED AIR' 


By A. W. Marris? 


ABSTRACT 


Employing a counter-flow figure-of-eight heat exchanger, direct measurements 
are made of the Nusselt modulus for radial heat transfer to air pressurized up to 20 
atmospheres for Reynolds numbers up to 1.20 X 10°. For each heat transfer 
determination a simultaneous friction factor measurement is made and it is found 
that the latter is independent of heat transfer. 

Results in reasonable agreement with the momentum transfer theory are ob- 
tained for Reynolds numbers less than 0.75 X 105, provided the ratio of the 
eddy diffusivities for heat and momentum is taken as unity. For such values of 
the Reynolds number, the same value of the heat transfer coefficient was obtained 
irrespective of whether the Reynolds number was obtained by having high pres- 
sure (density) and low velocity, or high velocity and low pressure. For higher values 
of the Reynolds number, however, the value of the heat transfer coefficient 
appeared to become dependent on the over-all heat transfer rate. 


1. INTRODUCTION 


In this study, using a figure-of-eight counter-flow heat exchanger of circular 
section, an attempt is made to determine directly the heat transfer coefficient 
for Reynolds moduli up to 1.20 X 10° by the device of increasing density rather 
than speed of flow. The work was undertaken to investigate the constancy or 
otherwise of the ratio of the eddy diffusivities for heat and momentum of the 
momentum transfer analogy. As pressurized air is used as a coolant it is also 
of interest to know whether the heat transfer coefficient is in fact independent 
of the heat transfer rate. 

It has been shown by Saunders (5) that the Prandtl modulus of air is sensibly 
independent of pressure and that the device of pressurization reduces errors 
caused by local heating due to compressibility. 

Since the theory expresses the Nusselt modulus in terms of the tube friction 
factor as well as Reynolds and Prandtl moduli, the former was measured 
directly for each heat transfer determination. 

2. PRINCIPAL NOMENCLATURE AND VALUES OF CONSTANTS OF 
APPARATUS 
L = length of stainless steel tube across which heat transfer takes place 
402.9 cm. 
r internal radius of s.s. tube = 0.850 cm. 
A; = ar. 
v, = fluid velocity in s.s. tube in cm. per sec. 
Sark 
k,.,. = thermal conductivity of stainless steel tube = 0.038 c.g.s. units at 
s* C; 
b = thickness of wall of s.s. tube = 0.1025 cm. 


1 Manuscript received September 21, 1958. 
Work carried out at the School of Engineering, University College Auckland, Auckland, 


New Zealand. 
2 Present address: Physics Department, University of British Columbia, Vancouver, B.C. 
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= quantity of heat transferred to center stream per second 


Re Pr (4ar,)SAT. 

Reynolds modulus = 2v,7 ,/v;. 

kinematic viscosity of air at mid-point of s.s. tube. 
Prandtl modulus. 

molecular conductivity of air. 

mean mixed inlet temperature. 

mean mixed outlet temperature. 

T” —T’. 


= 3(T7’+ T”) = T’ + 3AT = temperature at which property values 


are taken. 

(temperature of inside wall of s.s. tube) — (mean mixed temperature 
of stream). 

(temperature of inside wall of s.s. tube at point a distance x from 
entrance) — (7’ + (x/L)AT). 


= 67,+ «. 


= temperature drop across wall of s.s. tube = Qb/ks, 5.5. 


pressure drop across orifice gauge in dynes per sq. cm. 


= pressure drop down stainless steel tube in dynes per sq. cm. 
= Fanning friction factor for s.s. tube, defined by 6g = f(L/r,) pw’. 


area of orifice in sq. cm. 

fluid velocity in orifice in cm./sec. 

area of pipe in which orifice plate is inserted. 

fluid velocity in pipe in cm./sec. 

A,/Ao = 2.6405. 

A,/A, = 2.3026. 

density of air in pipe at orifice. 

mean density of air in s.s. tube = density at temperature 7,,. 
Po/ Pr- 

over-all pressure in apparatus in dynes per sq. cm. 
(6p + 369)/TI. 


= number of normal atmospheres pressure existing at mid-point of s.s. 


tube = (II + 6p + 44¢)/76.00 = II(1 + 2)/76.00. 


Nusselt modulus = ¢ ors m 2! Re Pr i 
SéT \ s 2L 6T 


Rev/ (f/8). 

= eddy diffusivity for heat, the effective increase in thermal diffusivity 
due to turbulence. 

eddy diffusivity for momentum, the effective increase in kinematic 
viscosity due to turbulence. 

(ex /éy )Pr. 


3. THE COUNTER-FLOW HEAT EXCHANGER 


The mechanism of the counter-flow heat exchanger is shown schematically in 
Fig. 1. The cold fluid enters the long center tube from the left-hand mixing tank 
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ieee 


Cold run In a = T+AT 


Radial Heat Flow 


Tube From Whose Hot Fluid In 
Inner Surface Heat 
Transfer Is Measured 














Fic. 1. Schematic representation of counter-flow heat exchanger. 


at temperature 7’. The hot fluid enters the annulus at the right and heat is 
transferred radially inwards across the center tube wall, so that the annulus 
fluid loses heat and the fluid under test in the center tube gains heat. This latter 
fluid leaves at a mixed temperature 7” on the right, 


[1] T” = T' + AT. 


By considering the heat flow across an element of surface area dS of the center 
tube it is easy to show (3) that if (d7) p is the difference of the mean mixed tem- 
peratures of the streams at any point P then 


[2] (dT) p/ (dT) = e-*”S 


where S is the surface area of the section PQ of the tube, U is the over-all heat 
transfer coefficient defined by 
[3] dQ = UdSdT, 


dQ being the quantity of heat flowing across dS in unit time, and uy is the differ- 
ence of the reciprocals of the thermal capacities of the fluids in the two streams, 
which is assumed constant along the tube length. 

When U is assumed constant along the length of the tube it follows that if the 
thermal capacities of the two streams are equal, that is » is zero, then dT is 
constant along the tube length. Such a state of affairs is of course only realizable 
in practice when the temperature differences are infinitesimally small. 

In the ideal case considered here, the mean temperature gradients of the two 
flows are parallel lines. In the diagram, AB represents that for inner stream and 
CD that for annulus stream. The temperature gradient along the inner tube EF 
must lie between these lines. 
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Writing 67 = (temperature of inside wail of inner tube) — (mean mixed flow 
temperature of center stream fluid), then 


[4] Nu = 3 Re Pr 2. 
4. DESCRIPTION OF APPARATUS 

The apparatus, built of mild steel except for the tube across which heat 
transfer was to take place and the orifice plate for measuring the flow velocity, 
was constructed for use with a compressor which, with appropriate adjustment 
of inlet and outlet valves, would simultaneously compress the air overall and 
drive it through the figure-of-eight circuit. 

Air from the compressor, which had passed through a moisture separator and 
filter tank, was led through 1-1/64 in. internal diameter mild steel tubing; 
through a vertical 2 ft. long temperature stabilizing water jacket; horizontally 
through an over-all pressure measuring unit comprising two Bourdon gauges set 
in a turned section contained between two flanges, the gauges graduated to read 
0-100 Ib. per sq. in. and 0-300 Ib. per sq. in. respectively; through a second 
horizontal water jacket; through the orifice gauge, the latter of the type German 
Standard Orifice 1930 (2), comprising a 1/8 in. thick stainless steel plate con- 
taining a 5/8 in. diameter orifice set between flanges in the 1-1/64 in. internal 
diameter tube; and thence to the left-hand ‘‘entrance’’ mixing tank of the heat 
exchanger. 

This cylindrical tank 4-7/8 in. in length by 2-7/8 in. internal diameter contained 
first a three-bladed propeller mounted to revolve freely in the stream, and then 
equally spaced about the circumference of the tank, three thermocouple py- 
rometers. These, of length 1/2 in., 1 in., and 1-1/2 in., were 1/8 in. diameter mild 
steel tubes of 1/16 in. bore flanged at one end. The thermocouple wires had been 
threaded through the tubes and the head of the thermocouple soldered to the 
tube end to seal the 1/16 in. bore. This pyrometer could then be inserted into 
its ‘‘saddle’’, and the flange screwed down hard against a rubber ring, rendering 
the whole pressure tight. 

From the left-hand mixing tank the stream entered a section which converged 
it into the 13 ft. long internally commercially polished 3/4 in. O.D., 17.00 mm. 
(average) 1.D., Staybrite F.D.P. stainless steel tube from the inside surface of 
which the heat transfer to the center stream was to take place. The s.s. tube 
fitted into this converging unit, and in the slot encircling the tube was set a 3/4 in. 
“OQ” ring to prevent flow between the center and annulus streams. Also in this 
section was the first pressure tapping for measuring the pressure drop down the 
s.s. tube, and thence its friction factor. 

A third unit was flanged to this converging unit, which comprised the exit 
for the annulus flow. 

After it had left the s.s. tube the fluid passed through an identical unit which 
gripped the right end of the s.s. tube, took the second friction factor pressure 
tapping, and diverged the flow into the right-hand “‘exit’’ mixing tank which, 
like the “‘entry”’ mixing tank, contained a second free propeller and three more 
thermocouples arranged in the same way. 
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The flow continued horizontally for a few feet and was then led down and 
around and back into the annulus, having been heated by means of six (and 
later eight) 65-ohm coils wound on the tube. These heating coils, capable of 
being connected in series or parallel as required, were wound from resistance 
wire ribbon of 6.10 ohms per foot, and of cross-sectional dimension 0.312 in. 
< 0.0035 in. The wire was insulated from the tube with asbestos, and then the 
coils were covered with more asbestos to minimize external heat losses. The 
maximum voltage employed for the heaters was 270 volts. 

After it had left the annulus tube the stream was then led downwards away 
from the working section (disturbances caused by the sudden expansion at the 
outlet could cause serious errors if the outlet were too close to the working 
section) ; through an open-shut valve (used in pressurizing the air prior to starting 
its bodily flow through the circuit); through a second Saunders diaphragm yalve 
placed so that it could be adjusted simultaneously with the valve on the branch 
pipe near the inlet; past a needle valve for fine adjustment; and finally let out 
into the atmosphere. 

The s.s. tube wall temperature was measured by three thermocouples soldered 
to the steel at distances }L, 3L, and ?L. The temperature of the outside wall was 
thus measured, and the required inside wall temperature was then obtained by 
extrapolation, since the wall thickness, the thermal conductivity of the stainless 
steel, and the heat flowing radially inwards through the tube wall were known 


quantities, 
To facilitate the affixing and inspection of these tube wall couples the outer 
tube of the working section was made up of four sections, which were flanged at 


the ends and connected to the annulus exit and entry units and the three tube 
wall thermocouple units by means of caps, which slid over the tube and screwed 
on to the appropriate unit, thereby squeezing the flanges against hard rubber 
pads. 

The three tube wall thermocouple units also contained three legs with feet of 
lignum vitae—these were turned and collared and affixed to the steel legs with 
‘‘Heldite’’—to support the s.s. tube. 

It was essential that both the junction itself and the leads to it, which neces- 
sarily crossed the annulus stream, be shielded from the stream. For this purpose 
a tube of bakelite impregnated paper was turned. The top end of the tube, 
which was drilled axially to take the thermocouple leads, was flanged and the 
lower end collared to take a soft rubber pad shaped to fit over the s.s. tube and 
surround the junction. The pad was attached to the tube with “‘Heldite’’. The 
flange end was screwed down against a hard rubber ring so that the lower end 
squeezed the softer pad against the s.s. tube, fulfilling the fourfold purpose of 
holding the s.s. tube rigidly against its supports; of rendering the installation 
pressure tight (in fact the thermocouple was not in the pressurized region); of 
making the soft rubber close over and clutch against the thermocouple wire 
immediately above the head, thus holding the thermocouple rigid and providing 
proof against inadvertent jerks; and finally of keeping the thermojunction and 
its leads free from extraneous heating or cooling created by the impingement of 
the annulus stream. 
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In this apparatus all the pipe and tank flanges, etc., were screwed to the pipes 
and tank sections, the pipe ends abutting against shoulders. Heldite was used as 
packing, and in cases where this gave leaks the threads were sweated together 
with solder. Again the corners and T-pieces were cut from blocks of metal, into 
which the pipes were similarly screwed. 

Rings of cork impregnated paper, gaskets smeared on both faces with Heldite, 
were used as sealings at the spigots of all the flanges. 

The working section of the apparatus, including the mixing tanks, was con- 
tained in a 20 ft. long pinex-lined wooden box, complete with lid. The box was 
built up of sections fitted together by tongue and slot joints. 

The nine thermocouples, three for the determination of the center stream 
inlet temperature 7”, three for the longitudinal gradient down the stainless steel 
tube, and three for the center stream outlet temperature, were constructed from 
gauge 21 (Imperial Standard) copper and constantan wires, No joints were made 
in the wires until they reached either the cold junction or the mercury switch. 
Since the cold junction was adjacent to the middle of the working section, the 
copper and constantan leads of the first and last three couples were 8 ft. long, 
while those of the three wall couples were 5 ft. long. 

While it was necessary that all the nine cold junctions be at the same tem- 
perature to within 1/50°C., it was only necessary that the cold junction tem- 
perature, giving the temperature level in the apparatus, be known to within say 
1°. It was thus decided that a cold junction consisting of distilled water con- 
tained in a vacuum flask, which was contained in a second larger flask, would be 
adequate. 

The nine junctions were coated with ‘“‘araldyte” wire insulating material 
(made by Ciba products) and bound together so that the heads were in the 
closest possible proximity. After the wires and a mercury-in-glass thermometer 
graduated to read to 1/20°C. had been threaded through the corks and set in 
place, the cork holes were filled up with paraffin-wax. 

The switch consisted simply of 24 mercury wells, bridged in pairs as required 
by the thick-tinned copper wires. 

A Cambridge Vernier Potentiometer reading to J, uv. (about 1/400°C.) was 
employed for the measurement of the thermocouple e.m.f.’s. 

The six mixing-tank thermocouples were removable and were calibrated at 
frequent intervals throughout the experiment. A facsimile of a tube-wall thermo- 
couple was constructed with a short length of the stainless steel tube and cali- 
brated. Likewise the orifice gauge was calibrated by timing the delivery of 
250 Ib. of water, and so also the pressure gauges. 


5. FORMULAE EMPLOYED 


The Fanning friction factor for the s.s. tube defined by 
5) 
may be expressed as 


[6] 
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where C is the compressibility factor of orifice gauge defined by 
[7] C = 1 — 0,346 6p/Tl 
and x is the orifice coefficient determined as a function of Reynolds number from 


the calibration. 
Neglecting pressure head loss and temperature difference between orifice 


gauge and entry mixing tank, one has 
[8] d= (1+ feAT)(1 +2) = 1 +8 + 4adT 
and taking the appropriate value of the density variation coefficient from the 
Critical tables, 
[9] A=1+n+4+1.9 X 10-*AT. 
When the numerical values of the apparatus constants are inserted, the friction 


factor is given by 


[10] , = 841.58 ( — 0.692 ie) (1 +n+1.9 X 10°AT) < 


For the Reynolds modulus one has the result 
; M* ie as 
[11] Re’ = s| cos iss (Cx + 
m — 1 Vi Pr 
where, utilizing Boyd’s (1) result of a 4% linear increase in viscosity for 25 atm. 
pressure, 
‘ 2 Ee s a 
[12] 710 a A323 X10 'N + 9.42 X 107°T,,) 
so that finally, on substituting the requisite numerical values, 
[13] Re? = 2.2452(1 — 0.692 6p/M1) x2{1 + 
+ (1.9AT — 3.2N — 9.42T,,) K 10-*} Nép X 105. 
The Nusselt modulus is given by 
: r AT 
4 Nu = =~ Re Pr =. 
Uy 2L 67 

The Prandtl modulus of air is independent of pressure but depends on the 
temperature according to 
[14] (Pr) 7, = 0.712(1 — 3.09 X 10-4T,,). 

In the experiment, 67, required to be the difference between the inside wall 
temperature of the heat transfer tube and the mean mixed flow temperature, is 
taken as the arithmetic mean of the three values of 67, corresponding to 
x = 1L,x = 3L,x = 3L, where 67, is defined as (temperature of inside wall of 
tube) — (T’ + (x/L) AT). 

After the temperature of the outside wall of the s.s. tube was measured it was 
necessary, to correct for the temperature drop across the wall of the tube. 

One has 
[15] iT; =3f? => G, 
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where 67,’ = measured quantity, (temperature of outside wall of s.s. tube) 
— (1" + («/L) AT), 


and e, = temperature drop across the tube wall at the point x. 






An average value e of e, is given by 


ee Cee 
[16] (“7 3 Re Pr AT iL. 












and for the Nusselt modulus one has 


onl — 3.09 X 10-*T)(1 + 7.78 Re af x 10) x10" 






[17] Nu = 0.7511 Re 





where 67 refers to the outside wall temperature as measured and is the mean of 





the three values. 










RESULTS 
The values of the reciprocals of the friction factor obtained for the stainless 
steel heat transfer tube are plotted in Fig. 2. To each point there corresponds a 
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Fic. 2. Experimental values of Fanning friction factor. Continuous curve is von Karman’s 
law for a smooth tube. 





Nusselt modulus determination, see below. The continuous line is von Karman’s 










curve 
[18] 1/f = (4.0 login Q + 1.4062)? 
where Q = Rev/(f/8) 






appropriate to a smooth tube. 

No variation in f with heat transfer rate could be detected. 

Fig. 3 shows the totality of values of the product Nuf for Reynolds numbers 
less than 0.75 X 105, and the continuous lines give the theoretical values (4) 
appropriate to Prandtl moduli of 0.55 and 0.75, these calculated on the basis of 
a linear temperature gradient along the tube length. 
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Fic. 3. (a) Totality of values of Nuf obtained for Reynolds numbers less than 0.75 x 105 
and over-all pressure between | and 20 atmospheres. Continuous lines are theoretical curves 
for P = 0.55 and P = 0.75 where P = kPr and k is the ratio of the eddy diffusivities for heat 
and momentum. 

(b) Subset from Fig. 3 (a) for over-all pressure between 1.00 and 1.30 atmospheres. 

(c) Subset from Fig. 3 (a) for over-all pressure between 1.30 and 3.30 atmospheres. 

(d) Subset from Fig. 3 (a) for over-all pressure between 3.30 and 6.60 atmospheres. 

(e) Subset from Fig. 3 (a) for over-all pressure between 6.60 and 20.00 atmospheres. 
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It is seen that up to Re = 0.75 X 10°, approximately the same value of Nuf 
is obtained irrespective of whether the value of the Reynolds number is obtained 
by velocity or by pressure (density), and that if the Prandt! modulus for air is 
taken as 0.70, then the value of the ratio ez/ey is unity. 

For Reynolds numbers above 0.75 X 10° it was found that, whether the tur- 
bulence be obtained by pressure or velocity, the value of Nuf was governed by 
another parameter, AJ’, the excess of the outlet mixing tank temperature over 
the inlet mixing tank temperature, a quantity which for a given Reynolds 
number is directly proportional to the heat transfer rate. 


(i) Re®(1,10.0,10)x10° 


20 24 2 32 36 40 
aT 


Fic. 4. Nuf against A7, representing heat transfer rate, for Reynolds numbers confined 
between given limits. 
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The graphs of Fig. 4 now indicate Nuf as a function of AT (of heat transfer 
rate) for Reynolds numbers confined in each case between definite limits, the 
variation in its value through Reynolds numbers varying between these limits 
being sufficiently small not to obscure the AT variation. On each of these graphs, 
for comparison, are given the corresponding theoretical values of Nuf (4) for 
Pr = 0.55 and Pr = 0.75 appropriate to the mean value of Re for the range 
considered. 
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VIBRATIONAL TRANSITION PROBABILITIES TO HIGH QUANTUM NUMBERS 
FOR THE NITROGEN FIRST AND SECOND POSITIVE BAND SYSTEMS 


By W. R. JARMAIN AND R. W. NICHOLLS 


As part of a study of the vibrational transition probability arrays for mole- 
cular band systems of specific astrophysical interest, an analytical method of 
computing the overlap integrals of the vibrational wave functions involved 
has been developed (2, 4, 8). Depending upon the way in which the method is 
applied, reliable results may be obtained without undue tedium either for 
bands the sum of whose quantum numbers does not exceed about 6, or for any 
band one of whose quantum numbers does not exceed about 3, the other 
quantum number remaining unrestricted. These bands usually include the 
vertex of the Condon parabola of the system concerned, and are also often 
just those bands which are of astrophysical interest. 

Prior to the development of this method, however, extensive tables had been 
obtained for certain band systems by the straightforward, though extremely 
tedious method of numerical integration of accurately computed Morse wave 
functions, over a much larger range of vibrational quantum numbers than has 
been previously reported (3, 5). These more extensive data have been found 


TABLE I 
PROBABILITIES FOR THE First Positive (B*IIz— A*%Z{%) BAND 
SYSTEM OF NITROGEN 


VIBRATIONAL TRANSITION 


es tee oe he bow boa Teh ot a a ie ee 
0 349 | 32) | 190} 08s | .036 | .01, | .00; | .00: 
40. | .002 | 10: | .17% | 145 | 085 | .04 | 015 | .00, | .00, | .00; 
2 205 | 21s | 115 | 002 | 07, | 125 | 11s 07, | O45 | .02; | .005| .005 | .00s 
3 050 | 30: | .039 | 16: | .03; | .009 | 06, | .09; | .08» | .07. 04, | .02, | 01s 
4 | 00, | .13, | 27: 00, | 11; | .O8 | .00, | .O1s | .06; 08, | .08) | .055 | .03, 
5 | 02; | 210] .18, | 04, | .04, | .10, | .03; | 00, 02; | .05; | .07; | 06; 
4 | 00; | .05» | 25» | .08; | .10; | .00, 08; | .07,| 01s | .00: | .025 05 
7 | 00, | 10; | -26, | .02) | 125} .005 | .049 | 076 03. | 00, | 005 
8 | | 01s | 16; | 249 | 00) 11s | .035 | Oo | .06s | 06, | .O1; 
9 | 002 | .03s | 20; | 185 | .015 | .07; | .06. | 00» | .03; | .06, 
_ | | | 00, | .054 | .22 | .13, 055 | .03; | 08; | .015 | 005 
un | | | | | | 00, | 08 | .267 | .060 | .08 | .016 | .075 | .03. 

| | | 

12 | Ole | 116 | 26; | 030 | 10; | 00 | 05. 
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useful in the interpretation of laboratory spectra where large vibrational 
quantum numbers are often encountered. They are therefore reported below. 
The entries in Tables I and II represent p(v’, v’’) = | fy’) v(v'’) dr |? where 
¥(v’) is the wave function of the upper vibrational level v’, and y(v’’) is the 
wave function of the lower vibrational level v’’. r is the internuclear sepa- 
ration. Table III shows the basic molecular data employed. 


TABLE II 


VIBRATIONAL TRANSITION PROBABILITIES FOR THE SECOND PosITIVE (C%II,— B*Ily) BAND 
SYSTEM OF NITROGEN 



































etal Be eae sielz s || nm | 12 
0 | 445 | .325| 146 | .055 | .016| .005 | .00, | | é 
39; | 010 | .20; | .195 | .115 | 047 | .020| .00. | .00; | .00, 
2 135 | .32 032 | 062 | .15s| .142 | 082 | 039 | .016 | .005 | .002 | .00, 
3 02: | .255| .16s | 116 | .00, | .092 | 141 | .10; | .057 | 02. | .015 | .002 | 002 
f | 00: | 070 | 304 04s | .15¢| 01s] -03: | .005] 115 | .07;| .049 | 01; | .O1 
TABLE III 


MOLECULAR CONSTANTS 








(The notation 1 is s the s same as s that used in Reference 8) 
Saas = a - - = => 
State we (cm.~!) | we e (cm.~ ” “ (at. wt. units)} re (8) a (cm. ~1) | K 
CHI 2035.1 | 17. 08 7, 00377 1.1482 2. 6636 x 108. 119. 15 
B11 1734.11 14.47 7.00377 1.2123 2.4517 X 108 119.84 
A’z 1460.37 13.891 7.00377 1.293 2.4021 xX 108 105.138 


| 
| 


Table I contains results for the First Positive System (B*II — A*Z) of Nz 
over the ranges of vibrational quantum numbers 0 < v’ < 12,0 < v” < 12. 

Table II contains results for the Second Positive System (C*II — B*II) of 
Nz over the ranges of vibrational quantum numbers 0<v’ <4, 0<v’’ <12. 
These values, of course, supersede earlier provisional results (6, 7), which were 
computed under simplifying approximations upon the form of vibrational 
wave functions. 

Figs. 1 and 2 plot these numerical data diagrammatically to indicate main 
and subsidiary Condon parabolae and also to show how. the oscillator strength 
for the electronic transition is shared among the bands. The trends suggested 
by similar figures for the dependence of the overlap integrals in v’ or v”’ 
progressions upon the varying quantum number have prompted some analyt- 
ical examination of the situation. A more detailed discussion of this work will 
be reported later. 

It should be realized that the numbers in Tables I and II are accurately 
computed.values of an approximation to the ‘true’ vibrational transition 
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Fic. 1. Block diagram of the vibrational transition probability array, together with 
Condon parabolae, for the First Positive band system of nitrogen. 
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Fic. 2. Block diagram of the vibrational transition probability array, together with 
Condon parabolae, for the Second Positive band system of nitrogen. 


probabilities because the electronic transition moment R,(r) has been assumed 
independent of r in the familiar transition integral Sve’) R.(r) v(v"’) dr. 
However, it may be briefly reported here that Turner and Nicholls have 
recently used the data of Table I together with computed values of the inte- 
grals fv’) v(v"’) r dr, fur’) y(v’’) r? dr for the same bands and experiment- 
ally measured values of their intensities to determine the dependence of R,(r) 
upon r by a method suggested by Fraser (1). Further details of this work will 
be published soon. 

In the interests of brevity, references to the work of other investigators on 
comparable problems in molecular transition probabilities have not been made 
here, but will be found elsewhere in more complete discussions (2, 4, 5). 

The studies reported in this note are part of a project supported by the Air 
Force Cambridge Research Centre with Contract AF 19(122)-—470. 
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